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CHAPTER  ONE 


INTRODUCTION 


INTRODUCTION 


In  the  past.  Geotechnical  Engineering  has  always  been  sort 
of  a  "black  science"  as  opposed  to  an  "exact  science", 
where  the  ability  to  predict  what  will  occur  beneath  the 
earth's  surface  was  often  based  priaarily  on  experience. 
Over  the  years  this  has  slowly  changed.  With  the  advent 
of  the  Finite  Eleaent  aethod  of  analysis,  Geotechnical 
Engineering  leaped  forward  towards  becoaing  a  aore  exact 
science.  The  ability  to  predict  is  still  based  on 
experience.  However  the  predictions  are  auch  better 
today . 

Finite  eleaent  analysis  has  a  couple  of  draw  backs  though. 
One  draw  back  is  that  aesh  generation  is  very  tedious  and 
tiae  consuaing.  The  second  draw  back  is  that  the  overall 
aatheaatical  analysis  is  also  very  tine  consuaing.  With 
the  invention  of  the  conputer  the  aatheaatical  analysis 
draw  back  is  pretty  such  avoided,  although  at  a 
substantial  nonetary  cost.  Perfornlng  a  finite  eleaent 
analysis  on  a  nainfrane  conputer  takes  a  lot  of  tine  and 
as  a  result  cost  a  lot  of  noney.  As  coaputers  have  become 
aore  sophisticated  and  quicker,  the  running  tiae  has 
decreased.  However  the  cost  per  unit  of  tine  has 
increased.  Therefore,  a  dynanic  finite  eleaent  analysis 
still  costs  a  fair  anount. 
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With  the  fast  coaputers  now  available  to  perfora  finite 
eleaent  analyses,  they  becaae  auch  aore  coaaon.  However 


they  still  have  the  first  drawback  of  the  tedious  aesh 
generation.  That  is  until  prograas  were  developed  that 
perfora  aesh  generation  theaselves. 

Coaputers  have  continued  to  becoae  aore  sophisticated. 

With  the  inception  of  the  386  generation  personal  coaputer 
(PC),  aany  of  thea  have  a  greater  capability  than  the 
aainfraae  coaputers  of  yesteryear.  In  addition  to  their 
large  capabilities  these  PCs  are  relatively  inexpensive. 

With  the  increasing  capabilities  of  these  PCs  there  there 
is  an  increasing  aaount  of  work  that  can  be  accoaplished 
with  thea.  As  a  result  over  the  last  forty  years  the 
engineering  field  has  gone  froa  the  very  tedious  hand 
developaent  of  a  finite  eleaent  aesh  and  its*  analysis,  to 
just  recently  being  able  to  use  a  PC  to  develop  and 
analyze  a  finite  eleaent  aesh. 

This  capability  of  using  a  PC  to  develop  and  analyze  a 
finite  eleaent  aesh  has  opened  a  whole  new  area  of 
analyzing  coaplex  probleas  by  the  practitioner.  Not  only 
has  it  opened  new  areas  but  it  has  substantially  decreased 
the  cost.  For  the  cost  of  perforaing  three  fairly  coaplex 
analyses  on  a  aainfraae  coaputer,  a  consulting  fira  can 
purchase  a  PC  systea  capable  of  perforaing  those  saae 
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analyses,  a  coaHercial  finite  element  analysis  program 
with  mesh  generation  and  post  processing  capability,  and 
have  personnel  trained.  As  a  result,  the  actual  cost  to 
the  consultant  for  performing  a  fourth  analysis  is  the 
personnel  cost  of  developing  the  mesh,  which  would  have 
occurred  any  way,  and  the  minimal  cost  of  computer  and 
program  maintenance. 


There  are  many  comprehensive  finite  element  analysis 
programs,  for  dynamic  and  static  problems,  available  on 
the  market  today.  Programs  such  as  STRUDL ,  NISA,  ANSYS , 
and  NASTRAN  are  readily  available,  to  Just  name  a  few,  and 
relatively  easy  to  use.  After  obtaining  an  office 
computer  system  and  finite  element  program  package,  a  new 
user  merely  has  to  sit  down  and  "use"  the  program  to  learn 
how  it  operates.  This  "use"  should  consist  of  first 
reading  a  few  of  the  basic  chapters  in  the  user's  manuals 
and  then  trying  to  develop  and  analyze  problems  with 
"known"  solutions.  Once  the  new  user  is  able  to  obtain 
analysis  responses  that  natch  the  "known"  solutions  they 
are  ready  to  go.  As  becoming  intimately  familiar  with  a 
comprehensive  computer  program  may  take  many  months,  it  is 
)  recommended  that  the  new  user  "cross-check"  the  program 

output,  with  another  method,  until  they  are  very 
comfortable  with  the  program  and  it»  performance. 


) 


Until  the  start  of  this  master's  research  problem  the 
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author  had  very  little  computer  experience  and  only  a  very 


basic  understanding  of  what  a  finite  element  analysis  was. 
This  research  problem  consisted  of  providing  a  "new  user" 
the  opportunity  of  learning,  from  scratch,  how  to  generate 
finite  element  meshes  and  analyze  them  using  the  NISA  II 
computer  program.  This  involved  loading  the  program  onto 
a  PC,  modeling  problems  that  had  solutions  readily 
available  (ie  verification  problems  obtained  with  the 
user's  manuals)  and  then  analyzing  problems  and  comparing 
the  response  to  theoretical  responses  or  those  provided  by 
a  different  finite  element  program  (GTSTRUDL). 

Tl»«-~foJLl awing-  oliapters  wi-11-  provide  an  individual  that  has 
no  finite  element  experience  and/or  NISA  II  program 
experience  with  a  more  clarified  user's  manual, 
highlighting  the  problem  areas  the  author  encountered,  and 
some  comparisons  between  NISA  answers  and  those  provided 
by  theory  or  another  finite  element  program  (GTSTRUDL). 

Appendix  A  provides  a  list  of  the  equipment  required  to 
operate  the  NISA  II  program  and  appendix  B  provides  an 
idea  of  the  time  required  to  perform  analyses  of  various 
sizes. 
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CHAPTER  TWO 
USER'S  MANUAL 


I 
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INTRODUCTION 


There  are  basically  three  sain  methods  of  developing  a 
finite  element  model  for  analysis.  Method  nne  consists  of 
using  a  finite  element  modeling  program  like  PATRAN .  Once 
the  model  has  been  created  NISA  II  can  preform  the 
analysis  after  the  model  has  been  Interfaced  with  the 
DISPLAY-PRE  program  via  a  neutral  file  system.  Method  two 
consists  of  using  a  text  or  line  editor,  like  EDLIN,  and 
writing  an  ASCII  file  following  the  format  In  Sections  5 
to  8  (Input  &  executive  commands,  model  data,  analysis 
data,  and  modal  dynamic  data  respectively)  of  the  NISA  II 
User's  Manual.  The  third  method  consists  of  using  the 
DISPLAY-PRE  program  portion  of  the  NISA  II  package. 

This  user's  manual  will  concentrate  on  how  to  develop 
finite  element  models  using  the  DISPLAY-PRE  program,  and 
then  how  to  perform  a  finite  element  analysis  using  the 
NISA  II  program.  Since  writing  an  all  Inclusive  users 
manual  Is  beyond  the  scope  of  this  project.  It  Is  the 
author's  intent  that  this  users  manual  clarify  and 
simplify  the  generation  and  analysis  of  finite  element 
models  using  the  NISA  II  documentation.  Therefore,  this 
chapter  should  be  used  concurrently  with  the  actual  NISA 
II  User  Manuals. 

The  DISPLAY-PRE  program  is  a  user  friendly  program  that, 


in  addition  to  providing  two  nodes  for  creating  nodels 
(■enu  node  and  coaaand  node),  autonatical ly  produces  on 
the  screen  for  viewing,  a  geoaetric  representation  of  the 
aodel  being  created.  The  first  node  of  creating  a  aodel 
is  the  aenu  aode  which  guides  the  user  through  the  various 
coaaands  by  way  of  a  pyraaid.  This  allows  the  user  to 
select  a  "aain"  coaaand  and  then  allows  the  selection  of 
an  option  froa  the  "aain"  coaaand's  sub-aenu.  The  second 
aode  is  the  coaaand  node.  This  allows  the  user  to  bypass 
the  aenu  aode  and  input  desired  coaaands  and  sub-connands 
directly  (for  soneone  who  is  faniliar  with  the  progran  the 
coaaand  node  is  the  quicker  nethod).  The  coaaands  for 
both  nodes  are  identical.  To  access  the  coaaand  node  once 
the  user  is  in  the  DISPLAY'PRE  progran  nerely  type  "c". 


COMMAND  FORMAT 

The  general  foraat  for  the  coaaand  input  consists  of: 

Entry  •  1  2  345 

variable  naae  function  input  output  data 

Nane:  This  is  the  "aain"  coaaand. 

Function:  This  is  the  desired  option  available  froa  the 

"aain"  coaaand. 

Input:  This  is  the  list  of  identities  that  the  function 
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uses  to  perform  the  command. 


Output:  This  is  the  list  of  generated  entitles. 

Data:  This  is  the  data  used  by  the  command  to  perform  the 

required  function. 


TYPICAL  COMMANDS 

The  following  commands  are  available  in  the  DISPLAY-PRE 
program  and  are  typical  of  those  used  to  create  a  finite 
element  model  for  use  in  Geotechnical  Engineering.  For 
the  complete  list  of  all  available  commands  refer  to  the 
DISPLAY  II  User's  Manual. 

Name  Description 

GRD:  This  stands  for  grid  which  is  defined  as  a 

point  in  3-0  space.  The  user  can  create,  add, 
delete,  plot,  search,  etc.  a  grid  using  this 
command.  Refer  to  pages  3.1-1  to  3.1-37  of  the 
DISPLAY  II  User's  Manual. 

LIN:  This  command  is  used  to  define  a  line.  Lines 

can  be  created  through  two  or  more  grids, 
plotted,  deleted,  etc.  using  this  command. 

Refer  to  pages  3.2-1  to  3.2-64  of  the  DISPLAY- 
PRE  User's  Manual. 
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PAT:  This  cosaand  stands  for  patch  and  is  used  to 

create  or  aanipulate  a  patch.  A  patch  is  a 
surface  defined  by  paraaetric  cubic  equations. 
Refer  to  pages  3.3-1  to  3.3-46  of  the  DISPLAY 
II  User's  Manual.  Refer  to  the  eigenvalue 
exaaple  to  see  a  aethod  of  creating  a  patch  and 
what  one  looks  like. 

HYP:  This  coaaand  stands  for  hyperpatch  and  is 

siailar  to  the  PAT  coaaand  only  it  is  used  to 
create  a  3-D  solid  entity.  Refer  to  pages  3.4- 
1  to  3.4-28  of  DISPLAY  II  User's  Manual. 

PEG:  This  coaaand  stands  for  finite  eleaent 

generation  and  is  used  to  create  the  finite 
eleaent  aesh.  Refer  to  pages  4.1-1  to  4.1-28 
of  the  DISPLAY  II  User's  Manual. 

FZM:  This  coaaand  stands  for  finite  eleaent 

generation  with  "topological  zooa" .  With  this 
coaaand  the  user  can  create  local  high  aesh 
densities  without  disturbing  the  rest  of  the 
aesh.  Refer  to  pages  4.2-1  to  4.2-36  of  the 
DISPLAY  II  User's  Manual. 

NOD:  This  coaaand  stands  for  node  and  can  be  used  to 


add,  delete,  search  for,  plot,  etc.  a  node  or 
list  of  nodes.  Refer  to  pa^es  4.3-1  to  4.3-30 
of  the  DISPLAY  II  User's  Manual. 


.E:  This  coaaand  stands  for  element  and  is  similar 

to  the  node  command.  Refer  to  pages  4.4-1  to 
4.4-60  of  the  DISPLAY  II  User's  Manual. 


0: 


This  command  stands  for  property  and  is  used  to 
specify  the  geometric  properties  of  a  property 
index.  For  ease  of  providing  the  specified 
properties  every  element  in  the  model  is 
assigned  a  property  index  number.  This  allows 
the  user  to  input  one  command  to  assign  the 
geometric  properties  to  a  group  of  elements 
Instead  of  Inputing  one  command  for  each 
element.  Refer  to  pages  4.5-1  to  4.5-8  of  the 
DISPLAY  II  User's  Manual. 


iT :  This  command  stands  for  material  and  is  used  to 

specify  the  actual  material  properties  (example 
Poisson's  ratio)  to  a  material  index.  Like  the 
property  index  each  element  is  assigned  a 
material  index.  Therefore  the  user  can  assign 
a  material  property  to  a  group  of  elements  with 
one  command.  Refer  to  pages  4.6-1  to  4.6-8  of 
the  DISPLAY  II  User's  Manual. 
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The  following  example  is  provided  to  show  how 
the  property  index  and  the  material  index  is 
used  to  provide  properties  to  various  items. 

The  item  to  be  modeled  consists  of  the 
following  materials:  a  concrete  beam  18  inches 
thick  and  E»4.074.280  psi  (f'c=5000  psi),  a 
concrete  beam  18  inches  thick  and  £>3,155,924 
psi  (f'c>3000  psi),  a  concrete  slab  (6  inches 
thick  and  £>4,074,280  psi  (f'c>5000  psi),  and  a 
concrete  slab  24  inches  thick  and  £>3,155,924 


psi  (f'c>3000 

pal)  . 

In  this  case 

the  material 

and  property 

index 

numbers  could 

be  assigned  as 

follows : 

Material 

£  (psi)  T 

(In) 

MID  * 

PID  ♦ 

beam 

4,074,280 

18 

1 

1 

beam 

3 , 155,924 

18 

2 

1 

slab 

4,074 , 280 

6 

1 

2 

slab 

3,155,924 

24 

2 

3 

FOR:  This  stands  for  force  and  allows  the  user  to 

manipulate  the  forces  in  the  model.  Note; 
Forces  are  only  applied  to  nodes.  Refer  to 
pages  4.7-1  to  4.7-14  of  the  DISPLAY  II  User's 
Manual  . 

PR£:  This  stands  for  pressure  and  allows  the  user  to 

manipulate  the  pressures  in  the  model.  Note: 
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Pressures  are  only  applied  to  elements.  Refer 
to  paees  4.8-1  to  4.8-10  of  the  DISPLAY  II 
User's  Manual. 

DIS:  This  stands  for  displacement  and  allows  the 

user  to  manipulate  displacement  constraints,  le 
boundary  conditions,  at  specified  nodes. 
Although  the  user's  manual  does  not  state  It 
the  author  determined  that  the  program  will 
accept  a  maximum  of  ten  variables  between 
slashes  "/“  In  entry  3  of  this  command.  If 
greater  than  ten  variables  are  required  to 
define  the  displacement  constraints,  multiple 
DISP  commands  are  required.  Refer  to  the 
various  verification  examples  and  pages  4.9-1 
to  4.9-14  of  the  DISPLAY  II  User's  Manual. 


2 . 7 


GENERAL  EXECUTIVE  COMMANDS 

DEC:  This  conmand  places  the  prosraa  in  the  mode  to 

accept  Executive  coamands.  Executive  coaaands 
tell  the  prosraa  what  type  of  analysis  to 
perfora,  what  files  to  save  for  post¬ 
processing,  etc.  Refer  to  page  5.1-1  of  the 
NISA  II  User's  Manual.  Appendix  C  is  a 
photocopy  of  table  5.2.3  on  page  5.2-7  of  the 
NISA  II  User's  Manual  and  provides  an 
alphabetical  listing  of  all  executive  coaaands 
available. 

I  INS:  This  is  an  executive  coaaand  that  allows  the 

insertion  of  data  lines. 

I  ANAL>:  This  tells  the  prograa  what  type  of  analysis 

to  perfora.  For  exaaple  "ANAL-STATIC"  tells 
the  program  to  perfora  a  static  analysis. 

FILE-:  This  tells  the  prograa  the  prefix  to  be  used 

when  saving  the  desired  data  files.  For 
exaaple  "FILE-SOl"  tells  the  prograa  that  the 
prefix  SOI  is  to  be  used  when  saving  any  binary 
data  file.  The  author  determined  that  using 

I  greater  than  five  alphanumeric  characters  in 

this  coaaand  will  cause  problems  with  the 
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prograa  trying  to  retrieve  the  internally 
generated  data  files  used  for  the  dynamic 
analysis  and  post  processing. 

Note  this  command  is  to  be  used  in  conjunction 
with  the  following  SAVE-  command. 

SAVE*:  This  tells  the  program  what  data  files  are  to 

be  saved.  For  example  "SAVE-26,27"  tells  the 
program  to  save  the  internally  generated  data 
files  26  and  27.  These  files  are  named 
S0126,DAT  and  S0127.DAT  respectively  and  are 
required  if  the  user  desires  to  perform  any 
post-processing.  Section  3  of  the  NISA  II 
User's  Manual,  especially  pages  3.15-1  and 
3.15-2,  describes  a  number  of  internally 
generated  data  files  the  program  uses  to 
perform  its  analysis.  All  of  these  files  are 
temporary  files  and  will  not  be  saved  unless, 
specifically  saved  with  this  command. 

QUl:  This  command  either  takes  the  prograa  out  of 

the  INS  mode  or  the  DEC  mode  depending  on  which 
mode  the  prograa  is  in  when  the  QUI  command  is 
input . 
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COMMAND  RULES 


The  general  rules  governing  coaaand  input  are: 

Coaaands  can  either  be  typed  in  capital  letters  or 
lower  case  letters. 

Coanas  separate  the  variables  (ie  name,  function, 

input,  output,  data).  Two  consecutive  coaaas  indicate  a 
null  (void)  entry.  The  nroeraa  will  autoaatical Iv  provide 
"DEFAULT"  values  to  null  entries,  if  applicable. 

Slashes  "/"  separate  aultlple  entries  in  one  variable 
(exaaple:  naae,  function,  input  fl/input  #2/input  *3, 
etc.).  Like  coaaas,  two  consecutive  slashes  indicates  a 
null  (void)  entry.  Again  the  nrograa  will  autoaatical Iv 
provide  "DEFAULT"  values  to  null  entries,  if  applicable. 

Units:  Any  unit  (exaaple  feet,  foot-pound,  aeter, 

etc.)  can  be  used  in  this  prograa.  The  only  requlreaent 
is  that  the  user  aaintalns  unit  consistency. 

For  a  aore  in  depth  description  of  the  general  rules 
refer  to  page  2.1-3  of  the  DISPLAY  II  User's  Manual. 


2.10 


PROGRAM  USAGE 


MODEL  GENERATION;  The  follotrlng  steps  are  required  to 
generate  a  finite  eleaent  aodel: 

Once  the  machine  is  on.  get  into  the  directory  and 
sub-directory  that  the  appropriate  files  will  be  stored 
in.  For  example  on  the  machine  the  author  used,  all  model 
generation  was  accomplished  in  the  sub-directory  "TEST" 
under  the  directory  "EMRCNISA"  on  the  C  drive.  Therefore, 
to  create  a  model  the  author  typed  in  "cd  emrcnisaXtest" 
once  the  C>  was  received.  This  changed  the  directory  from 
the  C  directory  to  the  sub-directory  "TEST"  of  the 
directory  "EMRCNISA". 

Now  type  "NISA286"  and  hit  the  return.  This  will 
automatically  access  the  main  menu  of  the  NISA  II  program, 
with  the  "DISPGDB"  option  highlighted.  Hit  return  again 
to  access  the  DISPLAY-PRE  program. 

The  program  is  now  ready  to  accept  the  commands 
required  to  generate  the  finite  element  model.  As 
discussed  above  this  can  be  accomplished  either  in  the 
user  friendly  menu  mode  or  the  command  mode. 

To  end  the  DISPLAY-PRE  session  type  in  the  command 
"END" . 

RUN  ANALYSIS;  To  perform  an  analysis  the  following  steps 
are  required: 

Access  the  directory  and  sub-directory  that  the  file 
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containing  the  model  to  be  analyzed  is  in,  then  type 
"NISA286"  and  hit  return. 

Now  move  the  cursor,  using  the  arrow  keys,  to  the 
desired  analysis  (Example:  STATIC/EIGN/BUCK  for  either  a 
static,  eigenvalue  or  buckling  analysis)  mode  and  hit 
return.  Note  if  a  dynamic  analysis  Is  desired,  an 
eigenvalue  analysis  must  be  performed  first  so  that  the 
natural  frequencies  of  the  model  can  be  used  in  the  actual 
dynamic  analysis  (for  more  info  refer  to  the  dynamic 
example ) . 

No  matter  what  type  of  analysis  is  to  be  performed 
the  program  will  ask  for  the  input  file  name  at  this  time. 
For  a  static,  eigenvalue,  or  buckling  analysis  type  the 
file  name  generated  with  the  "WRI , NISA , ( f i le  naae).NIS" 
command  and  hit  return.  For  a  dynamic  analysis  type  in 
the  file  name  that  was  created  using  a  text  editor. 
Remember  as  stated  above  an  eigenvalue  analysis  must  be 
accomplished  before  a  dynamic  analysis  can  be  completed. 

Now  the  program  will  ask  for  the  output  file  name. 
This  is  the  file  name  the  program  will  store  the  analysis 
in .  For  ease  of  file  maintenance,  it  is  recommended  to 
use  the  same  prefix  as  the  input  name  and  use  the  .out 
suffix  (ie  "file  name. OUT").  Note:  if  an  output  file  name 
is  used  that  is  currently  employed  in  the  originating 
sub-directory,  the  file  in  that  sub-directory  will  be 
written  over.  After  Inputing  the  output  file  name  hit  the 
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return  and  wait  for  the  analysis. 

To  end  an  analysis  prior  to  its  conpletion  hit  the 
control  key  "CTRL"  and  the  "Break"  key  s laul taneous ly . 

The  program  will  then  respond  with  "Strike  any  key  when 
ready".  At  this  tine  strike  any  key  and  the  program  will 
access  the  "main"  menu.  Now  move  the  cursor,  using  the 
arrow  keys,  to  "QUIT"  and  hit  return.  The  program  will 
now  return  back  to  the  originating  directory  and 
sub-directory . 

EXIT  PROGRAM:  To  exit  any  portion  of  the  DISPGDB  (DISPLAY 
II)  program,  type  end.  The  program  will  then  ask  if  a 
data  base  file  is  to  be  written.  Respond  with  either  yes 
"Y"  or  no  "N"  and  hit  return.  The  program  will  now  access 
the  "main"  menu.  Now  move  the  cursor,  using  the  arrow 
keys,  to  "QUIT"  and  hit  the  return.  The  program  will  now 
return  back  to  the  originating  sub-directory. 

Refer  to  "RUN  ANALYSIS"  above  for  the  method  of 
exiting  the  program  prior  to  the  completion  of  an  ongoing 
analysis . 

VIEW  A  PILE:  To  view  any  file  first  access  the  directory 
and  sub-directory  that  the  file  Is  in.  Then  use  standard 
DOS  commands  such  as  "type"  with  the  file  name  placed 
after  the  word  "type"  once  a  C>,  D>,  A> ,  etc.  for  the 
proper  directory  and  sub-directory  has  been  received. 
Example  Otype  sOl.std;  this  will  access  the  file  sOl.std 


In  the  C  drive  and  scroll  It  across  the  screen.  To  stop 
the  screen  fro«  scrolling  hit  the  control  key,  "CTRL"  and 
the  s  key  simultaneously.  To  start  the  scrolling  again, 
hit  any  key. 


PRINT  A  FILE:  To  print  a  file  add  a  ”>prn”  after  the 
above  discussed  type  command.  Example  Otype  s01.8td>prn 
The  fill  sOl.std  in  the  C  drive  will  now  be  printed. 
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HELPFUL  HINTS 


Based  on  the  author's  experience,  because  the  NISA  II  and 
DISPLAY  II  programs  are  so  user  friendly,  there  really  is 
a  limited  number  of  helpful  hints  that  can  be  provided  to 
the  unfamiliar  user.  Most  of  the  helpful  hints  are 
methods  of  avoiding  problems  when  developing  either  the 
DISPLAY  II  file  or  the  ASCII  file,  required  for  dynamic 
analysis.  Wherever  a  problem  or  a  limitation,  not  listed 
in  the  NISA  II  or  DISPLAY  II  User's  Manual,  was  discovered 
a  method  of  avoiding  this  problem  is  provided.  These 
methods  are  listed  wherever  they  are  appropriate,  whether 
it  be  in  the  general  command,  Executive  command,  or 
verification  problem  section  of  chapter  two. 

Additionally,  the  following  hints  are  provided: 

1.  Whether  the  user  is  experienced  or  a  novice  it  is 
recommended  that  the  user  sit  down  and  develop  the  model, 
to  be  analyzed,  with  paper  and  pencil  prior  to  sitting  at 
the  key  board.  The  main  reason  is.  it  is  much  easier  to 
erase  an  unwanted  or  incorrect  data  line  on  paper  than  it 
is  to  delete  one  once  it  has  been  Inputed  into  the  DISPLAY 
II  program.  In  fact,  some  command  lines  can  not  be 
deleted.  As  a  result  if  one  of  these  lines  is  incorrectly 
Inputed.  the  user  must  start  completely  over  inputing  the 
required  model  generation  commands  to  correct  the  model. 

2.  Keep  the  hard  disk  free  of  unwanted  clutter  and 
unnecessary  files.  Due  to  the  program  storing  numerous 


2.15 


teaporary  files  when  In  the  process  of  perforaing  an 
analysis  a  10,000  degree  of  freedoa  aodel  requires 
approxlaately  12  aegabytes  of  disk  space.  A  10,000  DOF 
systea  aay  seea  excessive  for  a  Geotechnical  Engineer  and 
therefore  why  worry  about  apace?  The  reason  this  "house 
keeping"  needs  to  be  accoaplished  is  because  a  1350  DOF 
systea,  which  is  easily  within  the  range  a  Geotechnical 
Engineer  would  use,  requires  approxlaately  2  aegabytes  of 
disk  space.  If  there  is  not  enough  disk  space  available 
two  things  will  happen:  1.  the  prograa  will  not  coaplete 
the  analysis  and  2.  because  the  prograa  does  not  coaplete 
its*  analysis  all  of  the  teaporary  files  created  reaaln  on 
the  disk  until  deleted.  The  author  at  one  point  atteapted 
a  1350  DOF  aodel  with  1.8  aegabytes  of  space  available  and 
ended  up  with  the  prograa  "crashing",  0  bytes  of  disk 
space  available,  and  approxlaately  20  teaporary  files  on 
the  disk  that  had  to  be  deleted.  When  the  systea 
"crashed"  these  teaporary  files  are  stored  with  the  prefix 
the  user  Inputed  with  the  Executive  coaaand  "File-"  and 
the  suffix  ".tap".  For  exaaple  if  the  Executive  coaaand 
"File-sOl"  was  used  and  the  systea  "crashed"  one  of  the 
teaporary  files  saved  on  the  disk  aay  be  "s0128.tap". 

Note:  all  the  internally  generated  files  are  nuabered 
26 . ,  27 . ,  28 . ,  32 . ,  etc . 
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EXAMPLE  PROBLEMS 


Because  the  NISA  II  Verification  Probleas  Manual  was 
written  prior  to  the  developaent  of  the  DISPLAY-PRE 
prograa.  those  probleas  do  not  provide  the  DISPLAY-PRE 
coaaands  required  for  aodel  eeneration.  As  a  result  the 
following  probleas  were  selected  froa  the  Verification 
Probleas  Manual,  as  good  representative  exaaples  of  what  a 
Geotechnical  Engineer  aay  encounter,  for  a  step  by  step 
description  of  how  to  develop  the  model  using  DISPLAY-PRE 
and  then  how  to  perfora  the  desired  analysis. 

In  order  to  coaplete  a  dynaaic  analysis  of  a 
structure,  facility,  aass  etc.  two  files  are  required  to 
"design"  the  finite  eleaent  aodel.  The  first  file 
contains  the  actual  geoaetric  and  aaterial  properties  of 
the  aodel  and  is  developed  using  the  DISPLAY  II  prograa. 
This  file  is  then  used  to  coaplete  an  eigenvalue  analysis 
of  the  aodel.  The  eigenvalue  analysis  deteraines  the 
natural  frequencies  of  the  systea,  which  is  then  used  in 
the  dynaaic  analysis.  The  second  file  is  written  using  a 
standard  text  or  line  editor,  like  EDLIN,  and  contains  the 
loading  characteristics  of  the  aodel.  Additionally,  this 
file  contains  the  coaaands  required  to  tell  the  prograa 
what  the  desired  output  is.  This  file  is  generated 
following  Section  8  "MODAL  DYNAMIC  DATA"  of  the  NISA  II 
User's  Manual. 

The  following  exaaples  consists  of  the  input  file 
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(the  actual  coaaands  Inputed  during  the  DISPLAY-PRE 
session) ,  the  NISA  file  (the  actual  file  that  is  created 
by  the  DISPLAY-PRE  pre-processing  prograa  and  is  used 
during  the  finite  element  analysis),  and  a  line  by  line 
description  of  the  input  file.  Additionally,  the  dynamic 
problems  contain  the  dynamic  input  file  created  using  the 
text  or  line  editor  and  a  line  by  line  description  of  this 
input  file. 

The  line  by  line  descriptions  discussed  above  do  not 
discuss  every  line  of  the  input  file.  However,  they  do 
discuss  every  different  line  of  input.  If  there  is  more 
than  one  line  of  input  performing  the  same  function,  for 
example  a  model  may  contain  five  lines  of  input  listing 
the  movement  constraints  the  model  will  have,  only  one  of 
the  input  lines  is  discussed. 
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NISA  Verincaiion  Manual 


VERIFICATION  PROBLEM  2.20  ( 

TITLE: 

Beam  on  a  tensionless  elastic  foundation 
ELEMENT  TYPE: 

2-D  beam  element  (NKTP  «  13.  NORDR  »  1) 

2-D  spring  element  (NKTP  =  18,  NORDR  »  1) 

2*D  gap  element  (NKTP  ®  42) 

PROBLEM: 

A  30  inch  long  beam  (3  inches  wide,  1  inch  deep)  rests  on  an  elastic  soil  foundation  (Fig.  2.20.1).  The 
soil  is  active  in  compression  only.  A  concentrated  load  of  80(X)  lbs.  is  applied  at  the  midspan  of  the 
beam.  The  foundation  modulus  is  5000  IbM^ 

The  midspan  displacement  and  the  contact  length  of  the  beam  are  calculated. 

PROPERTIES: 

Material: 

EX  m  30XE6ptl  (Modulus  of  dasddty  of  bent  maiefial) 

NUXY  «  0.0  (Pdsaon'siatio) 

CrasS'Secdon: 

IZZ  m  025  in*  (Secood  area  moment  about  beam  Z>axis) 

A  «  3  in*  (Cross-aecdcnal  aica) 

FINITE  ELEMENT  MODEL: 

Due  to  symmetry,  only  one  half  of  die  structure  is  modeled  uring  2>D  beam,  2-D  spring  and  2-D  gap  (or 
compiessioiwmly  spring)  elements  (Hg.  2.20.1).  The  beam  is  divided  into  cigto  equal  length  beam 
elements  with  two  nodes  per  elemeoL  Ilie  tenskmkss  ftmndstiOD  is  simulated  bjr  nine  (Ssoeie  equkBstnt 
springs.  The  two  springs  (near  die  midspan)  are  assumed  to  be  active  in  tendon  and  compiesdon.  These 
two  firings  support  die  beam  initially.  The  remaining  seven  compresdononly  firings  are  modeled  by 
gap  elements  with  negative  flexiUlity  coefficients  and  zero  gap  wid& 

The  spring  constant  (K)  is  taken  as  foundadon  modulus  dmes  beam-width  times  spacing  between  the 
springs  (K  «  5000*3*1.875  «  28125  IbAn).  The  flexibility  coeffident  (F)  for  a  typical  compression-only 
spring  is  the  reciprocal  of  the  spring  constant  (F  « -l/K). 
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BRD,ADD, 1 ,0/0/0 
tiRD,TRS, 1,2,15 
L1N,2GD, 1,2,1 
PEG, BAR, 1 ,,,8/13,1/1 
LIN.TRS, 1 ,2,/-10 
PEG, BAR, 2, , ,8/13 
ELE,DEL,9Ti6 
ELE, ADD, 9, 18/1 //2, 1/10 
ELE,ADD, 10, 18/l//3,2/l 1 
ELE, ADD, 1 1 ,42/l//4,3/12 
ELE,ADD, 12,42/l//4,4/13 
ELE,ADD, 13,42/l//4,5/14 
ELE,ADD, 14,42/l//4,6/15 
ELE,ADD, 15,42/l//4,7/16 
ELE,ADD, 16,42/l//4,8/17 
ELE,ADD, 17,42/l//5,9/18 
D ISP, ADD, 1 , 1 ,0/////0 
Dr3P,ADD,2T9, 1 ,0 
D1BP,ADD, 10T18, 1 ,0/0 
FORCE ,  ADD  ,1,1,  /•--4000. 

PROP , ADD , 1 , 3/0 , 25/0/0 
MAT , ADD , 1 , 30 . OEA/0 . O 
PROP, ADD, 2, 14062.5 
PROP, ADD, 3, 28 125 

PROP , ADD ,4,0. /-O. 35555E-4/0. / 1  - 0E-8/-90 
PROP , ADD , 5 , 0 . / -0 . 7 1 1 1 1 E-4 /O . / 1 . OE-8/ -90 
DEC 
INS 

ANAL=STATIC 
FILE«S01 
SAVE=26 , 27 

*T I TLE, ATTEMPT  OF  VERIFICATION  PROBLEM  S20.STD 

*13,1,0,1,0,0,0, 0 ,0,0,0. OOOOE-0 1,0. OOOOE-0 1 

*14,1, POINT  LOAD  AT  BEAM'S  CENTER 

QUI 

QUI 

WRI ,NISA,S01.NIS 
WRI , DBS, SOI . BIN 
END 
N 
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**  EXKC.UriVE 
ANAl  =^-!iTATIC 

FILE-^BOl 

BAVE-=26,27 

*TIT1  F 

ATTEMPT  OF  VERIFICATION  PROBLEM  S20.STD 
♦ELTYPE 

1,  13,  1 

2,  18,  1 

3,  42,  1 

♦RCTABLE 

1,  4 

3 . 0000 ,  0 . 2500 ,  O . 0000 ,  O . 0000 

2,  1 
14062.5 

3.  1 

28125 

4,  5 

0 . 0000  ,  -  .  35555E--4 , 0 .  ,  1 .  E-8 ,  ~90 

r-  pr; 

0 .  0000 ,  - .  71111  E  -'4 , 0 .  ,  1 .  E  -S ,  --90 
■^ELEMENT 


1  ^  » 

U  .  i, 

■0 . 0000 , 

0  •  ^ 

1 . 3750 , 

'  • .  0000 , 

0  •  0C>^>C»  f 

■  T 

1  ?  T  T 

3. 7500 , 

0 . 0000 , 

0 . 0000 , 

4 

1  »  1  1 

5 . 625u , 

0.0000, 

0 .  C>00C.» , 

t.j 

?  ?  f  9 

7 . 5000 , 

0.0000 , 

0 . 0C>00 , 

o 

1  1  ?  » 

7.3750, 

0 . 0000 , 

0 . 0000 , 

7 

1111 

1 1 . 2500 , 

0 . 0000 , 

0 . 0000 , 

Q 

1111 

13. 1250, 

0.0000, 

0 . 0000 , 

9 

1111 

1 5 . 0000 , 

0. 0000 , 

0 . 0000 , 

10 

1111 

0 . 0000 , 

-10.0000, 

0 . 0000 , 

11 

1111 

1 . 8750 , 

-10.0000, 

0.0000, 

12 

fill 

3.7500, 

-10.0000, 

0 . 0000 , 

13 

15  11 

5.6250, 

-10.0000, 

0 . 0000 , 

14 

115  1 

7 . 5000 , 

-10. 0000 , 

0 .  OOCiO , 

15 

1111 

9.3750, 

-10. 0000 , 

C> .  0000 , 

16 

1111 

1 1 . 2500 , 

- 1 0 . 0000 , 

0 » 0000 , 

1  '7 

1111 

13. 1250, 

- 1 0 . 0000 , 

'0 .  'OCiO'  j , 

13 

1111 

1 5 .  OOCiO , 

“ 1 0 . 0000 , 

0 . 0000 , 

^MATEFil  AL 


EX 

i 

,  ,  0,30000E+08 

NUXY  , 

,  ,  O.OOOOOE+00 

*L-DCAi5E 

0 ,  1 , 0 , 0 , 0 , 0  , 

0, 

0. OOOOE-01 , 0. OOOOE 

*LCTITLE 

POINT  LOAD  AT 

BEAM'S  CENTER 

♦SPDISP 

**  SET  ID  = 

1 

1,R0T2 

0.000000 

1  ,UX 

0.000000 

2,UX 

0.000000 

3,UX 

0.000000 

4,UX 

0.000000 

5,UX 

0.000000 

6,UX 

0.000000 

7,UX 

0.000000 

8,UX 

0.000000 

9,UX 

0.000000 

1 0 ,  UY 

0 . 000000 

10,L1X 

0.000000 

11  ,UY 

0 . 000000 

11  ,UX 

0 . 000000 

12, UY 

0 . 000000 

12, UX 

0 . 000000 

13, UY 

0. 000000 

13,UX 

0. 000000 

14, UY 

0.  000000 

14,  UX 

0 . 000000 

15, UY 

0 . 000000 

15,  UX 

0 . 000000 

16,UY 

0.000000 

16,  UX 

0 . 000000 

17, UY 

0. 000000 

17, UX 

0. 000000 

18, UY 

0 . oooooo 

13, UX 

0 . oooooo 

♦CFORCE 

**  SET  ID  = 

1 

1  ,FY 

1 

-0 . 40000E+04 

♦ENDDATA 
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VERIFICATION  PROBLEM  2.20 


COMNAND>GRO.ADD, 1 .0/0/0 

This  adds  a  GRID  (a  point  In  3-D  space),  calls  it 
grid  *1,  and  puts  it  at  coordinates  (0,0,0). 

COMMAND>GRD , TRS .1,2.15 

This  makes  grid  #2  by  "translating"  grid  #1  15  units 
in  the  X  direction  (actually  it  translates  grid  #1  15/0/0 
units  where  the  format  is  X/Y/Z;  since  there  are  no  "0''s 
typed  in.  they  are  automatically  assigned).  If  it  was 
desired  that  grid  #2  be  placed  15  units  in  the  Y 
direction,  the  input  would  be  either  "6RD , TRS , 1 , 2 . 0/15"  OR 
"GRD , TRS , 1 , 2 , /15" .  Default  values  are  "0". 

COMMAND>LIN, 2GD, 1 ,2,1 

This  creates  a  line  starting  at  grid  #1  and  ending  at 
grid  *2,  calling  it  line  *1.  The  format  is 
"LIN, 2GD. starting  grid  *, ending  grid  #,  line 


COMMAND>FEG,BAR, 1 ., ,8/13, 1/1 

FEG  stands  for  "Finite  Element  Generation".  This 
entry  creates  a  group  of  8  bar/beam  elements  (of  element 
type  13),  in  the  global  coordinate  system,  assigning  a 
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■aterlal  ID  *  of  1  and  a  property  ID  *  of  1,  and  a  zoom  of 
1.  Each  eleaent  la  defined  with  a  node  on  each  end  and 
both  nodes  and  eleaents  are  automatically  nuabered  by  the 
program.  This  is  accomplished  using  the  following  format: 
"PE6,BAR.3.4.5.6A/6B/6C,7A/7B/7C"  where: 

Entry  3  is  the  number  of  lines  to  be  "FEGED" 
Entry  4  is  the  list  of  generated  eleaents.  If 
left  blank  the  program  will  automatically 
number  the  eleaents. 

Entry  5  is  the  list  of  generated  nodes.  If 
left  blank  the  program  will  automatically 
number  the  nodes. 

Entry  6A  is  the  number  of  eleaents  to  be 
generated  for  each  line. 

Entry  6B  is  the  NKTP  #,  which  stands  for  the 
element  type  identification  number.  Each 
eleaent  type  has  a  different  NKTP  #.  For  an 
explanation  of  exactly  what  a  specific  NKTP  # 
is,  refer  to  the  element  library  section  of  the 
NISA  II  Users  Manual,  pages  4.0-1  to  4.120-5. 
Entry  6C  is  the  node  ♦  that  has  a  local 
coordinate  system  the  user  may  wish  to  use. 

This  entry  is  used  only  if  the  bean  is  not  to 
be  defined  in  the  global  coordinate  system  (ie 
the  Cartesian,  Cylindrical,  or  Spherical 
systems).  At  tines  it  is  more  convenient  to 
refer  to  items  with  a  local  coordinate  system. 
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for  example  a  force  or  displacement  constraint 
may  be  defined  with  directions  which  are  not 
aligned  with  the  global  coordinate  system.  For 
more  information  on  local  coordinate  systems 
refer  to  pages  3.2-1  to  3.2-4  of  the  NISA  II 
User's  Manual. 

Entry  7A  A  7B  are  the  material  and  property  ID 
#'s  respectfully.  This  set  of  #s  assigned  to 
one  or  more  elements.  Is  used  later  in  the 
Input  phase  for  assigning  actual  properties. 

If  left  blank  the  program  automatically  assigns 
a  material  ID  #  of  1  and  a  property  ID  #  of  1. 
These  numbers  do  not  have  to  be  the  same.  A 
group  of  elements  can  have  the  same  material 
properties  (Poisson's  ratio.  Young's  modulus, 
etc)  with  some  of  the  elements  having  one 
geometric  property  and  others  having  a 
different  geometric  property.  As  a  result  this 
group  of  elements  would  have  one  material  ID  * 
and  two  property  ID  #s. 

Entry  7C  is  the  zoom  factor,  ie  the  length  of 
the  last  element  to  the  length  of  the  first 
element.  If  all  elements  are  the  same  length 
then  the  zoom-1.  If  left  blank  the  program 
automatically  assigns  a  value  of  1.  For  more 
information  and  a  pictorial  representation 
refer  to  page  B-3  of  the  DISPLAY  II  User's 
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Nanua 1 . 


COMMAND>LIN.TRS , 1 , 2 . /-lO 

This  creates  a  second  line,  #2,  by  "translating"  line 
♦1  -10  units  In  the  Y  direction. 

COMMAND>FEG , BAR . 2 , . .8/13 

This  creates  a  group  of  8  bar/beam  eleaents  (of  NKTP 
type  13)  along  line  #2,  each  eleaent  is  defined  with  a 
node  on  each  end,  and  each  node  and  eleaents  Is 
autoaatical ly  numbered  by  the  prograa.  This  was  done 
solely  to  create  the  nodes  required  for  defining  the 
spring  eleaents.  Note,  the  eleaents  of  this  line  are  not 
needed  and  will  be  deleted  later.  Another  method  of 
generating  these  nodes  would  be  by  generating  each 
individual  node  using  the  command  "NOD , ADD , NODE  #,X/Y/Z 
(C00RDINATES)/L0CAL  COORDINATE  SYSTEM  *".  By  using  the 
PEG  command  the  nine  required  nodes,  for  the  required 
springs,  can  be  created  at  the  same  time.  By  using  the 
NOD  command  creating  the  nine  nodes  would  require  multiple 
NOD  commands. 

COMMAND>ELE , DEL . 9T1 6 

This  deletes  elements  9  thru  16.  Even  though  the 
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eleaents  are  deleted  the  nodes  still  remain. 


C0MMAND>ELE,ADD,9,18/1//2,1/10 

This  adds  an  element  #9  with  NKTP-18,  order»l, 
material  ID  #  °'0,  property  ID  #  >2  between  nodes  1  and  10, 
using  the  following  format: 


"ELE , ADD, 3,4A/4B/4C/4D,5A/5B, S"  where: 

Entry  3  is  the  element  ♦  to  be  added.  If  this  is 
left  blank  the  program  will  automatically  assign 
a  number. 

Entry  4A  is  the  NKTP  ♦  (see  entry  6B  under  the 
PEG  description  above). 

Entry  4B  is  the  order  #.  Some  elements  only 
have  one  order,  N0RDR«1,  others  have  many 
depending  on  the  element  type.  For  a  pictorial 
description  of  how  this  order  works/varies  refer 
to  section  4,  the  "elemert  library"  of  the  NISA 
II  User's  Manual;  specifically  NKTP  *s 
1,2,3,4,5,&  32  provide  good  examples  of  some  of 
the  element  types  with  varying  orders. 

Entry  4C  is  the  material  ID  ♦.  Since  these  are 
spring  elements  the  material  characteristics  are 
Included  with  the  property  characteristics  (see 
the  description  of  NKTP  #  18  in  the  NISA  II 


User's  Manual).  Therefore,  this  element  does 
not  need  a  material  ID  number.  Accordingly  It 
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1  8  "  0  "  . 

Entry  4D  is  the  property  ID  #.  (See  entry  7B 
under  the  PEG  description  above.) 

Entry  5A  &  SB  are  the  two  nodes  that  defines  the 
element . 

Entry  6  is  the  vector  number  of  a  previously 
defined  vector  used  in  establishing  a  local 
coordinate  system.  This  can  be  used  for 
elements  of  NKTP  =  38. 

COMMAND>DISP.ADD, 1 , 1 .0/////0 

This  constrains  node  1  against  movement  in  the  X 
direction  and  against  rotation  in  the  Z  direction,  and 
assigns  this  a  constraint  ID  #  of  l  with  the  following 
format:  "DI SP , ADD , 3 , 4 , 5A-5F" ,  where: 

Entry  3  is  the  list  of  nodes  with  the  assigned 
constraint.  If  it  was  desired  that  all  nodes 
from  1  to  10  be  constrained,  with  the  same 
constraint,  the  input  would  be  "ITIO". 

Entry  4  is  the  constraint  ID  #.  The  set  ID 
under  which  the  constraints  are  specified. 

A 

Entry  SA  to  5P  are  the  degree(s)  of  freedom  that 
are  constrained.  This  is  accomplished  by 
placing  a  "0"  in  the  degree  of  freedom  that  must 
be  constrained  with  the  following  format: 
"X/Y/Z/ROTX/ROTY/ROTZ" .  If  the  desired 
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constraints  are  In  the  Y  and  ROTX  directions  the 


input  would  be  either  "/0//0"  or  "/0//0//". 
Consecutive  slashes  indicate  a  null  entry  which 
the  program  identifies  and  either  ignores  or 
places  a  default  value  in  the  location.  In  this 
case  the  program  ignores  the  entry.  For  more 
information  on  null  entries  refer  to  the  general 
command  rules  previously  addressed  in  this 
chapter  or  to  pages  2.1-3  to  2.1-6  of  the 
DISPLAY  II  User's  Manual. 


COMNAND>FORCE . ADD ,1.1. /-4000 

This  entry  adds  a  force  of  -4000  pounds  in  the  Y 
direction,  at  node  1,  with  the  following  format: 

"FOR . ADD , 3 , 4 , 5" ,  where: 

Entry  3  is  the  node  that  the  force  is  applied 
to . 

Entry  4  is  the  force  set  ♦.  A  value  greater 
than  1  is  used  for  multiply  loading  conditions. 
Example:  suppose  the  response  of  a  model  to 
different  loading  conditions  is  desired,  like 

41 

different  load  amounts  at  the  same  node  and/or 
the  same  or  different  forces  at  different  nodes. 
To  do  this  one  force  condition  would  be 
identified  with  one  force  set  #  and  the  other 
force  condition  would  have  a  different  force  set 


#.  This  allows  a  wodel  to  be  analyzed  with 
different  loading  conditions  at  the  same  time 
instead  of  performing  one  analysis  with  one 
loading  condition  and  then  remodeling  the  model 
with  the  different  loading  condition.  In  this 
case  suppose  it  is  desired  that  the  systems 
response  be  determined  for  a  force  of  -4000  lbs 
and  7500  lbs,  both  in  the  Y  direction.  The 
input  would  be:  C0MMAND> "FORCE , ADD . 1 , 1 , /-4000 " 
AND  COMMAND>"FORCE,ADD.1,2./7500" . 

When  the  analysis  is  performed,  answers  will  be 
received  for  both  loading  conditions. 

Entry  5  is  the  force  value  in  the  desired 
direction  X/Y/Z/MX/MY/MZ .  Example:  if  the 
desired  load  condition  is  a  force  of  100  lbs  in 
the  Z  direction  and  a  moment  of  10  ft. lbs.  in 
the  Y  direction  the  input  would  be  either 
"//100//10"  or  "0/0/100/0/10/0".  If  a  position 
is  left  blank  the  program  will  automatically 
place  a  "0"  in  it. 


C0MMAND>PR0P,ADD, 1 ,37o. 25/0/0 

This  adds  a  cross  sectional  area  of  3,  a  moment  of 
Inertia  of  0.25,  a  transverse  shear  coefficient  of  0,  and 
a  local  y-coordlnate  for  stress  recovery  of  0  to  property 
ID  set  #1.  A  local  y-coordlnate  of  "0"  means  that  there 
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is  no  local  coordinate  systea,  rather  the  global 
coordinate  system  is  being  used.  This  conaand  has  the 
following  format  " PROP , ADD . 3A/3B/3C . . . , 4 "  ,  where; 

Entry  3  is  the  property  ID  set  that  the 
properties  will  be  assigned  to. 

Entry  4  is  the  list  of  properties  which  are 
defined  by  the  element  type  with  each  different 
property  being  separated  by  a  /.  (Refer  to  the 
specific  element  type  in  the  element  library 
section  of  the  NISA  II  Users  Manual  for  aore 
info  on  what  properties  are  listed  here). 

Entry  5  is  the  section  orientation  for  the  beam 
element . 


Note:  there  seems  to  be  a  couple  of  problems  with  this 
command.  First  the  input  field  is  a  maximum  of  10 
characters,  and  the  program  automatically  places  up  to 
four  "0"  after  the  decimal  point.  Therefore,  if  the  input 
is  a  five  figure  number  (ex  28125)  the  program  will  list 
••••••••••  in  the  (file  nane).NIS  file.  Second  the 

program  does  not  seem  to  accept  any  numbers  in  the  second 
or  fourth  fields  of  entry  number  3,  automatically 
assigning  its  own  numbers.  Therefore,  these  lines  must  be 
manually  changed  in  the  (file  name).NIS  file  using  a  text 
or  line  editor  like  EDLIN.  The  ( f 1 lename ) . NI S  is  discussed 
in  more  detail  later. 


COMNAND>NAT , ADD , 1 ,30.0E6/0.0 


This  adds  the  Baterlal  properties  of  Youngs  Modulus 
of  elasticity  >  30.0E6  psl  and  Poissons  ratio  =  0.0  to 
material  ID  set  #1  with  the  following  format 
"MAT,ADD.3,4A/4B/4C/. . . "  where: 

Entry  3  is  the  material  ID  set  ♦. 

Entry  4  is  the  list  of  material  properties. 
Like  the  property  characteristics  described 
above,  this  list  will  vary  for  each  element 
type . 


C0NMAN0>DEC 

This  places  the  program  in  the  mode  to  accept 
Executive  Commands.  For  more  info  on  executive  commands 
refer  to  the  previous  Executive  command  overview  and  pg 
7.10-1  of  the  DISPLAY  II  User's  Manual. 

C0MMAND>INS 

This  places  the  program  in  the  mode  to  allow  the 
insertion  of  data  lines. 

COMHAND>ANAL>STATIC 

This  tells  the  program  that  the  analysis  to  be 
performed  is  a  STATIC  analysis. 
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C0NNAND>F1LE>S01 


'  This  provides  the  file  name  prefix  SOI  for  the 

internally  generated  data  files. 

CONMAND>SAVE=26 . 27 

This  tells  the  program  to  save  files  26  and  27, 

^  internally  generated  and  named  files.  The  name  of  these 

files  will  be  S0126.DAT  and  S0127.DAT  respectively.  Note 
if  post  processing  Is  to  be  performed  files  26.dat  and 
^  27.dat  must  be  saved,  which  is  accomplished  using  this 

command.  Refer  to  the  previous  Executive  command  overview 
and  section  3  of  the  NISA  II  User's  Manual. 

I 

COMMAND>*TITLE, ATTEMPT  OF  VERIFICATION  PROBLEM  S20.STD 

This  tells  the  program  that  the  title  of  this 
analysis  is  "ATTEMPT  OP  VERIFICATION  PROBLEM  S20.ST0". 
which  the  program  will  then  print  at  the  top  of  the  answer 
output . 


C0MMAND>*I3 , 1,0, 1,0, 0,0,0, 0 ,0 .OOOOE-01 , 0 .OOOOE-01 

The  *13  is  the  same  as  *LDCASE  which  is  the  load  case 
"control  card"  for  the  desired  output.  Refer  to  page  7.1- 
2  of  the  NISA  II  Users  Manual  for  an  explanation  of  what 
each  number  means.  Note:  there  is  a  problem  with  this 
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coaiand  on  the  PC  version  of  88.7.  No  matter  what  entr 


laced  after  the  oroeraa  autoaatical ly  chanees 


"0, 1 ,3^0,0,4 .0"  on  the  •LDCASE  card  In  the  "(file 


C0MMAND>*I4.1, POINT  LOAD  AT  BEAM'S  CENTER 

The  *14  is  the  same  as  *LCTITLB  which  aeans  load  case 
title.  This  is  the  title  for  the  specific  load  case.  In 
this  example  the  title  "POINT  LOAD  AT  BEAM'S  CENTER"  is 
assigned  to  load  case  1.  If  there  were  multiple  load 
cases,  as  discussed  in  "FORCE , ADD .... "  above,  multiple 
load  case  title  cards  can  be  used,  one  for  each  load  case. 
In  the  description  for  entry  5  of  "FORCE . ADD ... "  two 
loadings  were  discussed,  one  force  and  one  moment.  In 
that  case  the  input  could  be: 

C0MMAND>*r4, 1, POINT  LOAD 
C0MMAND>*I4 . 2 .MOMENT 

or  what  ever  the  desired  titles  were.  The  program  will 
then  place  the  specified  title  at  the  heading  of  the 
answers  for  each  load  case. 
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CONNANO>QUI 


Quits  the  Insert  mode. 

COMNAND>QUI 

Quits  the  DEC  mode. 

COMMAND>WRI . N I QA , SO  1 . N I S 

This  takes  the  above  data  just  input  during  this 
session  and  writes  a  "NISA"  file  called  "SOl.NIS”  that  the 
program  uses  for  it's  analysis. 

COMNAND>NRI .DBS.SOl .BIN 

This  takes  the  above  data  and  writes  a  binary  data 
base  file  called  "SOI. BIN". 

COMMAND>END 

This  ends  the  session. 

A 

Data  Base  File;  At  this  time  the  program  will  ask  if  a 
data  base  file  DISP.BIN  is  to  be  written.  Since  this  has 
already  been  done  with  the  WRI . DBS , SOI . BIN  command  above, 
answer  no  "N". 
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The  Hodel  Is  now  ready  to  be  analyzed.  Refer  to  the 


previous  section  "Program  Usage"  for  guidance  on  how  to 
perform  the  analysis. 
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COMMAND  INPUT  DEMONSTRATION 


In  an  attempt  to  asaiat  the  unfamiliar  uaer  with 
understanding  what  the  various  DISPLAY-PRE  commands 
actually  do.  the  demonstration  below  provides  a  pictorial 
description  of  some  of  the  various  commands  required  In 
the  previous  verification  problem.  This  demonstration 
does  not  contain  all  required  commands  only  a  few  to 
provide  an  Idea  of  how  the  commands  work. 

GRO .ADO . 1 . 0/0/0  This  creates  a  point  at  coordinates 

(0,0,0) . 

(o.o.o) 

□ 

/ 


GRD.TRS , 1 , 2 . 15  This  creates  grid  *2  by  translating 

grid  »1  15  units  in  the  X  direction. 


(o,  o.  o) 
□ 

I 


(/r.o.o) 

□ 

z 
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LIN. 2GD, 1,2,1 


This  creates  line  *1  between  grids  #1 


and  #2. 


D 

I 


/ 


U 

z 


PEG  ,  BAR  ,  1 .  ,  ,  8/13 , 1 /I  This  cr-eates  8  elements  of  NKTP 

type  *13  along  line  #1,  and  assigns  MID-1  and  PlD^l.  Note 
there  is  a  node  created  at  the  end  of  each  element. 


LIN , TRS . 1 , 2 . /-lO  This  creates  line  *2  by  translating 

line  «1  -10  units  in  the  Y  direction. 


PEG . BAR , 2 . , , 8/ 1 3  This  creates  8  elements  of  NKTP  type 

13  along  line  $2  and  allows  the  program  to  automatically 
number  them.  Again,  as  mentlond  above,  each  element  has  a 


node  at  each  end. 


(«•,-»  o) 


ELE.DEL.9T16  This  deletes  elements  9  through  16, 

However,  nodes  10  to  18  remain. 


©©©  @©0®©© 


ELE.A00.9.18/1//2.1/10  This  adds  element  9.  of  NKTP 

type  18  and  PID  #2.  between  nodes  1  and  10. 


000  0©©®0® 


consists  of  basically  the  saae  coamand. 


'4- 


M  t'O  OO  00 


©  © 
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VERIFICATION  PROBLEM  5.5 


TITLE: 

Natural  frequencies  of  a  simply  supported  square  plate 
ELEMENT  TYPE: 

3-D  general  shell  element  (NKTP  ■  20,  NORDR  *  2) 

PROBLEM: 

The  first  five  natural  frequencies  of  a  simply  supported  square  plate,  with  dimensions  as  shown  in 
Fig.  S.5.1.  are  to  be  computed. 

PROPERTIES: 

Material: 

EX  8  10.5E6  psi  (Modulus  of  elasticity) 

NUXY  *  03333  (Poisson’s  ratio) 

DENS  «  2388E-4  lbs.-sec^M>^  (Musdeosity) 

Geometry: 

t  ■  1.0  in  (thkkness) 

FINITE  ELEMENT  MODEL: 

The  piatft  is  modeled  using  riween  j-D  general  shell  elements  (NKTP  *  20,  NC^WR  ■  2).  The  nodes  on 
the  periphery  of  the  plate  are  constrained  in  the  nonnal  qiz)  direction.  In-piaoe  mnslatioos  (UX  and  UY) 
for  aU  tire  nodes  are  caastfalaed  to  Obtain  offiy  the  beafingniodes. 

SOLUTION  PROCEDURE: 

Conventional  subspace  iteration  and  lumped  mass  matrix  fonnulation  are  used  to  extract  the  dgeanralues. 
RESULTS  &  COMPARISON: 

The  first  five  natural  frequencies  of  the  simply  supported  plate  obtaiired  from  NISA  are  compared  to  the 
theoretical  (Ref.  [1])  values  in  Table  53.1.  The  mode  shaqpe  ploo  are  shown  in  Hgs.  53.2  through  5.5.6. 

REFERENCE: 

1.  S.  Timoshenko,  D.  H.  Young,  and  W.  Weaver,  Jr.,  "Vibration  Froblms  in  Engineering".  4th  edition, 
John  WUey  &  Sons,  New  York  (1974). 


GF<D,  ADD,  1 ,0/0/0 

6RD,TR3, 1 ,2,4E3 

LIN,2GD, 1,2,1 

LIN,TRS, 1 ,2,0/48 

PAT,2LN, 1 ,2, 1 

PEG, QUA, 1 , , ,2/4/4/20 

PROP , ADD ,1,1/1/1/1/1/1/1/1 

MAT, ADD, 1 , 10.5E6/0.3333////////2.5a79E-4 

DISP, ADD, 1/9/57/65, 1 ,0/0/0/0/0 

DISP,ADD,2T8/5ST64, 1 ,0/0/0//0 

DISP , ADD , 10/ 14T1 5/23T24/28T29/37T38/42T43/51T52/56 , 
D I SP , ADD , 1 1 T 13/ 1 6T22/25T27/30T36/39T4 1 /44T50/53T55 , 
DEC 
INS 

ANAL==EIGEN 
FILE=£01 
SAVE=26,27 
RESE--=DN  ' 

AUTO^^ON 

£IGE=3UBS,C0NV 
MASS=LUMPED 
♦LDCASE, 1 

*E I GCNTL ,5,0,20,0,0.0,0.0,1. OE-  5 

♦MODEO, 1 ,1,1,1 

*EIGO, 1 , 1 ,0,3, 1 ,0, 1 

QUi 

QUl 

WRI,NISA,E01.NIS 
WRI ,DBS,E01 .BIN 
END 
N 


I 


,0/ 0/0/0 
,0/0 
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*■*  izxrMji  ivfi 
ANrtL-^l;;l'3k:N 
FJLE^^LOl 
SAVE  26,27 
Rt.SE  ON 
AUTO  ON 
EIGE=GLJBS,CONV 
NASS=LUMPED 
*ELTYPE 
1  »  20 , 
*RCTABLE 


1  ,  e 


1 . 0000 

,  1 . 0000 ,  1 . 

'0000 , 

1  -  0000 

,  l.OOOv, 

1 

^ELEMENT 

1, 

1 , 

1,  0 

1  , 

9 

11, 

1/, 

16.  15. 

1  ^  , 

1  , 

1  , 

1  ,  0 

4, 

■ 

12, 

19, 

*  f3  i,  A  /  if 

1  ,1  . 

i , 

1 , 

1  ,  0 

6 , 

7, 

13, 

21  . 

•“>  1  “»  i  C- 

m  1  /  • 

i  . 

1, 

1, 

1,  0 

7, 

8, 

9, 

14, 

1  1 

.i.  • • 

5, 

1, 

1, 

1,  0 

15, 

16, 

17, 

25, 

♦>1  , 

'^0  If  jC'i  f 

24, 

1, 

1  ) 

1,  0 

17, 

18, 

19, 

26, 

33 1 

32,  31, 

25, 

•  7, 

1, 

1, 

1,  0 

19, 

20, 

21, 

27, 

35, 

34 ,  33 , 

26, 

8, 

1  1 

1, 

1,  0 

21 , 

22, 

23, 

28, 

37, 

36 ,  35 , 

27, 

9, 

1, 

1, 

1,  0 

29, 

30, 

31, 

39, 

45, 

44,  43, 

38, 

10, 

1, 

1, 

1,  0 

31, 

32, 

33, 

40, 

47, 

46,  45, 

39, 

11, 

1, 

1, 

1,  0 

wO  f 

34, 

41, 

49, 

48,  47, 

40, 

12, 

1, 

1, 

1,  0 

35, 

36 , 

37, 

42, 

51, 

50,  49, 

41  , 

13, 

1, 

1, 

1,  0 

43, 

44, 

45, 

53, 

59, 

58,  57, 

52 , 

14, 

1, 

1, 

1,  0 

45, 

46, 

47, 

54, 

61, 

60 ,  59 , 

53 , 

15, 

1, 

1, 

1,  0 

47, 

48, 

49, 

,  55, 

63 , 

62,  61, 

54, 

16, 

1, 

1, 

1,  0 

49, 

50, 

51, 

56, 

65, 

64,  63, 

^  vJ  f 

♦NODES 

1 

9  9  9  9 

0.0000 , 

0 . 0000 , 

0. 0000, 

9  9  9  1 

6. 0000, 

0.0000, 

0.0000, 

3 

9  9  9  9 

12.0000, 

0 . 0000 , 

0.0000, 

4 

9  9  9  9 

1 8 . 0000 , 

0 . 0000 , 

0. 0000, 

5 

9  9  9  9 

24 . 0000 , 

0.0000, 

0.0000, 

6 

9  9  9  9 

30 . 0000 , 

0 . 0000 , 

0.0000, 

7 

19  9  9 

36.0000, 

0.0000, 

0.0000, 

8 

9  9  9  9 

42.0000, 

0 . 0000 , 

0.0000, 

9 

9  9  9  9 

48. 0000, 

0.0000, 

0.0000, 
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^  *  *  »  t  1  ♦ 

'■>  .  ( ''ll'j'.l  , 

6.  '..>000 , 

0 .  , 

M  ,,,, 

1 1' .  (;>oo(.) , 

6 . 0000 , 

0 . 0t>(j'0 , 

1*^  ,,,, 

24  .  <.‘>0(jCl  y 

6 . 0'I>00 , 

(■» .  , 

^  T  »  »  5 

36. 0000 y 

6 . 0000 , 

0 .  CjOOO  , 

M  ,,,, 

48 . 0000 , 

6 . 0000 , 

0 . 0000 , 

1^1 

**-*  »  1  1  1 

0 . 0000  y 

1 2 . 0000 , 

0. 0000, 

1  t  »  » 

6. 0000, 

1 2 . 0000 , 

0.0000, 

17  i<,. 

i 2 . 0000 , 

12.0000, 

0 . 0000 , 

IQ  » » » » 

1 3 . 0000 , 

1 2 . 0000 , 

0.0000, 

19  ,,,i 

24 . 0000  y 

12. 0000, 

0 . 0000 , 

20  ,  ,  ,  , 

30 . 0000 , 

1 2 .  '0000 , 

0 . 0000 , 

O  1 

■^  •  1  f  »  »  f 

36 . 0000 , 

1 2 . 0000 , 

0 .  OOO'O , 

— 1  1  f  » 

42.0000 , 

1 2 . 0000 , 

0 .  OO'OO , 

•*•■-•'  Till 

48. 0000 , 

12. 00 '00 , 

0 . 0'OOO , 

24  , ,  ,  . 

0 . 0000 . 

:L8.  OO'OO, 

0 .  OOO'O , 

1 

12. 0000, 

1 a . 0000 , 

0 .  OOOCi . 

*‘^A 

1  •  9 

24 . 0000 , 

1 3 . 0'OOO  , 

'0 . 0000 , 

'  -  rill 

36. 0000, 

1  3.  '0000  y 

0-  OO'OO, 

.iUW  «  ^  4  4 

43.0000 , 

1 3 .  '.J'JOO , 

0 .  O'OOO , 

'■  '.*  \ 

0 « 000C>  - 

24.0000, 

0. 0000 - 

9  19? 

6. 0000 , 

24 . 0'OOO , 

0 . 0000 , 

•^■'1  1  »  !  » 

12.0000, 

24 . 0000 , 

0 . 0000 , 

•-*-=“  1111 

1 a. 0000, 

24.0000, 

0 . 0000 , 

5*11 

24. OOoO, 

24 . 0000  y 

0.0000, 

-'•'^  1  1  »  » 

30 . 0000 , 

24 . 0000 , 

0 . 0000 , 

f  1  f  1 

36. 0000 , 

24.0000, 

0.0000, 

111? 

42.0000, 

24.0000, 

0.0000, 

37  1 t  f 1 

48.0000, 

24.0000, 

0.0000, 

38  ,,,, 

0.0000, 

30.0000, 

0.0000, 

39  ,,,, 

12.0000, 

30.0000, 

0.0000, 

40  t  « 1 1 

24.0000, 

30.0000, 

0.0000, 

41  1  1  y  1 

36.0000, 

30.0000, 

0.0000, 

42  ,,,, 

48.0000, 

30.0000, 

0.0000, 

43  y  y  y  y 

O.OOOOy 

36.0000, 

0.0000, 

44  ,,,, 

6.0000, 

• 0000  f 

0. 0000, 

45  y  y  5  y 

12.0000, 

36.0000, 

0.0000, 

46  y  y  y  I 

18.0000, 

36.0000, 

0.0000, 

47  yyyy 

24.0000, 

36.0000, 

0.0000, 

48  yyyy 

30.0000, 

36.0000, 

0.0000, 

49  yyy, 

36.0000, 

36.0000, 

0.0000, 

50  yyyy 

42. 0000, 

36 . 0':>00  , 

0.0000, 

’-'1  yyyy 

48. 0000, 

36 . 0000 , 

0. 0000 , 

e*o 

*-'-4-  fill 

0.0000, 

42.0000, 

0.0000, 

erir 

1111 

12.0000, 

42.0000, 

0.0000, 

54 

1111 

24 . 0000 , 

42.0000, 

0 . 0000 , 

55  yyyy 

36.0000,  ' 

42.0000, 

0.0000, 

56  yyyy 

48.0000, 

42.0000, 

0. 0000, 

57  yyyy 

0 . 0000 , 

48.0000, 

0.0000, 

58  yyyy 

6. 0000 , 

48.0000, 

0.  ooo':> , 

59  yyyy 

12.0000, 

48.0000, 

0. 0000 , 

60  yyyy 

1 a . 0000 , 

43.0000, 

0 .  O'OOO , 

61  yyyy 

24.0000, 

48.0000, 

0 . 0000 , 

62  yyyy 

30 . 0000 , 

48.0000, 

0.  0000 , 

63  yyyy 

36.0000, 

48.0000, 

0 . 0000 , 

64  yyyy 

42.0000, 

48.0000, 

0.0000, 

65  yyyy 

48.0000, 

40.0000, 

0 .  OOCiCi , 
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♦  MAT[:k  I  t-M 
t-X  ,  1  »  T 

NiJXY  ,  1  ,  . 

DI£NG  ,  1  ,  ,  o.?::, 

>12 

5,0, 20 , 0 , 0 . 0 , 0 . 0 , 1 . 0t:-5 
*SPDISP 

**  SET  ID  =  1 


]  OGOOL.-*  OB 
>33  .,.(;iE>00 
0, 25B7OE-03 


ID  = 

1 

1 ,ROTY, 

0.  OCiOOOO 

1  ,ROTX, 

0.  000000 

1,UZ  , 

0 . 000000 

1,UY  , 

0.000000 

1,UX  , 

0 . 000000 

2,R0TY. 

0. 000000 

2,UZ  , 

0 . 000000 

2,llY  , 

0.  (lO'jOOO 

2,UX  . 

0 .  (j'JOOOO 

3,RDry, 

0 .  OOOO’OO 

3,UZ  , 

0.  OOO'ji'.iO 

3 ,  UY  , 

0 . 000000 

3,u>:  , 

0 . 000000 

4 ,  f^O'VY . 

0 . 0(jCiC>'j0 

4,UZ  , 

0 . 000000 

4,UY  , 

0. 000000 

4,UX  , 

0 . 000000 

5,RQTY, 

0.000000 

5,UZ  , 

0.000000 

5,UY  , 

0.000000 

5,UX  , 

0.000000 

6,R0TY, 

0.000000 

6,UZ  , 

0.000000 

6,UY  , 

0.000000 

6,UX  , 

0.000000 

7,R0TY, 

0.000000 

7,UZ  , 

0.000000 

7,UY  , 

0.000000 

7,UX  , 

0.000000 

Q,ROTY, 

0.000000 

Q,U2  , 

0 . 000000 

e,UY  , 

0.000000 

B,UX  , 
9,R0TY, 
9,RaTX, 
9,UZ  , 
9,UY  , 
‘?,UX  , 

10,R0TX, 
10, U2  , 
10, UY  , 

10,  UX  , 

11,  UY  , 

11,  UX  , 

12,  UY  , 

12,  UX  , 

13,  UY  , 

13,  UX  , 

14, R0TX, 
14, UZ  , 
14, UY  , 
14, UX  , 


0.000000 
0.000000 
0.000000 
0.000000 
0 . 000000 
0.000000 
o.oooopo 
0.000000 
0.000000 
0. 000000 
0.000000 
0 . 000000 
0.000000 
0.000000 
0.000000 
0.000000 
0. 000000 
0.000000 
0.000000 
0. 000000 


lfi,RaTX 

0 . 000000 

Cl.  CiOCiCiOO 

15,  UV 

0 . oooooo 

i5,UX 

Cl.  OCiCiCiOO 

16,  UY 

0 . oooooo 

1 6 ,  UX 

« 

C'.  OOCiOCiO 

r/,uY 

0 . oooooo 

17, UX 

0. 0C>C>CiC>C» 

la.uY 

0 .  CjCiOOOCi 

1 8 ,  UX 

Cl.  CiCiCiCiCiC* 

1 9 ,  UY 

0 .  OCi'OCiOO 

19,UX 

0 .  CiCiOCiCiO 

20 ,  UY 

0 . oooooo 

20,  UX 

0 . oooooo 

21  ,UY 

0 .  OOOClOO 

21  ,UX 

0.  ChjCiOCiO 

22 ,  UY 

0 .  OCi'OOCiO 

22, UX 

0 . oooooo 

23, ROT X 

0 . oooooo 

23,  UZ 

0 . oooooo 

23,  UY 

0 . oooooo 

23,  UX 

0.000000 

24,R0TX 

0.000000 

24,  UZ 

0 . oooooo 

24,  UY 

0.000000 

24,  UX 

0.000000 

2S,UY 

0.000000 

25,  UX 

0.000000 

26,  UY 

0. oooooo 

26,  UX 

0.000000 

27,  UY 

0.000000 

27, UX 

0.000000 

28,R0TX, 

0. oooooo 

28,  UZ 

? 

0.000000 

28,  UY 

f 

0 . oooooo 

28, UX 

1 

0 . oooooo 

29,R0TX, 

0.000000 

29,  UZ 

1 

0.000000 

29,  UY 

1 

0 . oooooo 

29,  UX 

- 

0 . oooooo 

30,  UY 

0 . oooooo 

30,  UX 

0 .  C»00C»0C» 

31  ,UY 

» 

0 . oooooo 

31  ,UX 

V  * 

C>.  OOOCiCiO 

32,  UY 

» 

0 .  OOO'OOO 

32,  UX 

1 

C).  0000C>0 

33,  UY 

1 

0 . OOOOOO 

33,  UX 

0 . OOOOOO 

34,  UY 

0 .  C»00000 

34, UX 

f 

Cl.  C*C>0000 

35,  UY 

Cl .  oooooo 
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3£i,UX  , 

C*.  000000 

, UY  , 

0 . 000000 

36, UX  , 

0 . 000000 

37, ROT X , 

0 . 000000 

37, UZ  , 

<.i .  000000 

37, UY  , 

0 . 000000 

37, UX  , 

0 . 000000 

39 , ROTX , 

0 . 000000 

33. U2  , 

0. 000000 

33, UY  , 

0 . 000000 

33 , UX  , 

0 . 000000 

39, UY  , 

0 . 000000 

39, UX  , 

0 . 000000 

40 , UY  , 

0 . 000000 

40, UX  , 

0 . 000000 

41,UY  , 

0 . 000000 

41, UX  , 

0 . 000000 

42, ROTX, 

0 . 000000 

42, UZ  , 

0.000000 

42, UY  , 

0.000000 

42, UX  , 

0 . 000000 

43, ROTX, 

0.000000 

43, UZ  , 

0.000000 

43, UY  , 

0.000000 

43, UX  , 

0.000000 

44, UY  , 

0 . 000000 

44, UX  , 

0- 000000 

45, UY  , 

0.000000 

45, UX  , 

0,000000 

46, UY  , 

0.000000 

46, UX  , 

0.000000 

47, UY  , 

0 . 000000 

47, UX  , 

0 . 000000 

48, UY  , 

0.000000 

48, UX  , 

0.000000 

49, UY  , 

0.000000 

49, UX  , 

0.000000 

50, UY  , 

0.000000 

50, UX  , 

0 . 000000 

51  ,ROTX, 

0.000000 

51, UZ  , 

0 . 000000 

51, UY  , 

^  0 . 000000 

51, UX  , 

0. 000000 

52, ROTX, 

0. 000000 

52, UZ  , 

0 . 000000 

52, UY  , 

0.000000 

52, UX  , 

0. 000000 

53, UY  , 

0.000000 

53, UX  , 

0.000000 

54, UY  , 

0 . 000000 

54, UX  , 

0.000000 
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55 

,UV 

55 

,UX 

56 

,ROTX 

56 

,UZ 

56 

,UY 

5<^> 

,UX 

57 

,ROTY 

57 

,ROTX 

57 

.UZ 

57 

,UY 

57 

,u>; 

5fi 

,  rsD'i'Y 

5B 

,uz 

5S 

,UY 

53 

,UX 

59,r:GTy , 
59,  LIZ  . 
59 , UY  , 

59,  UX  , 
Z.O,ROTY, 

60,  UZ  , 
60, UY  , 

60,  UX  , 

61 , R0TY, 
61, UZ  , 
61, UY  , 

61,  UX  , 

62, RDTY, 
62, UZ  , 
62 , UY  , 

62,  uX  , 

63,  ROT Y, 
63, UZ  , 
63 , UY  , 

63,  UX  , 

64, R0TY, 
64, UZ  , 
64, UY  , 

64,  UX  , 

65, R0TY, 
65, ROT X, 
65, UZ  , 
65, UY  , 
65,  UX  , 

♦MODEO 

*EI0O 

1 , 0 ,  ,  1  ,  ,  1 
♦ENDDATA 


0 .  <j00000 
o.  ooc>ooo 
0 . 000000 
0 . 000000 
0 . 000000 
0 .  0i;>00L>c> 

Cl .  CiCiOCiCiO 
0 .  OOCiOOCi 
Cl .  OC>OOCiC» 
Cl.  000000 
0.  CiOOCiCm  i 
Cl .  C)  Cl  Cl  Cl  Cl  Cl 
C> .  Cl  Cl  Cl  00  Cl 
Cl .  CiOCiCiOO 
Cl .  CiOCiOOC) 

0 .  CiOCiOCiCi 
0 .  CiCiOCiOO 
0 .  CiCiCiCi'OC* 

Cl .  000000 
0.000000 
0 . 000000 
0 . 000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.  000000 
0.000000 
0.000000 
0.000000 
0.000000 
0 . 000000 
0.000000 
0 . 000000 
Cl .  OCiOCiCiO 
0.000000 
0 . 000000 
C) .  0C>C»C)0C> 
0 . 000000 
0. 000000 
0.  OOOOCIO 
0 . 000000 
0 .  OOCiOOC* 
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VERIFICATION  PROBLEM  5.5 


COMNAND>GRD,ADD, 1 ,0/0/0 

This  adds  a  grid  (a  point  in  3-D  space),  calls  it 
grid  #1,  and  puts  it  at  coordinates  (0,0,0). 

CONMAND>GRD , TRS .1,2,48 

This  nakes  grid  *2  by  "translating"  grid  #1  48  units 
in  the  X  direction.  Actually  it  translates  grid  *1  48/0/0 
units  where  the  format  is  X/Y/Z:  since  there  are  no  "0"s 
typed  in,  they  are  automatically  assigned.  If  it  was 
desired  that  *2  be  placed  48  units  in  the  Y  direction,  the 
input  would  be  either  "GRD ,TRS , 1 , 2 , 0/48"  or 
"GRD . TRS , 1 , 2 , /48" .  Default  values  are  "0". 

COMMAND>L IN , 2GD .1,2,1 

This  creates  a  line  starting  at  grid  #1  and  ending  at 
grid  *2,  calling  it  line  *1.  The  format  is 
"LIN . 2GD . starting  grid  4, ending  grid 

A 

C0MMAND>L IN, TRS. 1,2,0/48 

This  creates  line  *2  by  translating  line  *1  48  units 
in  the  Y  direction.  This  command  could  be  either  be  input 
as  is  or  as  "LIN . TRS , 1 , 2 , /48" .  The  format  is 
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LIN,TRS,3,4,5A/5B/5C"  where: 


Entry  3  In  the  line  to  be  translated. 

Entry  4  is  the  line  *  to  be  created. 

Entry  5  is  the  direction  (X/Y/Z)  entry  3  is  to 
be  translated  in.  Default  values  in  this  entry 
are  "0"8. 


COMMAND>PAT , 2LN , 1 , 2 , 1 

PAT  stands  for  patch,  which  is  defined  as  a  plate  of 
defined  size  in  2  0  space  or  a  surface  of  defined  size  in 
3  D  space..  This  entry  creates  patch  *1  between  lines  #1 
and  #2  with  the  following  forwat:  "PAT , 2LN , LINE  ID  #1, 
LINE  ID  #2.  PATCH  .  Refer  to  the  example  at  the  end  of 
verification  problea  break  down  for  an  example  of  how  a 
is  created  and  what  it  looks  like. 


COMMAND>FEG.QUA, 1 , , .2/4/4/20 

PEG  stands  for  "Finite  Element  Generation".  This 
entry  creates  a  group  of  16  quadrilateral  shell  elements, 

A 

4  elements  by  4  elements,  in  the  plate  created  by  the 
above  patch.  This  is  accomplished  with  the  following 
format;  " PEG , QUA , 3 . 4 , 6 , 6A/6B/6C/6D , 7A/7B/7C/7D  where: 


Entry  3  is  the  patch  ID  number  that  is  to  be 
"FEGED" . 

Entry  4  is  element  ID  list.  This  is  the  list  of 
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the  eleaent  numbers  to  be  created.  If  left 
blank  the  program  will  automatically  number  the 
elements  being  created. 

Entry  5  is  the  list  of  node  numbers  to  be 
created.  Like  the  element  entry  above  if  this 
entry  is  left  blank  the  program  will 
automatically  number  the  nodes. 

Entry  6A  is  the  order  of  the  elements  to  be 
created.  In  this  case  the  elements  are  second 
order  shell  elements,  of  NKTP  type  20.  and 
contain  8  nodes  per  element.  Refer  to  the 
Element  Section  of  the  DISPLAY  II  User's  Manual. 
Entry  6B  and  6C  are  the  number  of  elements  to  be 
generated  along  the  parametric  directions  Cl  and 
C2.  Parametric  directions  are  Internally 
assigned  by  the  program  to  two  edges  of  a  patch 
and  three  edges  of  a  hyperpatch  and  are  used  by 
the  program  for  defining  mesh  density  of  the 
patch.  For  example  suppose  it  is  desired  to 
"break"  a  line  1/3  of  the  way  from  the  end. 

Where  is  the  end  and  where  is  the  start?  The 
program  uses  parametric  directions  to  show  where 
the  start  (where  the  directions  Intersect)  and 
end  is.  Refer  to  pages  2.2-7  to  2.2-10  of  the 
DISPLAY  II  User's  Manual. 

Entry  6D  is  the  NKTP  ♦,  which  stands  for  the 
elements  type  identification  number.  Each 
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element  type  has  a  different  NKTP  ♦.  For  an 
explanation  of  exactly  what  a  specific  NKTP  * 
is,  refer  to  the  element  library  section  of  the 
NISA  II  User's  Manual,  pages  4.0-1  to  4.120-5. 
Entry  7A  and  7B  are  the  material  and  property  ID 
*'s  respectfully.  This  set  of  #8,  assigned  to 
one  or  more  elements,  is  used  later  in  the  input 
phase  for  assigning  act.ual  properties.  If  left 
blank  the  program  automatically  assigns  a 
material  ID  #  of  1  and  a  property  ID  #  of  1.  In 
this  case  all  the  elements  have  the  same 
material  characteristics  and  property 
characteristics  and  therefore  all  have  the  same 
MID  »  and  the  same  PID  #. 

Entry  7C  and  7D  are  the  zoom  values  to  be 
assigned  along  the  parametric  directions 
discussed  in  entry  6B  and  6C  above.  If  left 
blank  the  program  assigns  zoom  -  1  to  both. 

C0MMAND>PR0P,ADD. 1 ,1/1/1/1/1/1/1/1 

A 

This  adds  a  thickness  of  1  unit  to  all  eight 
nodes  of  the  elements  with  PID  =•  1  with  the  following 
format  "PROP , ADD , 3 , 4A/4B/4C . . . ”  where: 

Entry  3  is  the  property  index  #  (PID). 

Entry  4  is  the  list  of  properties  which  are 
defined  by  ttie  element  type  with  each  different 
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property  separated  by  a  /.  Refer  to  the 


specific  element  type  In  the  element  library 
section  of  the  NISA  II  User's  Manual  for  more 
Info  on  what  properties  are  listed  here. 

COMMAND>MAT.ADD, 1 . 10 . 5E6/0 . 3333/ ////// /2 . 5879E-4 

This  adds  the  material  properties  of  modulus  of 
elasticity  (EX)  =  10.5E6  psl,  Poisson's  ratio  (NUXY)  = 
0.333  and  mass  density  »  2.5879E-4  lbs. -sec  squared  per 
Inches  to  the  fourth  power  to  all  the  elements  with 
material  ID  «1.  This  is  accomplished  with  the  following 
format:  MAT , ADD , 3 , 4A/4B/4C . . . "  where: 

Entry  3  Is  the  material  ID  #. 

Entry  4  is  the  list  of  material  properties. 

Like  the  property  characteristics  described 
above,  this  list  will  vary  for  each  element 
type  with  each  separate  property  being  separated 
by  a  slash  "/". 

C0MMAND>DISP, ADD, 1V9/57/65, 1 .O/O/O/O/O 

This  constrains  nodes  1,9,57,  and  65  in  the  X,  Y,  Z, 
Rotation  in  X  and  Rotation  in  Y.  and  assigns  a  constraint 
ID  #  of  1  with  the  following  f  lat: 

"DISP.ADD,3A/3B/3C.  .  .  ,4,5A-5F*'  .  where: 

Entry  3  is  the  list  of  nodes  to  be  constrained. 
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Entry  4  is  the  constraint  ID  #.  The  set  ID 
under  which  the  constraints  are  specified. 

Entry  5A  to  5P  are  the  degree(s)  of  freedom  that 
are  constrained.  This  is  accomplished  by 
placing  a  "0"  in  the  degree  of  freedom  that  must 
be  constrained  with  the  following  format: 
"X/Y/Z/ROTX/ROTY/ROTZ" .  Consecutive  slashes  or 
constrained  directions  not  listed,  indicate  a 
null  entry  which  the  program  either  ignores  or 
places  a  default  value  in  that  place.  In  this 
case  the  program  ignores  the  null  entries.  For 
more  information  on  null  entries  refer  to  the 
general  command  rules  previously  addressed  in 
this  chapter  of  to  pages  2.1-3  to  2.1-6  of  the 
DISPLAY  II  User's  Manual. 

COMMAND>DEC 

This  places  the  program  in  the  mode  to  accept 
Executive  Commands.  For  more  info  on  executive  commands 
refer  to  the  previous  Executl'’e  command  overview  and  page 
7.10-1  of  the  DISPLAY  Ti  User's  Manual. 

C0MMAND>INS 

This  places  the  program  in  the  mode  to  allow  the 


insertion  of  data  lines. 


COMNANO>ANAL>‘EIGEN 

This  tells  the  proeraa  that  the  analysis  to  be 
performed  is  an  EIGENVALUE  analysis. 

C0MMAND>FILE-E01 

This  provides  the  file  name  prefix  EOl  for  the 
internally  generated  binary  data  files. 

C0MMAND>SAVE»26 , 27 

This  tells  the  program  to  save  the  internally 
generated  binary  data  files  26  and  27.  The  names  of  these 
files  mill  be  E0126.dat  and  E0127.dat  respectl vely .  Note 
this  command  must  be  used  to  save  the  data  files  26  and 
27,  which  are  required  for  post-processing  or  dynamic 
analysis . 


COMMAND>RESE>>ON 

This  command  tells  the  program  to  look  at  the 
numbering  of  the  elements  and  resequence  them  (ie  renumber 
then  internally)  to  make  the  finite  element  analysis  more 
efficient.  The  input  of  this  command  is  not  required  as 
the  program  will  automatically  perform  this  function. 

When  creating  a  finite  element  model  there  is  definite 
"pecking  order"  of  how  elements  and  nodes  should  be 
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numbered.  The  analysis  of  a  model  that  does  not  follow 


*  this  order  will  take  longer  than  one  that  does.  This 
command  has  the  program  check  the  model  to  see  If  It 
conforms  to  the  "pecking  order".  If  the  model  does  not 

*  then  the  program  will  Internally  change  the  element 
numbering  so  that  It  does.  The  change  will  not  effect  the 
user  as  It  Is  transparent  to  the  user  and  the  element  IDs 

*  that  the  user  uses  are  not  changed.  For  more  Info  refer 
page  5.3-5  of  the  NISA  II  User's  Manual 


► 


COMMAND>AUTO-ON 


This  command  activates  the  automatic  spurious  normal 
^  rotation  constraint  for  shells.  This  Input  Is  not 

specifically  required  as  the  default  mode  Is  for  this  to 
be  on.  Nodes  In  basic  shell  theory  have  5  degrees  of 
freedom.  When  transferring  a  model  from  the  non-global 
system  to  a  global  system  rotation  In  the  theta  Z 
direction  may  occur.  Additionally  rotation  perpendicular 
to  the  theta  Z  direction  may  occur.  If  this  happens  then 
the  shell  (le  model)  will  not  be  defined  and  therefore 
cannot  be  solved.  « Having  the  AUTO  CONSTRAINT  on  will 
constrain  this  perpendicular  rotation,  thereby  allowing 
the  model  to  be  defined. 
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CONNAND>EIGE-SUBS , CONV 


This  coaaand  tells  the  prograa  that  the  eigenvalue 
extraction  method  is  to  be  conventional  subspace.  This 
command  produces  some  effect  on  the  speed  of  the  analysis. 
If  the  eigenvalues  are  far  apart  then  having  the  method  be 
accelerated  subspace  will  cause  the  program  to 
aggressively  shift  the  frequencies  in  order  to  speed  up 
the  convergence  onto  the  natural  frequencies.  Refer  to 
pages  3.5-1  to  3.5-4  and  5.3-6  in  the  NISA  II  User's 
Manual . 


COMMAND>MASS-LUMPED 

This  command  tells  the  program  to  activate  the  lumped 
matrix  formulation,  ie  it  tells  the  program  to  solve  the 
stiffness  matrix  diagonally  instead  of  consistent  with  the 
way  the  matrix  was  developed.  Refer  to  page  3.5-5  of  the 
NISA  II  User's  Manual. 

COMMAND>*LDCASE . 1 

This  is  the  load  case  "control  card"  for  the  desired 
output.  Refer  to  page  7.1-2  of  the  NISA  II  User's  Manual 
for  an  explanation  of  what  the  nuaber(s)  following  the 
card  heading  mean.  Note  when  performing  an  Eigenvalue 
analysis  this  card  is  not  applicable.  However,  using  this 
card  was  the  only  method  the  author  found  of  "saving"  the 
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dl splaceient  constraints  that  were  input  previousl 
■ore  difficult  method  of  inputing  the  displacement 
constraints  would  be  to  input  them  directly  into  the  NIS 
file  using  a  text  or  line  editor,  like  EDLIIN,  Therefore 
the  user  should  either  use  this  command  to  save  the 
displacement  constraints  and  then  delete  the  two  specific 
load  case  lines  in  the  NIS  file  by  using  a  text  or  line 
editor,  like  EDLIN,  or  input  the  displacement  constraints 
directly  into  the  NIS  file.  (The  NIS  file  will  be 
explained  below.)  The  load  case  lines  that  need  to  be 
deleted  consists  of  one  line  containing  *LDCASE  and  a 
second  line  containing  the  string  of  numbers  input 
immediately  after  the  word  *LOCASE  above. 

COMMAND>*EIGCNTL ,5,0,20,0,0.0,0.0,1 .OE-5 

This  command  tells  the  program  what  the  control 
parameters  are  for  the  eigenvalue  extraction.  Refer  to 
pages  7.1-5  to  7.1-6a  of  the  NISA  II  User's  Manual  for  an 
explanation  of  what  each  number  means. 

A 

C0MMAND>*M0DE0 , 1 , 1 , 1 , 1 

This  command  tells  the  program  what  modes  the 
analysis  is  required  on.  The  format  is  "MODEO , 2 , 3 , 4 , 5 " 
where : 

Entry  2  is  the  first  mode  number  to  be 


analyzed . 

Entry  3  is  the  last  aode  nuaber  to  be  analyzed. 
Entry  4  Is  the  increaent  between  the  first  and 
last  aodes  that  are  to  be  analyzed. 

Entry  5  is  the  output  ID  ♦  that  refers  to  the  ID 
*  listed  in  the  output  coaaand  *EIGO. 

For  aore  info  refer  to  page  7.1-15  of  the  NISA  II  User's 
Manual . 


COMMAND>*EIGO .1,1, 0,3, 1,0,1 


This  coaaand  tells  the  prograa  what  output  options 
to  provide.  For  an  explanation  of  what  each  nuaber  aeans 
refer  to  pages  7.6-2  to  7.5-4  of  the  NISA  II  User's 
Manual . 

Note;  there  seeas  to  be  a  problea  with  the  MODEO  and  EIGO 

coaaands  being  accepted  by  the  DISPLAY-PRE  prograa  and 

then  saved  in  the  NIS.  file.  As  a  result  the  above 

coaaands  aust  be  inputed  directly  into  the  ".NIS"  file 

using  either  a  text  or  line  editor.  This  is  accoaplished 

with  the  following  foraat: 

line  x:  *N0DE0 

line  x-t-1;  aodeo  data 

line  x+2:  •EIGO 

line  x-t-S:  eigo  data 

where  line  x-*-4  is  the  ENDDATA  line 
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The  EIGO  line  must  be  the  first  output  line  listed.  If 
any  other  output  is  requested  (Example  REGIONS,  the 
regions  for  averaged  nodal  stresses  printout;  STRSFIL,  the 
nodal  stress  filtering  data;  etc.;  it  must  be  listed  after 
the  EIGO  data  line.  Refer  to  Section  7.5  OUTPUT  CONTROL 
DATA  of  the  NISA  II  User's  Manual. 

COMMAND>QUI 

Quits  the  insert  node. 


COMMAND>QUI 

Quits  the  DEC  mode. 

COMNAND>WRI . N I SA . EOl . N IS 

This  takes  the  above  data  Just  input  during  this 
session  and  writes  a  "NISA"  file  called  "EOl.NIS"  that  the 
program  uses  for  it's  analysis. 

A 

COMMAND>WRI .DBS, EOl .BIN 

This  takes  the  above  data  and  writes  a  binary  data 


base  file  called  "EOl. BIN 


CONNAND>ENO 


This  ends  the  session. 

Data  Base  File:  At  this  tiae  the  prograa  will  ask  if  a 
data  base  file  DISP.BIN  is  to  be  written.  Since  this  has 
already  been  done  with  the  WRI , DBS , EOl . BIN  coaaand  above, 
answer  no  "N". 

The  node!  is  now  ready  to  be  analyzed.  Refer  to  the 
previous  section  "Proirran  Usage"  for  guidance  on  how  to 
perform  the  analysis. 


2.63 


COMMAND  INPUT  DEMONSTRATION 


In  an  atteapt  to  assist  the  unfaalliar  user  with 
understanding  what  the  various  DISPLAY-PRE  conaands 
actually  do,  specifically  the  creation  of  a  patch,  the 
deaonstration  below  provides  a  pictorial  description  of 
some  of  the  various  commands  required  in  the  previous 
verification  problem.  This  deaonstration  does  not  contain 
all  required  commands  only  a  few  to  provide  an  idea  of  how 
the  coaaands  work. 


GRD, ADD, 1 ,0/0/0  This  creates  a  point  at  coordinates 

(0,0,0) . 


□ 

I 


GRD , TRS , 1 , 2 , 48  This  creates  point  (grid)  42  by 

translating  grid  #1  43  units  in  the  X  direction. 


□ 

I 


□ 

t 
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PAT . 2LN , 1 . 2 . 1  This  creates  patch  tl  between  lines  *1 

and  #2. 


PEG. QUA.l ,, ,2/4/4/20  This  creates  a  group  of  16 

quadrilateral  shell  eleaents  of  NORDR  2,  4  eleaents  by  4 
eleaents,  in  the  patch  created  above.  Due  to  space 
lialtation  not  all  the  nodes  are  nuabered. 
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t  ranuau  uynuaic  Analysis 


NISA  Veriflcition  Manual 


VERIFICATION  PROBLEM 


6.3 


( 


TITLE: 

Transient  response  of  a  beam  subjected  to  distributed  load 
ELEMENT  TYPE: 

2-D  beam  element  (NKTP  •  13.  NORDR  •  1) 

2-D  translational  spring  element  (NKTP  s  18) 

PROBLEM: 

A  beam  supported  by  two  springs  of  stiffhess  K/2  -  n^VTL?  at  its  ends,  as  shown  in  Fig.  6.3.1.  is 
subjected  to  a  distributed  load  given  by: 

p(x)  =  1-2  Ixl/L  for  -IV2  S  X  S  U2 

The  variation  of  the  load  with  time  is  shown  in  Rg.  63.2.  Lateral  motion  of  the  beam  is  constrained.  The 
natural  fiequendes  and  time  history  response  of  the  mid-point  of  the  beam  are  to  be  computed. 

PROPERTIES: 

Material: 

EX  «  3.0x10^  pal  CModnhttofdaitfdty) 

NUXY  »  03  (Poisson’s  nio) 

DENS  »  0.000728  Ibs.-sec^yin'^  (Mass  density) 

Cross-Section: 

A  ■  4.0  in^  (Cron  sectional  iwa) 

I2Z  «  13333  in^  (Moment  ttfitntia  about  beam  2^axis) 

FINITE  ELEMENT  MODEL: 

The  finte  elenent  model  is  shown  in  Hg.  633.  Since  the  tpfBtd  load  and  the  system  are  synunettic,  only 
one  half  of  tire  beam  is  modeled.  Symmetry  boundary  conations  (UX  <■  0.0,  ROTZ  «  OD)  are  applied  at 
node  1.  Lateral  motion  is  constrained  (UX  «  0.0)  at  nodes  2  and  3.  and  node  4  is  fixed  (UX  «  0.0. 
UY-0.0). 


SOLUTION  PROCEDURE: 

Condstent  mass  matrix  formulation  and  subspace  iteration  technique  is  used  in  evaluating  the  natural 
frequencies.  Fim  two  symmetric  modes  are  used  for  the  transient  dynamic  analysis. 

RESULTS  &  COMPARISON: 

The  natural  frequencies  of  the  first  two  symmetric  modes  are  compared  to  the  theoretical  values  in 
Table  6.3.1.  The  peak  transient  response  of  the  mid-poim  of  the  beam  is  compared  to  the  theoretical  value 
in  Table  63.2.  The  displacemem  history  of  the  mid-poim  of  the  beam  is  sho^  in  Hg.  6.3.4. 
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?igur«  6.3.1  Bmb  with  SistributAd  Load 


GRD,ADD, 1 ,0/0/0 

GRD,TRS, 1,2,40 

LIN, 2GD, 1,2,1 

FEG,BAR, 1 ,,,2/13,1/1 

NDD,ADD,4,40/-4 

ELE,ADD,3,18/l//2,3/4 

DISP,ADD, 1 , 1 ,0/////0 

DISP,ADD,2T3,1 ,0 

DISP,ADD,4, 1 ,0/0 

PROP , ADD , 1 , 4 . / 1 . 3333 

MAT ,  ADD  ,  1 , 3 . 0E7  / 0 . 3 /  /  /  /  /  /  /  / 7 . 2aE--4 

PROP , ADD , 2 , 3 . 805E3 

DEC 

INS 

ANAL=EIGENVALUE 
AUTOCONSTRA I NT^OFF 
REBEQUENCING  OF  ELEMENTB=ON 
FILE=ETD01 
SAVE=26,27 

EIGEN  EXTRACTIQN=SUBSPACE 

MASS  MATRIX  F0RMULATI0N=C0N3ISTENT 

*T1TLE, ATTEMPT  OF  VERIFICATION  PROBLEM  ET03.STD 

*12,5,-1,20,0,0.0,0.0, 1 .OE-5 

*13,1 

QUI 

QUl 

WRI,NISA,ETD01.NIS 
WRI , DBS, ETD01.BIN 
END 
N 
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**  executive: 

ANAL=EIGENVALUE 
AUTOCONSTRA I  NT~OF- F 
RESEOIJENCINB  OF  ELEMENTS=Or: 
FILE^ETDOl 


SAVE=26,27 

EIGEN  EXTRACT I ON^SUBSP ACE 

MASS  MATRIX  FORMULATIQN=CONSISTENT 

*TITLE 

ATTEMPT  OF  VERIFICATION  PROBLEM  ET03,STD 
*ELTYPE 

1,  13,  i 

2 ,  13,  i 
■»^R'';TABLE 

1  ,  2 

4.0000,  1.3333, 

2,  1 


3B00. 0000 , 
♦ELEMENT 


1 , 

1 , 

1  5 

1,  0 

1 , 

1 , 

1, 

1,  0 

Xm  • 

•Tf 

0 , 

2  , 

2,  0 

T 

♦NODES 

^  »  5  ^  ^ 

0 .  OOCiCi , 

0 .  C»CiC>C> , 

0 .  OOC>Ci 

n 

^  115  5 

20 .  'OCiOCi , 

0 .  CiOOCi , 

0 .  CKi'OCi 

1111 

4Ci .  OCiCiO , 

0.0000, 

0 . 0000 

^  1111 

40 . 0000 , 

-4.0000, 

0.0000 

♦MATERIAL 

EX  ,  1  ,  ,  0.30000E+08 

NUXY  ,  1,5  0.30000E+00 

DENS  ,  1,5  0.72800E-03 

♦  12 


5,-1 , 20 , 0 ,  Cl .  0 

♦SPDISP 

♦  ♦  SET  ID  =■- 

1 , RQTZ , 
1  ,UX  , 

2, UX  , 

3, UX  , 
4  ,  UY  . 

4, LIX  , 

♦ENDDATA 


0 .0,1. OE-5 
1 

0 .  CiCi'OCiCi'O 
<■) ,  CiCiCi'OCi'j 
C> .  C>C>C>C>00 
Cl .  OOOCi'OC) 
C>«  CiCiCiCiCi’O 
O.  Cl  0  Cl  Cl  Cl  Cl 
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VERIFICATION  PROBLEM  6.3 


C0MMAND>6RD. ADD. 1 , 0/0/0 

This  adds  a  grid  (a  point  in  3-D  space),  calls  it 
grid  #1,  and  puts  it  at  coordinates  (0,0.0). 

COMMAND>GRD , TRS ,1.2,40 

This  aakes  grid  #2  by  "translating"  grid  #1  40  units 
in  the  X  direction.  Actu  .ly  it  translates  grid  ♦!  40/0/0 
units  where  the  format  la  X/Y/Z:  since  there  are  no  "0"a 
typed  in.  they  are  automatically  assigned.  If  it  was 
desired  that  #2  be  placed  40  units  in  the  Y  direction,  the 
input  would  be  either  "GRD , TRS , 1 , 2 , 0/40"  or 
"GRD.TRS,1.2./40".  Default  values  are  "0". 

COMNANDVLIN.EGO, 1 .2. 1 

This  creates  a  line  starting  at  grid  #1  and  ending  at 
grid  *2,  calling  it  line  *1.  The  format  is 
"LIN , 2GD . starting  grid  #. ending  grid 

A 

COMMAND>FEG.BAR, 1 ,,.2/13,1/1 

PEG  ‘:ands  for  "Finite  Element  Generation".  This 
entry  creates  a  group  of  2  bar/beam  elements  (of  element 
*ype  13).  in  the  global  coordinate  system,  assigning  a 


2.71 


Material  ID  #  of  1  and  a  property  ID  #  of  1 ,  and  a  zoom  of 


1.  Each  element  is  defined  with  a  node  on  each  end  and 

I 

both  nodes  and  elements  are  automatically  numbered  by  the 
program.  This  Is  accomplished  using  the  following  format: 
"PE6,BAR,3,4 ,5.6A/6B/6C,7A/7B/7C"  where: 

I 

Entry  3  is  the  number  of  lines  to  be  "FEGED". 
Entry  4  is  the  list  of  generated  elements.  If 
left  blank  the  program  will  automatically  number 
the  elements  . 

Entry  5  is  the  list  of  generated  nodes.  If  left 
^  blank  the  program  will  automatically  number  the 

nodes . 

Entry  6A  is  the  number  of  elements  to  be 
generated  for  each  line. 

I 

Entry  6B  is  the  NKTP  #,  which  stands  for  the 
element  type  identification  number.  Each 
element  type  has  a  different  NKTP  #.  For  an 

I 

explanation  of  exactly  what  a  specific  NKTP  * 
is,  refer  to  the  element  library  section  of  the 
NISA  II  User's  Manual,  pages  4.0-1  to  4.120-5. 
Entry  dC  Is  the  node  #  that  has  a  local 
coordinate  system  the  user  may  wish  to  use. 

A 

This  entry  is  used  only  if  the  beam  is  not  to  be 
defined  in  the  global  coordinate  system  (le  the 
Cartesian.  Cylindrical,  or  Spherical  systems). 

At  tines  It  Is  more  convenient  to  refer  to  items 
with  a  local  coordinate  system,  for  example  a 
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force  or  displaceKent  constraint  Bay  be  defined 
with  directions  which  are  not  aligned  with  the 
global  coordinate  system.  For  sore  information 
on  local  coordinate  systems  refer  to  pages  3.2-1 
to  3.2-4  of  the  NISA  II  User's  Manual. 

Entry  7A  &  7B  are  the  material  and  property  ID 
*'s  respectfully.  This  set  of  #s ,  assigned  to 
one  or  more  elements,  is  used  later  in  the  input 
phase  for  assigning  actual  properties.  If  left 
blank  the  program  automatically  assigns  a 
material  ID  #  of  1  and  a  property  ID  #  of  1. 
These  numbers  do  not  have  to  be  the  same.  A 
group  of  elements  can  have  the  same  material 
properties  (Poisson's  ratio.  Young's  modulus, 
etc.)  with  some  of  the  elements  having  one 
geometric  property  and  others  having  a  different 
geometric  property.  As  a  result  this  group  of 
elements  would  have  one  material  ID  #  and  two 
property  ID  #s. 

Entry  7C  is  the  zoom  factor,  le  the  length  of 
the  last  element  to  length  of  the  first  element. 
If  all  elements  are  the  same  length  then  the 

A 

zoom  «  1.  If  left  blank  the  program 
automatically  assigns  a  value  of  1.  For  more 
information  and  a  pictorial  representation  refer 
to  page  B-3  of  the  DISPLAY  II  User's  Manual. 
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COMMAND>NOO,ADD,4 ,40/-4 

This  coaaand  adds  a  node  at  coordinates  (40, -4,0)  and 
calls  It  node  *  4.  This  Is  accoapllshed  with  the 
following  format:  "NOD . ADD , 3 . 4A/4B/4C/4D"  where: 

Entry  3  Is  the  node  #  to  be  generated. 

Entry  4A.  4B.  &  4C  are  the  (X/Y/Z)  coordinates 
of  the  node. 

Entry  4C  Is  the  local  coordinate  systea  in  which 
the  coordinates  are  to  be  specified.  If  this  Is 
left  blank  then  the  node  Is  in  the  global 
systea . 


COMMAND>ELE , ADD .3,18/l//2,3/4 

This  adds  an  element  #3  with  NKTP  -  18,  order  -  1. 
property  ID  #  >  2  between  nodes  3  and  4,  using  the 
following  format:  "ELE . ADD , 3 , 4A/4B/4C/4D , 5A/5B , 6 "  where: 

Entry  3  is  the  element  ♦  to  be  added.  If  this 
Is  left  blank  the  program  will  automatically 
assign  a  number. 

Entry  4A  is  the  NKTP  ♦  (see  entry  6B  under  the 
PEG  description  above). 

A 

Entry  4B  Is  the  order  #.  Some  elements  only 
have  one  order,  NORDRxl,  others  have  many 
depending  on  the  element  type.  For  a  pictorial 
description  of  how  this  order  works/varies  refer 
to  section  4,  the  "element  library"  of  the  NISA 
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II  User's  Manual:  specifically  NKTP  #s 
1.2,3,4,8.&  32  provide  good  examples  of  some  of 
the  element  types  with  varying  orders. 

Entry  4C  is  the  material  ID  ♦.  This  ♦  is  used 
later  in  the  input  phase  to  assign  actual 
material  properties  to  this  element. 

Entry  4D  is  the  property  ID  #.  (See  entry  7B 
under  the  PEG  description  above.) 

Entry  5A  &  SB  are  the  two  nodes  that  define  the 
element.  Most  elements  require  two  nodes,  one 
at  each  end.  However  since  this  particular 
element  is  a  mass  concentration  it  only  requires 
one  node. 

Entry  6  is  the  vector  number  of  a  previously 
defined  vector  used  in  establishing  a  local 
coordinate  system.  This  can  be  used  for 
elements  of  NKTP  -  11,  12,  &  38. 

COMMAND>DISP , ADD, 1 , 1 ,0/////0 

This  constrains  node  *1  in  the  X  and  ROTZ  directions 
and  assigns  a  constraint  ID  *  of  1  with  the  following 

A 

format  :  ”DISP . ADD , 3A/3B/3C . . . , 4 . 5A-5P "  ,  where: 

Entry  3  is  the  list  of  nodes  to  be  constrained. 
Entry  4  is  the  constraint  ID  ♦.  The  set  ID 
under  which  the  constraints  are  specified. 

Entry  5A  to  5P  ''re  the  degree(s)  of  freedom  that 
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are  constrained.  This  is  acconpllshed  by 
placing  a  "0"  in  the  degree  of  freedom  that  must 
be  constrained  with  the  following  format: 
"X/Y/Z/ROTX/ROTY/ROTZ” .  Consecutive  slashes  or 
constrained  directions  not  listed,  indicate  a 
null  entry  which  the  program  either  ignores  or 
places  a  default  value  in  that  place.  In  this 
case  the  program  ignores  the  null  entries.  For 
more  Information  on  null  entries  refer  to  the 
general  command  rules  previously  addressed  in 
this  chapter  of  to  pages  2.1-3  to  2.1-6  of  the 
DISPLAY  II  User's  Manual. 


COMMAND>PROP.ADD, 1 .4./1 .3333 

This  adds  a  cross  sectional  area  of  4.0  inches 
squared  and  a  moment  of  inertia  on  1.3333  inches  to  the 
forth  to  all  the  elements  with  property  ID  ♦  1.  This  is 
accomplished  with  the  following  format 
"PR0P,ADD,3,4A/4B/4C. . . "  where: 

Entry  3  is  the  property  index  ♦  (PID). 

Entry  4  is  the  list  of  properties  which  are 
defined  by  the  element  type  with  each  different 
property  separated  by  a  /.  Refer  to  the 
specific  element  type  in  the  element  library 
section  of  the  NISA  II  User's  Manual  for  more 
info  on  what  properties  are  listed  here. 
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COMMAND>MAT.ADD, 1 . 3 . 0E7/0 . 3////////7 . 28E-4 


This  adds  the  Material  properties  of  Modulus  of 
elasticity  (EX)  3.0E7  psi,  Poisson's  ratio  (NUXY)  =  0.3, 
and  mass  density  (DENS)  >  7.28E-4  lbs. -sec  squared  per 
inches  to  the  forth.  This  is  accoMplished  with  the 
following  format:"  MAT  .  ADD  ,  3 , 4A/4B/4C/ .  .  . ’’  where: 

Entry  3  is  the  material  ID  st  ♦. 

Entry  4  is  the  list  of  material  properties. 

Like  the  property  characteristics  described 
above,  this  list  will  vary  for  each  element 
type. 


COMMAND>DEC 

This  places  the  program  in  the  mode  to  accept 
Executive  Commands.  For  more  info  on  executive  commands 
refer  to  the  previous  Executive  command  overview  and  page 
7.10-1  of  the  DISPLAY  II  User's  Manual. 

C0MMAND>INS 

This  places  the  program  in  the  mode  to  allow  the 
insertion  of  data  lines. 
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CONNAND> ANAL-EIGENVALUE 


This  tells  the  pT'ograa  that  the  analysis  to  be 
perforaed  is  an  EIGENVALUE  analysis. 

COMMAND>AUTOCONSTRAINT=OFF 

This  conmand  deactivates  the  automatic  spurious 
normal  rotation  constraint.  Since  these  elements  are 
bar/beam  elements  this  command,  which  basically  applies  to 
shell  elements,  is  not  really  applicable.  However,  the 
following  explanation  is  provided  for  what  this  command 
does  if  it  is  activated.  Nodes  in  basic  shell  theory  have 
5  degrees  of  freedom.  When  transferring  a  model  from  the 
non-global  system  to  a  global  system  rotation  in  the  theta 
Z  direction  may  occur.  Additionally  rotation 
perpendicular  to  the  theta  Z  direction  may  occur.  If  this 
happens  then  the  shell  (ie  model)  will  not  be  defined  and 
therefore  cannot  be  solved.  Having  the  AUTO  CONSTRAINT  on 
will  constrain  this  perpendicular  rotation,  thereby 
allowing  the  model  to  be  defined. 

COMNAND>RGSEQUENCING  OF  ELENENTS-ON 

This  command  tells  the  program  to  look  at  the 
numbering  of  the  elements  and  to  resequence  them  (le 
renumber  them  internally)  to  make  the  finite  element 
analysis  more  efficient.  This  command  is  not  really 
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required  as  the  he  prograa  is  set  up  to  autoaatical ly 
perfora  this  function.  When  creating  a  large  finite 
element  aodel  this  process  (ie  RESE-ON)  is  usually 
necessary  because  there  is  a  definite  "pecking  order"  of 
how  elements  and  nodes  should  be  numbered.  The  analysis 
of  a  model  that  does  not  follow  this  order  will  take 
longer  than  one  that  does.  This  command  has  the  program 
check  the  model  to  see  if  it  conforms  to  the  "pecking 
order".  If  the  model  does  not  then  the  program  will 
internally  change  the  element  numbering  so  that  it  does. 
The  change  will  not  effect  the  user  as  it  is  transparent 
to  the  user  and  the  element  IDs  that  the  user  uses  are  not 
changed.  For  more  info  refer  page  5.3-5  of  the  NISA  II 
User's  Manual. 

C0MMAN0>FILB«ETD01 

This  provides  the  file  name  prefix  ETDOl  for  the 
internally  generated  binary  data  files. 

C0MMAND>SAVE-26 , 27 

This  tells  the  program  to  save  the  internally 
generated  binary  data  files  26  and  27.  The  names  of  these 
files  will  be  ETD0126.dat  and  BTD0127.dat  respectively. 
Note  this  command  must  be  used  to  save  the  data  files  26 
and  27.  which  are  required  for  post-processing  or  dynamic 
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analysis  . 


COMMAND>EIGEN  EXTRACT  1 0N> SUBSPACE 

This  coaaand  tells  the  prograa  that  the  eigenvalue 
extraction  Method  is  to  be  accelerated  subspace.  The 
accelerated  selection  of  the  subspace  type  is  the  default 
setting  and  therefore  the  accelerated  did  not  have  to  be 
specified.  This  coanand  produces  some  effect  on  the  speed 
of  the  analysis.  If  the  eigenvalues  are  far  apart  then 
having  the  Method  be  accelerated  subspace  Mill  cause  the 
progran  to  aggressively  shift  the  frequencies  in  order  to 
speed  up  the  convergence  onto  the  natural  frequencies. 
Refer  to  pages  3.5-1  to  3.5-4  and  5.3-6  in  the  NISA  II 
User's  Manual. 

COMMAND>NASS  MATRIX  FORMULATION>°CONSISTENT 

This  coanand  tells  the  prograa  tt  activate  the 
consistent  Matrix  fornulation.  ie  it  tells  the  progran  to 
solve  the  stiffness  Matrix  consistent  with  the  way  the 
Matrix  was  developed.  Using  a  consistent  nass  Matrix 

A 

fornulation  usually  gives  more  accurate  results  than  the 
other  Method  available,  which  is  lunped.  Refer  to  page 
3.5-5  of  the  NISA  II  User’s  Manual. 
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COMMAND>*TITLE .  ATTEMPT  OF  VERIFICATION  PROBLEM  ET03.STD 

This  coaaand  tells  the  program  that  the  title  of  this 
analysis  is  "ATTEMPT  OF  VERIFICATION  PROBLEM  ET03.STD". 
which  will  be  printed  at  the  top  of  the  answer  output. 


COMMAND>*I2 , 5,-l,20,0,0.0.0.0.1.0E-5 

This  command  te’ls  the  program  what  the  control 
parameters  are  for  the  eigenvalue  extraction.  Refer  to 
pages  7.1-5  to  7.1-6a  of  the  NISA  II  User's  Manual  for  an 
explanation  of  what  each  number  means. 

C0MMAND>*I3 . 1 

*13  is  the  same  as  *LDCASE  which  is  the  load  case 
"control  card"  for  the  desired  output.  Refer  to  page  7.1- 
2  of  the  NISA  II  User's  Manual  for  an  explanation  of  what 
the  nunber(s)  following  the  card  heading  mean.  Note  when 
performing  an  Eigenvalue  analysis  this  card  is  not 
applicable .  However,  using  this  card  was  the  only  method 
the  author  found  of  "saving"  the  displacement  constraints 
that  were  input  previously.  A  more  difficult  method  of 
inputing  the  displacement  constraints  would  be  to  input 
them  directly  into  the  NIS  file  using  a  text  or  line 
editor,  like  EDLIIN.  Therefore  the  user  should  either  use 
this  command  to  save  the  displacement  constraints  and  then 
delete  the  two  specific  load  case  lines  in  the  NIS  file  by 
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using  a  text  or  line  editor,  like  BDLIN,  or  input  the 
displacement  constraints  directly  into  the  NIS  file.  (The 
NIS  file  will  be  explained  below.)  The  load  case  lines 
that  need  to  be  deleted  consists  of  one  line  containing 
*LDCASE  and  a  second  line  containing  the  string  of  numbers 
input  immediately  after  the  word  *LDCASE  above. 

COMMAND>QUI 

Quits  the  insert  mode. 

CONMAND>QUI 

Quits  the  DEC  mode. 

COMMAND>HRI ,NISA, ETDOl .NIS 

This  takes  the  above  data  just  input  during  this 
session  and  writes  a  "NISA"  file  called  "ETDOl. NIS"  that 
the  program  uses  for  it's  analysis. 

COMMAND>WRI .DBS, ETDOl .BIN 

This  takes  the  above  data  and  writes  a  binary  data 
base  file  called  "ETDOl. BIN" 


CONNAND>END 


This  ends  the  session. 


Data  Base  File:  At  this  tine  the  prograa  will  ask  if  a 
data  base  file  DISP.BIN  is  to  be  written.  Since  this  has 
already  been  done  with  the  WRI , DBS , ETDOl . BIN  conaand 
above,  answer  no  "N". 

The  Bodel  is  now  ready  to  have  an  eigenvalue  analysis 
performed  on  it.  Once  this  analysis  has  been  cowpleted 
a  dynaaic  analysis  can  be  preformed.  Remember  that  files 
26.dat  and  27.dat  must  be  saved  during  the  eigenvalue 
analysis  in  order  to  be  used  to  complete  the  dynamic 
analysis .  Refer  to  the  previous  section  "Program  Usage" 
for  guidance  on  how  to  perform  the  analysis. 
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anal=transie:nt 

DAMP=VISCOUS 
FILE=ETD01 
STARTING  TIME=0.0 
DELTA  TIME^O.OOOl 
ENDTIME=0.06 
♦TITLE 

ATTEMPT  OF  VERIFICATION  PROBLEM  TD03.STD 
♦MODESECTION 
1,2, 1,1 
♦DAMPING 
1,0. 0 

♦TIMEFUNCTION 

LINEARLY  INCREASING  LOAD 

1 00 1 , 3 , 0 , 1 . 0 

0 .0,0.0,0.02,1. 0 , 0 . 020 1,0. 0 
♦ARRIVAL TIME 
101,0. 0 
♦DPRESSURE 
1,1001,101,1,1 
1 .0,0.5 

2,1001,101,1,1 
0. 5,0.0 
♦RSET 
1001, 1,UY 
♦HISTORY 

0.0,0.08,1,1,1001 

♦ENDDATA 


VERIFICATION  PROBLEM  6.3.  DYNAMIC  INPUT 


ANAL>=TRANSIENT 

This  tells  the  prograa  that  the  analysis  to  be 
performed  is  a  transient  dynamic  analysis. 

DAMP=VISCOUS 

This  tells  the  program  that  the  model  has  viscous 
damping  in  it. 


FILE-ETDOl 

This  tells  program  that  the  prefix  ETDOl  is  to  be 
used  for  the  Internally  generated  binary  data  file.  This 
prefix  must  be  the  same  as  that  used  for  the  eigenvalue 
portion  of  the  problem.  The  reason  for  this  is  because 
the  dynamic  analysis  must  use  the  26.dat  and  27.dat  files 
generated  by  the  eigenvalue  analysis,  for  the  its' 
analysis.  This  command  provides  the  prefix  name  to  look 
for  the  required  data. 


STARTING  TIME-O.O 

This  tells  the  program  that  the  starting  time  of  the 
forcing  function  is  0.0. 
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DELTA  TlNE-0.0001 


This  tells  the  progran  to  analyze  the  response  of  the 
systen  froR  the  starting  tiae  to  the  ending  tiae  in  0.0001 
tiae  increaents.  This  value  aust  be  less  than  the 
difference  between  the  starting  tiae  and  the  ending  tiae 
or  the  prograa  will  not  provide  a  response. 


ENDTIMEx>0.08 

This  tells  the  prograa  that  the  ending  tiae  of  the 
forcing  function  is  0.08. 

Proa  this  point  on  all  coaaands  will  consists  of  aore  than 
one  line  of  data.  As  a  result  a  indicates  the  start 

of  a  new  set  of  data. 

•TITLE 

This  is  a  two  line  data  group.  This  line  tells  the 
prograa  that  the  following  line  of  input  will  be  the 
analysis  title. 


ATTEMPT  OF  VERIFICATION  PROBLEM  TD03.STD 

I 

This  is  the  title  provided  to  the  analysis  which  the 
prograa  will  then  print  at  the  top  of  the  answer  output. 
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•NODESECTION 


This  is  a  multi  line  data  group.  This  line  tells  the 
prograa  the  following  lineCs)  of  input  will  contain  the 
data  specifying  what  wodes  are  to  be  analyzed.  This 
coaaand  is  always  required,  refer  to  page  8.4-2  of  the 
NISA  II  User's  Manual. 

1.2, 1.1 

This  tells  the  prograa  that  nodes  1  and  2,  with  an 
increaent  of  1,  are  to  be  analyzed  with  damping  ID  #1*8 
damping  value.  This  is  accomplished  with  the  following 
format:  Entry  1,2, 3, 4  where: 

Entry  1  la  the  first  mode  to  be  analyzed. 

Entry  2  is  the  last  mode  to  be  analyzed. 

Entry  3  is  the  increment  between  nodes  to  be 
analyzed.  Suppose  the  modes  to  be  analyzed  were 
nodes  1,  5  and  10  then  entry  1  would  be  1,  entry 
2  would  be  10  and  entry  3  would  be  2. 

Entry  4  is  the  damping  ID  #  that  corresponds  to 
the  damping  parameters  set  previously.  The 
damping  ID  #  allows  a  system  to  be  analyzed  for 
various  damping  values  without  remodeling  the 
whole  system  each  tine.  Suppose  it  was  desired 
to  analyze  this  problem  with  two  different 
damping  values,  say  30%  and  70%  of  critical.  To 
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do  this,  two  daaplns  data  groups  would  be 


required.  Part  of  the  daapine  data  input  is  an 
arbitrary  ID  #.  Let's  let  the  30\  daaping  have 
an  ID  #  of  1030'  and  the  70k  have  an  ID  *  of  170 
(note  these  nuabers  are  totally  arbitrary).  Now 
to  find  the  response  of  this  systea  with  those 
two  values  input: 


♦MODESECTION 
1,2,1, 1030 
1,2,1, 170 

♦TIMEPUNCTION 

This  is  a  aulti  line  data  group.  This  line  tells  the 
prograa  that  the  following  lines  of  input  will  provide  the 
data  used  to  specify  the  tiae  function.  This  coaaand  is 
always  required  for  transient  dynaaic  analysis.  Refer  to 
page  8.5-3  of  the  NISA  II  User's  Manual. 


LINEARLY  INCREASING  LOAD 


This  is  the  title  used  to  describe  the  tiae  function. 
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1001,3.0,1.0 


This  tells  the  program  that  three  points  will  be  used 
to  define  the  forcing  function  in  terms  of  abscissa  and 
ordinates,  has  a  scale  of  1.0,  and  a  tine  function  ID  #  of 
1001.  This  is  accomplished  with  the  following  format: 
Entry  1 , 2 . 3 , $ , 4 , 5 , 6  where: 

Entry  1  is  the  time  function  IP  ♦.  This  is  an 
arbitrary  number  that  matches  the  number  used  in 
the  •GROUND,  *DPRESSURE  and/or  •DCFORCE  data 
lines  discussed  later. 

Entry  2  is  the  number  of  points  used  to  define 
the  time  function. 

Entry  3  provides  the  definition  method.  0»  the 
time  function  is  defined  in  terms  of  abscissa 
and  ordinates.  I’-  the  time  function  is  defined 
at  equal  intervals. 

tells  the  program  to  skip  all  entries  between 
the  last  one  inputed  and  entry  4. 

Entry  4  is  the  scale  factor  used  for  the  time 
function.  Basically  there  are  two  numbers  that 
can  be  input  here.  The  correct  number  depends 
on  the  definition  of  the  time  function.  The 

A 

easiest  way  to  explain  this  is  with  an  example. 
Therefore,  let  the  time  function  <■  200  sine  400 
theta.  Now  if  the  amplitude  was  based  on  200 
sine  400  theta  (ie  theta  -  0.6236  radians  and 
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the  force  >  173.2)  then  the  scale  factor  is  1. 

If  the  aaplitude  was  based  on  sine  400  theta  (ie 
theta  >  0.5236  radians  and  the  force  0.866) 
then  the  scale  factor  -  200.  However  any  number 
can  be  used  as  the  scale  factor  as  long  as  the 
user  maintains  consistency  through  out  the 
modal . 

Entry  5  is  the  starting  time  for  the  load 
function . 

Entry  6  is  the  time  increment  used  for  the  time 
function.  If  entry  3  is  "0"  this  entry  is  not 
used . 

0.0,0.0.0.02, 1.0,0.0201.0.0 

This  defines  the  time  function  by  using  a  maximum  of 
three  pairs  of  numbers  per  line,  with  the  total  number  of 
pairs  equal  to  the  number  listed  in  entry  2  above.  The 
pairs  of  numbers  consists  of  the  first  number  referring  to 
the  time  at  point  x  and  the  second  number  of  the  pair  is 
the  function  value  at  point  x. 

*ARRIVALTIME 

This  is  two  line  data  group  that  tells  the  program 
that  the  following  line  of  input  will  specify  arrival 
times  used  in  ground  motion,  concentrated  forces  and/or 


pressure  data. 


101,0.0 

This  tells  the  prograa  that  arrival  time  ID  #  101  has 
an  arrival  time  of  0.0.  This  is  accomplished  with  entry  1 
being  an  arbitrary  number  (referenced  in  the  *GR0UND , 
*DPRESSURE  and/or  *DCF0RCE  data  lines)  assigned  as  the  ID 

#  and  entry  2  being  the  arrival  time. 

•DPRESSURE 

This  is  a  multi  line  data  group.  This  line  tells  the 
program  that  the  following  lines  of  input  will  contain 
data  specifying  the  pressure  loading  of  elements. 

1 , 1001 .101,1,1 

This  tells  the  program  that  there  is  one  additional 
line  of  input  required  to  complete  the  pressure  loading 
data,  that  the  pressure  is  acting  on  face  1  of  element  1. 
and  the  time  function  ID  *  is  1001  and  the  arrival  time  ID 

*  is  ICI.  This  is  accomplished  with  the  following  format: 
Entry  1 , 2 , 3 , 4 , 5 , 6 , 7 , $ , 8  where: 

Entry  1  is  the  element  that  the  pressure  is 
acting  on. 

Entry  2  is  the  tine  function  ID  ♦. 
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Entry  3  Is  the  arrival  tlae  ID  ♦.  Note:  The 
paraaeters  that  entries  2  and  3  refer  to  depend 
on  the  type  of  dynamic  analysis  being  performed. 
For  example  if  a  frequency  response  analysis  is 
to  be  performed  entry  2  refers  to  the  amplitude 
spectrum  ID  *  and  entry  3  refers  to  the  phase 
spectrum  ID  #.  Refer  to  pages  8.6-2  and  8.6-8 
of  the  NISA  II  User's  Manual. 

Entry  4  is  the  face  number  that  the  pressure 
acts  on. 

Entry  5  is  the  number  of  additional  lines 
required  to  complete  the  pressure  loading  data. 
Entry  6  la  the  last  element  number  that  the 
pressure  is  acting  on. 

Entry  7  is  the  increment  between  entry  1  and 
entry  6  elements  the  pressure  is  acting  on. 

$  tells  the  program  to  skip  all  entries  between 
the  last  one  inputed  and  entry  8. 

Entry  8  is  the  uniform  pressure  load  intensity. 
This  entry  is  used  only  if  entry  5-0. 


1 .0,0.5 

This  tells  the  program  that  the  pressure  on  the 
element  referenced  above  (element  #1)  varies  from  a 
pressure  of  1.0  on  the  left  node  to  a  pressure  of  0.5  on 
the  right  node  of  the  element.  Refer  to  page  8.6-9  of  the 
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NISA  II  User's  Manual. 


•RSET 


This  Is  a  Bultl  line  data  group.  This  line  tells  the 
prograa  that  the  following  llne(s)  of  input  will  define 
the  nodal/elewent  response  quantities  that  are  to  be 
provided  by  the  prograa.  This  coaaand  is  required.  Refer 
to  page  8.7-2  of  the  NISA  11  User's  Manual. 


1001 , 1 ,UY 

This  tells  the  prograa  that  the  response  at  node  1, 
in  the  Y  direction,  is  to  be  provided  with  a  response  ID  t 
of  1001.  This  is  accoaplished  with  the  following  foraat: 
Entry  1 , 2 , 3 , 4 , 5 , $ , 6  where: 

Entry  1  is  an  arbitrarily  assigned  response  ID 

*. 

Entry  2  is  the  node  or  eleaent  #  to  be  analyzed. 
Entry  3  is  the  response  to  be  specifically 
provided . 

Entry  4  is  the  last  node  or  eleaent  ♦  to  be 
analyzed.  If  the  response  at  only  one  node  or 
eleaent  is  required  then  this  entry  aust  be  the 
saae  as  entry  2.  If  not  the  prograa  will  state 
that  there  is  an  error  in  this  data  line  and  not 
provide  a  response. 
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Entry  5  la  the  node  or  eleaent  increment.  If 
the  response  at  only  one  node  or  element  is 
required  then  with  entry  4  equaling  entry  2  this 
entry  must  be  one  "1".  If  not  the  program  will 
state  that  there  is  an  error  in  this  data  line 
and  not  provide  a  response. 

$  This  tells  the  program  to  skip  all  entries 
between  the  last  one  Inputed  and  entry  6. 

Entry  6  is  the  layer  number  for  composite 
elements  . 

Note:  Additional  input  lines  can  be  used  to  specify 
responses  for  different  nodes /  elements  and/or  different 
responses.  For  example  suppose  the  responses  at  node  1  in 
both  the  X  and  Y  directions  was  desired.  Then  there  would 
be  two  lines  of  input,  one  with  entry  3  containing  UX  and 
one  with  entry  3  containing  UY. 

•HISTORY 

This  is  a  two  line  data  group.  This  line  tells  the 
program  that  the  following  line  of  input  will  request  the 
response  history  computation  for  printout  and/or  post 
processing.  The  NISA  II  User's  Manual  states  that  this 
command  group  is  optional.  However  the  author  determined, 
by  running  various  analyses,  that  this  command  is  not 
optional .  If  this  command  and/or  the  ‘RSET  command  is  not 
used  then  the  program  ''*11  have  a  "fatal"  error  and  will 


provide  the  following  response  "FATAL  ERROR",  "NO  RESPONSE 
OUTPUT  HAS  BEEN  REQUESTED". 


0.0.0,08.1.1. 1001 


This  tells  the  prograa  that  the  tlwe  history,  uaximxiu 
values  and  tiae  of  occurrence  of  Baxiaun  values  are  to  be 
printed  starting  at  tlae  0.0,  ending  at  time  Is  0.08,  with 
increments  of  1,  and  that  the  response  ID  #  is  1001.  This 
is  accomplished  with  the  following  format: 

Entry  1 , 2 , 3 , 4 , $ , 5-14  where: 

Entry  1  is  the  starting  time.  This  entry  must 
equal  the  value  placed  in  the  "STARTING  TIME" 
Executive  coaaand. 

Entry  2  is  the  ending  time.  This  entry  must 
equal  the  value  placed  in  the  "ENDTINE" 

Executive  command. 

Entry  3  is  the  history  calculation  increment. 
This  increment  must  be  smaller  than  the 
difference  in  time  between  the  start  time  and 
the  end  time.  If  it  is  not  the  program  will 
not  provide  a  response. 

Entry  4  is  the  printout  option.  0  -  response 

4» 

history  calculated  but  not  printed.  1  -  time 
history,  maximum  values  and  time  of  occurrence 
of  maxlauB  values  are  printed.  2  -  only  maximum 
values  and  times  of  occurrence  are  printed. 
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$  tells  the  program  to  skip  all  entries 
between  the  last  one  inputed  and  entry  5. 

Entry  5  to  14  are  the  response  set  ID  ts  that 
the  time  history  response  Is  required.  These 
numbers  should  be  referenced  In  the  *RSET  Input. 


♦ENDDATA 

This  tells  the  program  that  this  is  the  last  line  of 

data . 


Once  this  file  Is  completed  and  the  eigenvalue  analysis 
has  been  performed  the  model  is  ready  to  have  the  dynamic 
analysis  performed.  To  perform  the  dynamic  analysis 
access  the  directory  and  sub-directory  that  the  above  file 
is  saved  in,  type  "NISA2d6"  and  hit  the  return.  Now  move 
the  cursor,  using  the  arrow  keys,  to  the  "DYNAMIC" 
selection  and  hit  the  return.  At  this  time  the  program 
will  ask  for  the  input  file  name.  Type  in  the  name  of  the 
file  above  and  hit  the  return.  Now  the  program  will  ask 
for  the  output  file  name.  Type  in  the  desired  output  name 
and  hit  return.  Now  wait  for  the  analysis. 
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Frequency  Response  Analysis 


NISA  VeriAcadon  Manual 


VERIFICATION  PROBLEM  8.2 


TITLE: 

Frequency  response  analysis  of  a  2-DOF  system 

ELEMENT  TYPE: 

2*D  mass  element  (NKTP  a  25) 

2-D  translational  spring  (NKTP  a  18) 

PROBLEM: 

A  2DOF  system  (Fig.  8.2.1)  is  subjected  to  a  harmonic  force  acting  at  node  2.  Damping  is  neglected. 

The  response  of  the  system  is  to  be  computed  for  an  exciting  frequency  range  of  0.0  to  7.S  HZ  at  intervals 
of  0.25  HZ. 

PROPERTIES: 

Material: 

Ttmslatinnal  mass  ■  200  tb.  flee.^An. 

•IVMflttotMl  ipring  «  0.5  IBAb. 

FINITE  ELEMENT  MODEL: 

The  fiidte  element  model  conrists  of  4  active  nodes  and  S  dements  (2  mass  demaws  and  3  qpdng 
elemews).  ZerodqilacenieotsaiepieaciibedatthacBdpalWB. 

SOLimON  PROCEDURE: 

An  eigenvalue  analysis  using  lumped  mass  fotmulattai  and  convewional  subapace  iteradon  technique  is 
carried  out  first,  fiom  which  2  eigenvalnes  and  aifBBvccton  me  obtained.  The  fteqoency  leaponse  nm  is 
pcrfiOODiCd  sod  V6^^QOS6  (flOEtpUtlldB  Hid  fhlW 

RESULTS  &  COMPARISON: 

The  natural  fiequendes  computed  by  NISA  are  compared  with  the  dteoredcal  resultt  in  Table  8.11.  The 
(Bsdacement  response  results  ate  presented  in  TKble  8.12  dong  with  theoretical  results.  The  amplitude 
and  phase  spectra  for  the  displacement  at  node  2  are  shown  in  Figs.  82.2  and  8.13. 
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figure  8.2.1  Two  Degrees  of  Freedon  System 


e.a  i.see  a.eoe  ^.sse  G.uoa  7  see 

FRtOUEnCY  CH2) 


Figure  8.2.2  Response  saiplitude  Spectrum  for  X'Displscement 
St  node  2 


GRD.ADD, 1 ,0/0/0 
BRDjTRS, 1 ,2,3 
LIN,2BD, 1 ,2,1 
FEB, BAR, 1 ,,,3/18,1/1 
ELE , ADD , 4,25/ 1 /2/2 , 2 
ELE , ADD , 5 , 25/ 1 /2/2 , 3 
DISP,ADD, 1/4, 1 ,0/0 
DISP,ADD,2/3, 1 ,/0 
PROP, ADD, 1,0.5 
PROP, ADD, 2, 200. 

DEC 

INS 

ANAL=EIGEN 

FILE=EF01 

SAVE=26,27 

RESE=OFF 

AUTO=ON 

EIGE=SUBS,CONV 
MASS=LUMPED 
*LDCASE, 1 

*EI GCNTL ,2,0,20,0,0.0,0.0,1. OE-05 

QUI 

QUI 

WRI ,NISA,EF01.NIS 
WRI, DBS, EF01.BIN 


**  EXECUTIVE-; 

ANAL=EIGEN 

EILE=EF01 

SAVE=26,27 

RESE=OFF 

AUT0=0N 

EIGE=SUBS,CONV 

MASS=LUMPED 

*ELTYPE 

1,  18,  1 

2 ,  25 ,  1 

*RCTABLE 

i  ,  1 

200 . 0000 

2,  1 

0 . 5000 
••^ELEMENT 


1 1  1 , 

1  ’  2  , 

1, 

1  ,  0 

2 ,  1  , 

2,  3, 

1, 

1,  0 

-'1  1  5 

3,  4, 

1, 

1,  0 

2, 

2 

2, 

2,  0 

s’  2, 

3 , 

♦NODES 

2, 

2,  0 

1  1  »  1  » 

0 . 0000 , 

0.0000, 

0.0000 

2  1  1  1  » 

1 . 0000 , 

0.0000, 

0. 0000 

•>  Iff) 

2 . 0000 , 

0.0000, 

0.0000 

^  1  1  i  i 

♦  12 

2,0,20,0,0.0,0.0 

♦SPDISP 

♦♦  SET  ID  = 

1,UY  , 

1, ux  , 

2, UY  , 

3, UY  , 

4, UY  , 
4,UX  , 

3.0000, 

,  1 .  OE-05 

1 

0 . 000000 

0 . 000000 

0 . 000000 

0 . 000000 

0 . 000000 

0 . 000000 

0.0000, 

0.0000 

■^ENDDATA 
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VERIFICATION  PROBLEM  8.2 


COMMAND>GRD,ADD. 1 .0/0/0 

This  adds  a  grid  (a  point  in  3-D  space),  calls  it 
grid  *1,  and  puts  it  at  coordinates  (0,0,0). 

CONMAND>GRD . TRS .1,2,3 

This  makes  grid  *2  by  "translating"  grid  *1  3  units 
in  the  X  direction.  Actually  it  translates  grid  *1  3/0/0 
units  where  the  format  is  X/Y/Z:  since  there  are  no  "0"s 
typed  in,  they  are  automatically  assigned.  If  it  was 
desired  that  *2  be  placed  3  units  in  the  Y  direction,  the 
input  would  be  either  "GRO . TRS, 1 , 2 , 0/3"  or 
"GRD , TRS , 1 . 2 , /3" .  Default  values  are  "0". 

C0MNAND>L IN , 2GD .1,2.1 


This  creates  a  line  starting  at  grid  #1  and  ending  at 
grid  #2.  calling  it  line  #1.  The  format  is 
"LIN , 2GD , starting  grid  «, ending  grid  . 

C0MNAND>FEG,BAR,1, , ,3/18,1/1 

PEG  stands  for  "Finite  Element  Generation".  This 
entry  creates  a  group  of  3  bar/beam  elements  (of  element 
type  18),  in  the  global  coordinate  system,  assigning  a 
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■aterlal  ID  *  of  1  and  a  property  ID  #  of  1,  and  a  zoom  of 
1.  Each  element  is  defined  with  a  node  on  each  end  and 
both  nodes  and  elements  are  automatically  numbered  by  the 
program.  This  is  accomplished  using  the  following  format: 
”PBG,BAR,3.4,6,6A/6B/6C,7A/7B/7C”  where: 

Entry  3  is  the  number  of  lines  to  be  "FEGED". 
Entry  4  is  the  list  of  generated  elements.  If 
left  blank  the  program  will  automatically  number 
the  elements. 

Entry  5  is  the  list  of  generated  nodes.  If  left 
blank  the  program  will  automatically  number  the 
nodes . 

Entry  6A  is  the  number  of  elements  to  be 
generated  for  each  line. 

Entry  6B  is  the  NKTP  ♦,  which  stands  for  the 
element  type  identification  number.  Each 
element  type  has  a  different  NKTP  •.  For  an 
explanation  of  exactly  what  a  specific  NKTP  * 
is,  refer  to  the  element  library  section  of  the 
NISA  II  User's  Manual,  pages  4.0-1  to  4.120-5. 
Entry  6C  is  the  node  ♦  that  has  a  local 
coordinate  system  the  user  may  wish  to  use. 

This  entry^  is  used  only  if  the  bean  is  not  to  be 
defined  in  the  global  coordinate  system  (le  the 
Cartesian,  Cylindrical,  or  Spherical  systems). 

At  times  it  is  more  convenient  to  refer  to  items 
with  a  local  coordinate  system,  for  example  a 
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force  or  displacement  constraint  may  be  defined 
with  directions  which  are  not  aligned  with  the 
global  coordinate  system.  For  more  information 
on  local  coordinate  systems  refer  to  pages  3.2-1 
to  3.2-4  of  the  NISA  II  User's  Manual. 

Entry  7A  &  7B  are  the  material  and  property  ID 
t's  respectfully.  This  set  of  «s ,  assigned  to 
one  or  more  elements,  is  used  later  in  the  input 
phase  for  assigning  actual  properties.  If  left 
blank  the  program  automatically  assigns  a 
material  ID  #  of  1  and  a  property  ID  *  of  1 . 
These  numbers  do  not  have  to  be  the  same.  A 
group  of  elements  can  have  the  same  material 
properties  (Poisson's  ratio,  Young's  modulus, 
etc.)  with  some  of  the  elements  having  one 
geometric  property  and  others  having  a  different 
geometric  property.  As  a  result  this  group  of 
elements  would  have  one  material  ID  #  and  two 
property  ID  *s . 

Entry  7C  is  the  zoom  factor,  le  the  length  of 
the  last  element  to  length  of  the  first  element. 
If  all  elements  are  the  same  length  then  the 
zoom  ^  1.  '*  If  left  blank  the  program 
automatically  assigns  a  value  of  1.  For  more 
information  and  a  pictorial  representation  refer 
to  page  B-3  of  the  DISPLAY  II  User's  Manual. 
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COMNANO>ELE.ADD.4.25/l/2/2,2 


i 


This  adds  an  eleaent  #4  with  NKTP  -  25,  order  >=  1, 
■aterial  ID  *  -  2,  property  ID  >  -  2  at  node  2,  using  the 
following  format:  •'ELE,ADD,3.4A/4B/4C/4D,5A/5B,6"  where: 

Entry  3  is  the  element  #  to  be  added.  If  this 
is  left  blank  the  program  will  automatically 
assign  a  number. 

Entry  4A  is  the  NKTP  #  (see  entry  6B  under  the 
PEG  description  above). 

Entry  4B  is  the  order  *.  Some  elements  only 
have  one  order,  N0RDR’>1,  others  have  many 
depending  on  the  element  type.  For  a  pictorial 
description  of  how  this  order  works/varies  refer 
to  section  4,  the  "element  library"  of  the  NISA 
II  User's  Manual:  specifically  NKTP  #8 
l,2,d,4,5,&  32  provide  good  examples  of  some  of 
the  element  types  with  varying  orders. 

Entry  4C  is  the  material  ID  #.  This  *  is  used 
later  in  the  input  phase  to  assign  actual 
material  properties  to  this  element. 

Entry  4D  is  the  property  ID  ♦.  (See  entry  7B 
under  the  PEG  description  above.) 

Entry  5A  A*  5B  are  the  two  nodes  that  define  the 
element.  Most  elements  require  two  nodes,  one 
at  each  end.  However  since  this  particular 
element  is  a  mass  concentration  it  only  requires 
one  node. 
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Entry  6  Is  the  vector  nuaber  of  a  previously 
defined  vector  used  in  establishing  a  local 
coordinate  systea.  This  can  be  used  for 
eleaents  of  NKTP  -  11,  12,  A  38. 

CONMANO>DISP.ADD. 1/4, 1 .0/0 

This  constrains  nodes  1  and  4  in  the  X  and  Y 
direction  and  assigns  a  constraint  ID  *  of  1  with  the 
following  foraat:  "DISP , ADD , 3A/3B/3C . . . . 4 , 5A-5F " ,  where: 

Entry  3  is  the  list  of  nodes  to  be  constrained. 
Entry  4  is  the  constraint  ID  ♦.  The  set  ID 
under  which  the  constraints  are  specified. 

Entry  5A  to  SF  are  the  degree(s)  of  freedoa  that 
are  constrained.  This  is  accoapllshed  by 
placing  a  "0"  in  the  degree  of  freedoa  that  aust 
be  constrained  with  the  following  foraat: 
"X/Y/Z/ROTX/ROTY/ROTZ" .  Consecutive  slashes  or 
constrained  directions  not  listed,  indicate  a 
null  entry  which  the  prograa  either  ignores  or 
places  a  default  value  in  that  place.  In  this 
case  the  prograa  ignores  the  null  entries.  For 
aore  infoTaation  on  null  entries  refer  to  the 
general  coaaand  rules  previously  addressed  in 
this  chapter  of  to  pages  2.1-3  to  2.1-6  of  the 
DISPLAY  II  User's  Manual. 
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CONNAND>PROP , ADO ,1,0.5 

This  adds  a  translational  spring  constant  of  0.5 
lb. /in.  to  all  the  elements  with  property  ID  #  1.  This  is 
accomplished  with  the  following  format 
"PR0P.ADD,3.4A/4B/4C. . . "  where: 

Entry  3  is  the  property  index  #  (PID). 

Entry  4  is  the  list  of  properties  which  are 
defined  by  the  element  type  with  each  different 
property  separated  by  a  /.  Refer  to  the 
specific  element  type  in  the  element  library 
section  of  the  NISA  II  User's  Manual  for  more 
info  on  what  properties  are  listed  here. 


C0NNAND>0BC 

This  places  the  program  in  the  mode  to  accept 
Executive  Commands.  For  more  info  on  executive  commands 
refer  to  the  previous  Executive  command  overview  and  page 
7.10-1  of  the  DISPLAY  II  User's  Manual. 

COMMAND>INS 

This  places  thb  program  in  the  mode  to  allow  the 
insertion  of  data  lines. 


► 
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CONMAND>ANAL-EIGEN 


This  tells  the  prograa  that  the  analysis  to  be 
performed  is  an  EIGENVALUE  analysis. 

C0MMAND>FILE*EF01 

This  provides  the  file  name  prefix  EFOl  for  the 
internally  generated  binary  data  files. 

C0MMAND>SAVE>26 ,27 

This  tells  the  program  to  save  the  Internally 
generated  binary  data  files  26  and  27.  The  names  of  these 
files  mill  be  BF0126.dat  and  EF0127.dat  respectively. 

Note  this  command  must  be  used  to  save  the  data  files  26 
and  27,  which  are  required  for  post-processing  or  dynamic 
analysis . 

COMNAND>RESE-OFF 

This  command  tells  the  program  not  to  look  at  the 
numbering  of  the  elements  nor  to  resequence  them  (ie 
renumber  them  Internally)  to  make  the  finite  element 
analysis  more  efficient.  The  program  is  set  up  to 
automatically  perform  this  function.  Since  this  model  is 
so  small  the  process  is  not  really  necessary.  However 
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when  creating  a  large  finite  element  aodel  this  process 
(ie  RESE*ON)  is  usually  necessary  because  there  is 
definite  "pecking  order”  of  how  elements  and  nodes  should 
be  numbered.  The  analysis  of  a  model  that  does  not  follow 
this  order  will  take  longer  than  one  that  does.  This 
command  has  the  program  check  the  model  to  see  if  it 
conforms  to  the  "pecking  order".  If  the  aodel  does  not 
then  the  program  will  internally  change  the  element 
numbering  so  that  it  does.  The  change  will  not  effect  the 
user  as  it  is  transparent  to  the  user  and  the  element  IDs 
that  the  user  uses  are  not  changed.  For  more  info  refer 
page  5.3-5  of  the  NISA  II  User's  Manual 


COMNAND>AUTO-ON 

This  command  activates  the  automatic  spurious  normal 
rotation  constraint  for  shells.  This  input  is  not 
specifically  required  as  the  default  mode  is  for  this  to 
be  on.  Nodes  in  basic  shell  theory  have  5  degrees  of 
freedom.  When  transferring  a  model  from  the  non-global 
system  to  a  global  system  rotation  in  the  theta  Z 
direction  may  occur.  Additionally  rotation  perpendicular 
to  the  theta  Z  direction  nay  occur.  If  this  happens  then 
the  shell  (ie  model)  will  not  be  defined  and  therefore 
cannot  be  solved.  Having  the  AUTO  CONSTRAINT  on  will 
constrain  this  perpendicular  rotation,  thereby  allowing 
the  model  to  be  defined. 
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COMMAND>EIGE=SUBS.CONV 


This  coanand  tells  the  prograa  that  the  eigenvalue 
extraction  method  is  to  be  conventional  subspace.  This 
coaaand  produces  soae  effect  on  the  speed  of  the  analysis. 
If  the  eigenvalues  are  far  apart  then  having  the  Method  be 
accelerated  subspace  will  cause  the  program  to 
aggressively  shift  the  frequencies  in  order  to  speed  up 
the  convergence  onto  the  natural  frequencies.  Refer  to 
pages  3.5-1  to  3.5-4  and  5.3-6  in  the  NISA  11  User's 
Manual . 


COMMAND>MASS=LUMPED 

This  coMMand  tells  the  program  to  activate  the  lumped 
matrix  formulation,  ie  it  tells  the  program  to  solve  the 
stiffness  matrix  diagonally  Instead  of  consistent  with  the 
way  the  matrix  was  developed.  Refer  to  page  3.5-5  of  the 
NISA  II  User's  Manual. 

COMMAND>*LDCASE , 1 

This  is  the  load  case  "control  card"  for  the  desired 
output.  Refer  to  page  7.1-2  of  the  NISA  II  User's  Manual 
for  an  explanation  of  what  the  nunber(s)  folljwing  the 
card  heading  mean.  Note  when  performing  an  Eigenvalue 
analysis  this  card  is  not  applicable.  However,  using  this 
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card  was  the  only  aethod  the  author  found  of  "saving"  the 
dlsDlaceaent  constraints  that  were  Input  previously.  A 
■ore  difficult  method  of  inputing  the  displacement 
constraints  mould  be  to  input  them  directly  Into  the  SIS 
file  using  a  text  or  line  editor,  like  EDLIIN.  Therefore 
the  user  should  either  use  this  command  to  save  the 
displacement  constraints  and  then  delete  the  two  specific 
load  case  lines  in  the  NIS  file  by  using  a  text  or  line 
editor,  like  EDLIN,  or  input  the  displacement  constraints 
directly  into  the  NIS  file.  (The  NIS  file  will  be 
explained  below.)  The  load  case  lines  that  need  to  be 
deleted  consists  of  one  line  containing  *LDCASE  and  a 
second  line  containing  the  string  of  numbers  input 
immediately  after  the  word  *LDCASE  above. 


CONMAND>*EIGCNTL, 2 ,0.20.0,0.0,0.0, 1 .OE-5 

This  command  tells  the  program  what  the  control 
parameters  are  for  the  eigenvalue  extraction.  Refer  to 
pages  7.1-6  to  7.1-6a  of  the  NISA  II  User's  Manual  for  an 
explanation  of  what  each  number  means. 


C0NNAND>QU1 


Quits  the  insert  mode. 
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CONNAND>QUI 


Quits  the  DEC  Mode. 

COMMAND>WRI , N I SA . EPO 1 . N I S 

This  takes  the  above  data  Just  input  during  this 
session  and  writes  a  "NISA"  file  called  "EFOl.NIS”  that 
the  program  uses  for  it's  analysis. 

COMNAND>WRI .DBS.EFOl .BIN 

This  takes  the  above  data  and  writes  a  binary  data 
base  file  called  "EF01.BIN" 

COMMAND>END 

This  ends  the  session. 


Data  Base  File:  At  this  time  the  program  will  ask  if  a 
data  base  file  DISP.BIN  is  to  be  written.  Since  this  has 
already  been  done  with  the  HRI , DBS , EFOl . BIN  command  above, 
answer  no  "N" . 

The  model  is  now  ready  to  have  an  eigenvalue  analysis 
performed  on  it.  Once  this  analysis  has  been  completed 
a  dynamic  analysis  can  be  preformed.  Remember  that  files 


26.dat  and  27.dat  ■ust  be  saved  during  the  eigenvalue 


analysis  in  order  to  be  used  to  complete  the  dynamic 
analysis .  Refer  to  the  previous  section  "Prograa  Usage" 
for  guidance  on  hoN  to  perform  the  analysis. 
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ANAL.=FREQUENCY 

DAMP= VISCOUS 

FILE=EF01 

SAVE*26,27 

FLOW»0.0 

FUPP=10.0 

BENF=OFF 

INPH=»LINEAR 

INPO»LINEAR 

MRES=OFF 

♦TITLE 

ATTEMPT  OF  VERIFICATION  PROBLEM  8. 
♦MODE 
1,2,1 
♦ADDF 

0 . 0 , 7 . 5 , 30 
♦SPEC 

AMPLITUDE  SPECTRUM 
1 00 1 , 0 , 3 , 1 . 0 , 1 .  O 
PHASE  SPECTRUM 
1 002 , 0 , 3 , 1 . 0 , 0 . 0 
♦DCFO 

2, FX, 1001 ,1002, $,200.0 
♦RSET 

1001 ,2, UX, 3, 1 
♦SPOUT 

1001,0,1,1001 

♦EMDDATA 


VERIFICATION  PROBLEM  8.2.  DYNAMIC  INPUT 


ANAL-FREQUENCY 

This  tells  the  progra*  that  the  analysis  to  be 
performed  is  a  Frequency  Response  analysis. 


DAMP-VISCOUS 

This  tells  the  program  that  the  model  has  viscous 
damping  in  it. 


I 


i 


FILE-EFOl 

This  tells  program  that  the  prefix  EFOl  is  to  be  used 
for  the  internally  generated  binary  data  files.  This 
prefix  must  be  the  same  as  that  used  for  the  eigenvalue 
portion  of  the  problem.  The  reason  for  this  is  because 
the  dynamic  analysis  must  use  the  26.dat  and  27.dat  files 
generated  by  the  eigenvalue  analysis,  for  the  its' 
analysis.  This  command  provides  the  prefix  name  to  look 
for  the  required  data. 

SAVE-26,27 

This  command  tells  the  program  to  save  the  internally 
generated  data  files,  26.dat  and  27.dat.  Note  this 
command  la  not  required,  as  the  dynamic  program 
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autoaatical ly  saves  these  two  files. 


PLOW-0.0 

This  tells  the  prograa  that  the  lower  liait  of  the 
exciting  frequencies  is  0.0  cycles  per  tine. 


FUPP-10.0 

This  tells  the  program  that  the  upper  limit  of  the 
exciting  frequencies  is  10.0  cycles  per  time. 

Note:  In  a  Frequency  Response  analysis  the  program 
determines  the  response  of  the  system  at  various  exciting 
frequencies  ranging  from  the  lower  limit  (FLOW  value)  to 
the  upper  limit  (FUPP  value). 


GENF-OFF 

This  command  tells  the  program  to  suppress  the 
generation  of  frequency  points  between  nodes.  A  system 
with  more  than  one  degree  of  freedom  has  more  than  one 
natural  frequency.  “The  first  node  (or  fundamental  node) 
refers  to  the  node  associated  with  the  lowest  natural 
frequency.  During  a  frequency  response  analysis  the 
program  determines  the  response  of  the  system  to 
frequencies  varying  between  those  set  with  the  upper  and 
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lower  limit  commands.  With  this  command  set  to  off 
(GENF-OFF)  only  the  response  to  the  natural  frequencies 
occurring  between  the  upper  and  lower  limits  will  be 
analyzed.  However  if  this  command  is  on  (GENF»ON)  then 
the  program  will  also  determine  the  system's  response  to 
various  other  frequencies  equally  spaced  between  the 
natural  frequencies  within  the  upper  and  lower  limits. 

INPH-LINEAR 

This  command  tells  the  program  that  the  effect  of  any 
specified  phase  angle  is  computed  along  a  linear  ordinate 
and  abscissa.  The  other  choices  would  be  to  have  the 
ordinate  and/or  the  abscissa  be  logarithmic. 

INPO«LINEAR 

This  command  tells  the  program  that  the  effect  of  any 
specified  power  spectral  density,  amplitude  spectrum  or 
shock  spectrum  is  computed  along  a  linear  ordinate  and 
abscissa.  The  other  choices  would  be  to  have  the 
ordinate  and/or  the  abscissa  be  logarithmic. 


NRES-OFF 

This  tells  the  program  not  to  (le  to  suppress  modal 
response)  printout  the  displacement,  velocity  and 
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acceleration  responses  for  all  aodes.  As  a  result  only 
those  nodes  specified  will  have  responses  printed. 


From  this  point  on  all  coanands  will  consists  of  nore  than 
one  line  of  data.  As  a  result  a  indicates  the  start 

of  a  new  set  of  data. 

•TITLE 

This  is  a  two  line  data  group.  This  line  tells  the 
program  that  the  following  line  of  input  will  be  the 
analysis  title. 

ATTEMPT  OP  VERIFICATION  PROBLEM  8.2 

This  is  the  title  provided  to  the  analysis  which  the 
progran  will  then  print  at  the  top  of  the  answer  output. 


•MODE 

This  is  a  nulti  line  data  group.  This  line  tells  the 
progran  that  the  following  line(s)  of  input  will  contain 
the  data  specif ying 'what  nodes  are  to  be  analyzed.  This 
connand  is  always  required,  refer  to  page  8.4-2  of  the 
NISA  II  User's  Manual. 


I 
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1.2.1 


This  tells  the  proeraa  that  aodes  1  and  2.  with  an 
increaent  of  1  between  the*,  are  to  be  analyzed  with  zero 
"0”  dawping.  This  is  accowplished  with  the  following 
format:  Entry  1,2, 3, 4  where: 

Entry  1  is  the  first  mode  to  be  analyzed. 

Entry  2  is  the  last  mode  to  be  analyzed. 

Entry  3  is  the  increment  between  modes  to  be 
analyzed.  Suppose  the  nodes  to  be  analyzed  were 
nodes  1 .  5  and  10  then  entry  1  would  be  1 ,  entry 
2  would  be  10  and  entry  3  would  be  2. 

Entry  4  is  the  damping  identification  *  that 
corresponds  to  the  damping  parameters  set 
previously.  The  damping  ID  *  allows  a  system  to 
be  analyzed  for  various  damping  values  without 
remodeling  the  whole  system  each  tine.  In  this 
example  there  is  zero  ”0”  damping.  Therefore, 
there  is  no  damping  data  input.  However, 
suppose  it  was  desired  to  analyze  this  problem 
with  two  different  damping  values,  say  30%  and 
70%  of  critical.  To  do  this  two  damping  data 
groups  would  be  required.  Part  of  the  damping 
data  input  is  an  arbitrary  ID  #.  Let's  let  the 
30%  damping  have  an  ID  *  of  1030  and  the  70% 
have  an  ID  •  of  170  (note  these  numbers  are 
totally  arbitrary).  Now  to  find  the  response  of 
this  system  with  those  two  values  input: 


2.118 


•MODE 


1,2,1, 1030 
1.2,1,170 


•ADDF 


This  is  a  two  line  data  group.  This  line  tells  the 
program  that  the  following  line  of  input  will  contain  the 
data  specifying  the  exciting  frequencies  that  are  to  be 
added  to  the  system. 


0.0,7.5,30 


This  tells  the  program  that  the  system  is  to  be 
analyzed  for  30  exciting  frequencies  equally  spaced 
between  0.0  cycles  per  time  and  7.5  cycles  per  time.  This 
is  accomplished  with  the  following  format:  Entry  1,2,3 
where : 

Entry  1  is  the  value  of  the  starting  frequency. 
Entry  2  is  the  value  of  the  ending  frequency. 
Entry  3  is  the  number  of  equal  increments 
between  entry  1  and  2  that  are  to  be  analyzed. 

•SPEC 


This  is  a  multi  line  data  group.  This  line  tells  the 
program  that  the  following  lines  of  input  will  provide  the 
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spectruB  definition.  This  is  accoaplished  with  aultiple 
two  line  series  where  the  first  line  provides  the  spectrua 
description  and  the  second  line  provides  the  data. 

AMPLITUDE  SPECTRUM 

This  tells  the  prograa  that  the  following  line  of 
data  defines  the  aaplitude  spectrua. 

1001,0,3,1.0.1.0 

This  tells  the  prograa  that  the  aaplitude  spectrua 
has  a  constant  value  of  1.0,  a  scale  factor  of  1.0  and  an 
ID  «  of  1001.  This  is  accoaplished  with  the  following 
foraat:  Entry  1,2, 3, 4, 5  where: 

Entry  1  is  an  arbitrarily  assigned  #. 

Entry  2  is  the  nuaber  of  points  used  to  define 
the  spectrua.  If  the  spectrua  is  constant  this 
entry  is  ignored.  If  the  spectrua  is  not 
constant  this  entry  aust  be  greater  than  or 
equal  to  two. 

Entry  3  is  a  nuaber  between  0  and  3  which  is 
used  to  define  the  input  spectrua.  "0"  - 
spectrua  is  a  function  of  frequency 
( cycles/tlae ) ;  "1”  -  spectrua  is  a  function  of 
frequency  ( radlans/tiae ) ;  "2"  -  spectrua  is  a 
function  of  period:  ”3"  -  spectrua  is  constant. 


2 . 120 


data . 


Entry  5  la  the  constant  apectrua  value. 

Note  if  the  apectrua  la  not  constant  then  the  next  llne(8) 
of  Input  would  define  the  apectrua.  This  would  be 
accoaplished  with  each  line  containing  a  aaxlaua  of  three 
pairs  of  nuabers.  The  first  nuaber  of  the  pair  is  the 
frequency/period  value.  The  second  nuaber  is  the  apectrua 
value  that  corresponds  with  the  frequency/period  value. 

The  total  nuaber  of  palra  auat  equal  the  value  placed  in 
entry  2  above . 

PHASE  SPECTRUM 

This  tells  the  prograa  that  the  following  line  of 
data  defines  the  phase  spectrua. 

1002 ,0,3, 1.0. 0.0 

This  tells  the  prograa  that  the  phase  spectrua  has  a 
constant  value  of  0.0,  a  scale  factor  of  1.0  and  has  an  ID 
#  of  1002.  Refer  to  the  above  description  for  foraat  and 
entry  discussions.  ' 

For  every  line  of  Aaplltude  Spectrua  data  there  has  to  be 
a  corresponding  line  of  Phase  Spectrua  data.  The  first 
entry  of  the  data  lines  aust  be  different.  For  exaaple 
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suppose  there  are  two  dynawic  forces  being  applied  to  a 
■odel.  One  of  the  forces  is  out  of  phase  with  the  other 
force  by  90  degrees.  This  would  be  wodeled  as  follows: 
Amplitude  Spectrum 
1001 .0,3,1,! 

Amplitude  Spectrum 
1002 . 0 , 3 , 1  , 1 
Phase  Spectrum 
1003,0,3, 1 .0 
Phase  Spectrum 
1004,0, 3,1,1.571 

The  description  "Amplitude  Spectrum"  or  "Phase  Spectrum" 
must  proceed  every  line  of  data.  Also  the  Phase  Spectrum 
data  line  1003  corresponds  to  Amplitude  Spectrum  data  line 
1001  and  1004  corresponds  with  1002.  Refer  to  the  above 
description  for  an  explanation  of  what  each  line  of  data 
means . 

•DCPO 

This  is  a  multi  line  data  group.  This  line  tells  the 
program  that  the  following  line(8)  of  input  will  contain 
the  data  used  to  specify  a  concentrated  nodal  force. 

2. FX, 1001, 1002, $.200.0 

This  tells  the  program  there  is  a  concentrated  force 
of  200.0  units  in  the  X  direction  at  node  2  for  both  the 
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amplitude  and  phase  spectruas.  This  is  acconplished  with 
the  following  format:  Entry  1 , 2 . 3 , 4 , 5 , 6 , $  ,  7  where: 


I 


I 


I 


» 


Entry  1  is  the  node  number  the  force  or  moment 
is  applied  to. 

Entry  2  is  the  force  or  moment  direction. 

Entries  3  A  4  are  the  ID  #8  that  corresponds  to 
the  amplitude  and  phase  spectrums. 

Entry  5  is  the  last  node  that  the  force  or 
moment  is  applied  to. 

Entry  6  is  the  increment  between  the  first  and 
last  nodes  this  force  or  moment  will  be  applied 
to . 

$  This  tells  the  program  to  skip  all  entries 
between  the  last  one  inputcd  and  entry  7. 

Entry  7  is  the  value  of  the  force  multiplier. 
Basically  there  are  two  numbers  that  can  be 
input  here.  The  correct  number  depends  on  the 
spectrum  definition.  The  easiest  way  to  explain 
this  is  with  an  example.  Therefore,  let  the 
forcing  function  »  200  sine  400  theta.  Now  if 
the  amplitude  definition  was  based  on  200  sine 
400  theta  (ie  theta  =  0.5236  radians  and  the 
force  =  173.2)  then  the  force  multiplier  is  1. 
However  if  the  amplitude  was  based  on  sine  400 
theta  (le  theta  =  0.5236  radians  and  the  force  = 
0.866)  then  the  force  multiplier  »  200. 

However,  any  number  can  be  used  as  the  force 


■ultiplier,  as  long  as  the  user  Maintains 


consistency  through  out  the  modal. 


•RSET 


This  is  a  multi  line  data  group.  This  line  tells  the 
program  that  the  following  line(s)  of  input  will  define 
the  nodal/element  response  quantities  that  are  to  be 
provided  by  the  program.  This  command  is  required  and 
there  must  be  a  minimum  of  one  line  of  *RSET  data  for  each 
line  of  AMPLITUDE  SPECTRUM.  Refer  to  page  8.7-2  of  the 
NISA  II  User's  Manual. 


1001, 2, UX, 3.1 

This  tells  the  program  that  the  response  at  nodes  2 
and  3,  in  the  X  direction,  is  to  be  provided  with  a 
response  ID  «  of  1001.  This  is  accomplished  with  the 
following  format:  Entry  1 , 2 , 3 , 4 . 5 , $ . 6  where: 

Entry  1  is  an  arbitrarily  assigned  response  ID 
* . 

Entry  2  is  the  node  or  element  #  to  be  analyzed. 
Entry  3  is  the  response  to  be  specifically 
provided . 

Entry  4  is  the  last  node  or  element  ♦  to  be 
analyzed . 

Entry  5  is  the  node  or  element  increment 


(between  the  entry  2  node/eleaent  and  entry  4 
node/eleaent ) . 


4  This  tells  the  program  to  skip  all  entries 
between  the  last  one  inputed  and  entry  6. 

Entry  6  is  the  layer  nuaber  for  composite 
elements . 

Note:  Additional  input  lines  can  be  used  to  specify 
responses  for  different  nodes/eleaents  and/or  different 
responses.  For  example  suppose  the  responses  at  node  2  in 
both  the  X  and  Y  directions  was  desired.  Then  there  would 
be  two  lines  of  input,  one  with  entry  3  containing  UX  and 
one  with  entry  3  containing  UY. 


•SPOUT 

This  is  a  multi  line  data  group.  This  line  tells  the 
program  that  following  line(8)  of  data  will  provide  the 
required  output  responses. 


1001,0, 1 , 1001 

This  tells  the  program  that  the  output  spectr, 
response  maxima  and<  the  frequencies  at  which  they  occur 
are  to  be  printed,  with  a  zero  "0"  phase  angle  shift,  for 
response  10  set  1001.  This  output  will  have  an 
identification  nuaber  1001,  which  is  referenced  during 
post  processing.  This  is  accomplished  with  the  following 


I  format: 

t 


Entry  1,2, 3, $.4-13  where: 

Entry  1  is  the  output  identification  ♦  to  be 
referenced  during  post  processing. 

Entry  2  is  the  phase  angle  shift. 

Entry  3  is  the  printout  option:  "0"  -  output 
spectra  is  calculated  but  only  the  response 
■axiaa  and  the  frequencies  at  which  they  occur 
are  printed:  "1"  »  output  spectra,  response 
■axiaa  and  the  frequencies  at  which  they  occur 
are  printed.  "2"  =  the  response  is  printed  at 
the  specified  phase  angle  shift. 

$  tells  the  prograa  to  skip  all  entries  between 
the  last  one  inputed  and  entry  4. 

Entry  4  to  13  are  the  response  set 
identification  noBber(s)  (up  to  ten)  for  which 
the  output  is  required.  This  (these)  nuBber(s) 
correspond  to  the  response  output  listed  in  the 
RSET  coaaand. 


*ENDDATA 

This  tells  the  program  that  this  is  the  last  line  of 
data.  < 

Once  this  file  is  completed  and  the  eigenvalue  analysis 
has  been  performed  the  model  is  ready  to  have  the  dynamic 
•  analysis  performed.  To  perform  the  dynamic  analysis 
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access  the  directory  and  sub-directory  that  the  above  file 
is  saved  in.  type  "NISA286”  and  hit  the  return.  Now  move 
the  cursor,  using  the  arrow  keys,  to  the  "DYNAMIC" 
selection  and  hit  the  return.  At  this  tiae  the  prograa 
will  ask  for  the  input  file  naae.  Type  in  the  naae  of  the 
file  above  and  hit  the  return.  Now  the  prograa  will  ask 
for  the  output  file  naae.  Type  in  the  desired  output  naae 
and  hit  return.  Now  wait  for  the  analysis. 
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CHAPTER  THREE 

PROGRAM  UTILIZATION  AND  COMPARISON 


INTRODUCTION 


Chapter  Two  provided  a  user's  aanual  for  the  NISA  II 
program  as  well  as  a  detailed  step  by  step  explanation  of 
how  to  generate  a  finite  element  model  and  analysis  it  for 
Static,  Eigenvalue,  Transient  Dynamic,  and  Frequency 
Response  conditions.  This  "user's  manual"  will  now  be 
utilized  to  model  and  analyze  two  different  problems,  a 
rigid  frame  on  a  slab  on  springs  and  a  circular  foundation 
on  a  halfspace. 


I 
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FRAME  WITH  SLAB  ON  SPRINGS 

The  following  example  consists  of  a  Rigid  frame  on  a 
concrete  slab  on  springs.  This  example  was  generated 
using  DISPLAY  II  and  then  analyzed  for  both  a  static 
response  and  a  frequency  response  using  the  NISA  II 
program.  This  analysis  is  presented  in  the  following 
order:  Model  sketch,  DISPLAY  II  input,  input  explanation, 

frequency  response  input,  input  explanation,  and  response 
output.  A  GTSTRUDL  response  is  currently  being  analyzed 
and  will  be  placed  into  the  Appendix  D  of  this  report,  for 
a  comparison  with  NISA  II.  when  it  is  completed. 
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This  exaaple  consists  of  a  rigid  fraae  on  a  concrete  slab, 
on  springs.  Machinery  located  in  the  center  of  the  rigid 
frame  causes  a  sinusoidal  dynamic  force  of  60,000  #  in  the 
X  direction  and  30,000  *  in  the  Z  direction.  These  two 
forces  are  evenly  distributed  between  two  points,  such 
that  Fx«30,000  #  and  Fz^lS.OOO  *  at  each  point,  are 
occurring  at  a  frequency  of  400  rpm,  and  Fz  is  90  degrees 
out  of  phase  with  Fx .  The  weight  of  the  machinery  is 
distributed  at  the  top  of  the  rigid  frame  as  follows; 


lOK  2flK  lOK 


lOK  20K  lOK 


Material  properties  are  as  follows: 
columns  (cl)=  W18xll9  beams  (bl)>=:  W36x230 

concrete  slab  3ft  thick  (b2)-  W30xl73 

poissons  ratio  >  0.17 
E  -  3.6E6  psl 
gamma  >  150  pcf 

Springs  (kips  per  inch) 

K1  -  586  K2  -  1170  K3  -  2340  K4  -  2500  K5  -  5000 
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GRD , ADU , i , 0 / 0 / 0 
GRD, TUB. 1 ,2,240 
LIN,2GD, 1 ,2 ' 1 
LIN,TRS, 1 ,2,/iaO 
PAT,2LN, 1,2,1 
VEW , ABB , -BO/0/ 1 0 
PLO 

PAT, IRS. 1 ,2,// -20 
l-EG  ,  QUA  /2/2/20 , 1  /  1 

PEG. QUA, 2. , , 1 /2/2/20  1/1 
ELE,DEI.,5T3 
NDI) ,  ADD  ,  19,  /  -  20 
NUD , ADD , 20 , 1 20 / “20 
NOD, ADD, 21 ,240/ “20 
ELE,ADD,5, 17/1//2, 1/10 
ELE, ADD,6, 17/l//3,2/l 1 
ELE,ADD,7, 17/ 1 //2 , 3/ 12 
ELE,ADD,8, 17/l//3,4/13 
ELE,ADD,9, 17/l//4,5/14 
ELE,ADD,10, 17/1 //3, 6/ 15 
ELE,ADD, 11 , 17/l//2,7/16 
ELE,ADD, 12, 17/17/3,8/ 17 
ELE,ADD,13, 17/1 //2, 9/18 
ELE, ADD, 14, 17/1/ /5, 1/19 
t ADD ,15,17/1/ /&, 2/20 
ELE, ADD, 16, 17/1 //5, 3/21 
NOD, ADD, 22, //ISO 
NOD, ADD, 23, 120// 180 
NOD , ADD , 24 , 240/ / 1 80 
NOD , ADD , 25 , / 1 80/ 1 80 
NOD, ADD, 26, 120/180/180 
NOD , ADD , 27 , 240/ 180/1 80 
ELE, ADD, 17, 12/1/2/7,1/22 
ELE, ADD, 18, 12/1/2/7,3/24 
ELE, ADD, 19, 12/1/2/7,7/25 
ELE, ADD, 20, 12/1/2/7,9/27 
ELE, ADD, 21 , 12/1/2/8,22/23 
ELE , ADD , 22 , 1 2/ 1 /2/8 , 23/24 
ELE , ADD ,23,12/1 /2/8 , 25/26 
ELE, ADD, 24, 12/1/2/8,26/27 
ELE, ADD, 25, 12/1/2/8,22/25 
ELE, ADD, 26, 12/1/2/8,23/26 
ELE, ADD, 27, 12/1/2/8,24/27 
ELE , ADD , 28 , 25/ 1 / / 1 0 , 22 
ELE , ADD , 29 , 25/ 1 / / 1 1 , 23 
ELE , ADO , 30 , 25/ 1 / / 1 O , 24 
ELE , ADD , 31 , 25/ 1 / / 1 0 , 25 
ELE , ADD , 32 , 25/ 1 / / 1 1 , 26 
ELE , ADD , 33 , 25/ 1 / / 1 0 , 27 
PRO, ADD, 1 .36/36/36/36 


MAT, ADD,  l,3.6t6/0.  1 /////////!'. 

PRD,  ADI 

PRO , ADD , 3 , 1 . 1 7E6 

PRO ,  ADD ,  4  -  ? .  34l£6 

PR0,ADD,5,2.5E6 

F'R0,ADD,6,5.0E6 

PRO, ADD, 7, 35. 1/253/2190 

PRO , ADD ,8,67.6/940/1 5000 

PRO , ADD , 9 , 50 . a / 598 / 8200 

M AT ,  ADD ,  2 , 2 .  9f:; 7  / 0 . 3 /  /  /  /  /  /  /  /  7 . 332F,- 4 

PRO, ADD, 10,25.08 

PRO.ADD, 11 ,51.76 

DISP,ADD,1T3, 1 ,0///0/0/0 

D1BP,ADD,4T9, 1 ,0/0//0/0/0 

DISP,ADD, 10T21 , 1 ,0/0/0/0/0/0 

D1 3P , ADD , 22T27 , 1 . / / /O/O/O 

FOR , ADD ,1,1,/ /-306 . 4 

DEC 

XN3 

ANAL-3T A  1  I C 
FILE^==GTSTAT 
SAVE=26 , 27 

♦TITLE, ANALYSIS  OF  CE6193  GTSTRUDL  EXAMPL 
♦LDCASE , 1 

♦EIGCNTL,! ,0,20, 0,0. 0,0. 0,1. OE-5 

♦MQDEQ, 1,20, 1,1 

*EIG0,1 ,1,0, 3, 1,0,1 

QUI 

QUI 

WRI , NISA , GTSTAT. NIS 
WRI ,NISA,QTEIGE.NIS 
WRI , DBS, GTSTAT. BIN 
WRI, DBS, GTEIGE.BIN 
END 
N 
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The  two  previous  pages  contain  the  coaaands  reqaired  to 
generate  this  nodel  using  the  DISPLAY  II  prograa.  With 
the  exception  of  the  following  clarifications/explanations 
all  coaaands  are  explained  in  Chapter  Two. 

COMMAND>VEW , ABS .-80/0/10 

This  command  tells  the  program  to  change  the  angle 
the  user  is  viewing  the  nodel  with.  This  actual  connand 
will  rotate  the  nodel  clockwise  80  degrees  about  the  X 
axis  and  10  degrees  counterclockwise  about  the  Z  axis 
(counterclockwise  angles  are  positive).  This  "new" 
viewing  angle  will  now  be  used  for  all  future  nodel 
"plots" . 

C0NMAN0>PL0 

This  connand  tells  the  progran  to  erase  the  current 
screen  picture  and  to  replot  the  nodel.  The  previous 
connand  changed  the  viewing  angle  and  this  new  "plot"  will 
be  with  the  new  viewing  angle. 

4 

C0MMAND>F0R. ADD, 1 ,l.//-386.4 

The  explanations  for  this  connand  provided  in  Chapter 
Two  Indicate  that  this  connand  tells  the  progran  that 
there  is  a  force  on  386.4  units  in  the  -Z  direction  on 
node  #1.  However  upon  exaninlng  the  figure  on  page  3.4  it 
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Is  noted  that  node  *1  has  no  such  force.  Therefore  why  ii 
this  coaaand  added?  It  was  discovered  while  developing 


this  aodel  that  If  a  static  analysis  is  to  be 


considering  lust  the  aass  of  the  aodel  then  the  coaaand 


"BODYFORCE"  aust  be  used.  However  the  DISPLAY  II  prograa 


does  not  accept  this  coaaand.  As  a  result  the  author  used 
the  coaaand  "FOR , ADD , . . . ”  for  two  purposes: 

1.  As  a  reminder  to  add  the  "BODYFORCE"  coaaand 


to  the  ".NIS"  file  created  by  the  DISPLAY  II 


prograa . 

2.  As  an  easy  aethod  of  deteraining  the  correct 
location  of  Inserting  the  "BODYFORCE"  coaaand  to 
the  ".NIS"  file.  This  is  accoaplished  by 
changing  the  "FORCE"  data  lines  in  the  ".NIS" 
file  to  "BODYFORCE"  data  lines  using  a  line  or 
text  editor.  Refer  to  page  7.3-11  of  the  NISA 
II  User's  Manual  for  the  correct  "BODYFORCE" 
f  oraat . 


CONNAND>*MODEO ,1,20,1.1 

This  coaaand  and  the  "*EIGO"  coaaand  were  both 
discussed  in  Chapter  Two.  However  this  is  a  realnder  that 
the  DISPLAY  II  prograa  does  not  accept  these  coaaands  and 
that  they  aust  be  inserted  into  the  ".NIS"  file  using  a 
line  or  text  editor. 
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COMNAND>WRI , . . . 


All  the  examples  discussed  in  Chapter  Two  had  one 
"NRl , NISA, . . . "  and  one  " WRI , DBS , . . . "  command.  This  model 
has  two  of  each.  This  may  not  be  required;  however,  since 
the  .NIS  file  for  a  static  analysis  contains  a  few 
different  commands  than  the  .NIS  file  for  an  eigenvalue 
analysis,  the  author  felt  that  the  easiest  way  to  perform 
both  analyses  was  to  let  the  DISPLAY  II  program  create  two 
.NIS  files,  one  for  the  static  analysis  and  one  for  the 
eigenvalue  analysis.  Once  the  two  files  are  created,  each 
nay  then  be  easily  tailored  using  a  text  or  line  editor, 
for  the  specific  analysis. 


3.9 


ANAL=FREDUENCY 

DAMP=VISCOUS 

FILE=BTEIGE 

GENF=OFF 

♦TITLE 

FREDUENCY  RESPONSE  OF  CE6193  GTSTRUDL  EXAMPLE 
♦MODE 
1,1,1, 1001 
♦DAMP 
1  C)0 1  ,  .  1 
♦ADDF 
41  .BBS 
♦  SPEC 

AMPLITUDE  SPECTRUM 
1 00 1 , 0 , 3 , 1 . 0 ,1.0 
AMPLITUDE  SPECTRUM 
1 002 , 0 , 3 , 1 , 1 
PHASE  SPECTRUM 
1 003 , 0 , 3 , 1 . 0 , 0 
PHASE  SPECTRUM 
1004,0,3, 1,1.571 
♦DCFQ 

23,FX , 1001 , 1003, «, 30000.0 
26, FX, 1001 , 1003,^,30000.0 
23, FZ, 1002,1004,^,15000.0 
26,F2, 1002, 1004,*, 15000,0 
♦RSET 

1001 .1 , U2,9,1 

1001 .22, UZ,27,1 

1002. 1, UZ, 9,1 

1002. 22,  UZ, 27,1 
♦SPOUT 

1001 ,0,0,1001 
1002,0,0,1002 
♦ENDDATA 


FREQUENCY  RESPONSE  INPUT  EXPLANATION 


All  the  coaaands  listed  in  the  frequency  response  input 
have  been  explained  in  Chapter  Two.  Accordingly  a  line  by 
line  description  is  not  provided.  However  the  following 
general  description  is  provided. 

This  Model  consists  of  a  rigid  frame,  on  a  concrete 
slab,  on  springs.  There  are  two  forces  being  applied  to 
nodes  23  and  26,  one  in  the  X  direction  and  one  in  the  Z 
direction.  The  Z  direction  force  is  90  degrees  out  of 
phase  with  the  X  direction  force.  An  analysis  of  the 
displaceaent  is  the  Z  direction  for  nodes  1  to  9  and  22  to 


27  is  required. 


CIRCULAR  FOUNDATION  ON  A  HALFSPACE 


The  folloMing  exaaple  consists  of  a  circular  foundation  on 
a  halfspace.  A  halfspace  has  no  boundaries  along  three 
edges,  typically  along  the  X  directions  and  one  of  either 
the  Y  or  Z  directions.  Since  an  infinite  length  cannot 
really  be  modeled,  a  halfspace  is  modeled  by  making  the 
halfspace  boundaries  many  times  larger  than  the  diameter 
of  the  foundation  on  the  halfspace.  In  this  model  the 
foundation  has  a  diameter  of  18  feet  and  the  boundaries 
are  located  234  feet  from  the  foundation's  center. 
Additionally  the  use  of  symmetry  is  used  in  this  model.  As 
a  result  only  one  half  of  the  model  is  analyzed.  This 
example  mas  generated  using  DISPLAY  II  and  then  analyzed 
for  a  frequency  response  using  the  NISA  II  program.  The 
NISA  II  analysis  mas  then  compared  to  theoretical  analyses 
of  the  same  model.  This  analysis  is  presented  in  the 
following  order:  Model  sketch,  DISPLAY  II  input,  input 
explanation,  frequency  response  input,  input  explanation. 

Appendix  E  contains  the  analyses  comparisons,  hand 
calculations,  graphs,  and  computer  printouts  used  for  the 
theoretical  analyse^. 
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C5RD,AUD,  1  ,0/0/0 
GRD,TRS, 1 ,2,54 
GRD,TRS, 1,3,234 
LIN, 2QD, 1,2,1 
LIN,2BD,2,3,2 
LIN,TRS, 1 ,3,/180 
LIN,2GD, 1 ,5,/234 
LIN, DEL  2GD 
LIN, DEL, 4 
LIN,TRS,l,5,/234 
LlN,TRS,2,4,/ieO 
LIN,TRS,2,6,/234 
GRD,ADD, 10,9/234 
LIN,2GD,6, 10,7 
LIN,TRS,7,B,/6 
PAT , 2LN ,1,3,1 
PAT , 2LN ,2,4,2 
PAT , 2LN , 5 , S , 3 
PAT,2LN,4,6,4 
PAT,2LN,7,S,5 

FEB , QUA , 1 , , , 2/B/3/3 , 1 / 1 /B/0 . 3333 
FEG,QUA,2, , ,2/3/3/3, 1/1/3/0.3333 
FEG,QUA,3, , ,2/B/S/3, 1/1/8/0. 12S 
FEG,QUA,4, , ,2/3/8/3,l/l/3/0. 12S 
FEG , QUA , 5 , , , 2/3/2/3 , 2/2/3/ 1 
N0D,HER,1T484 
N 

MAT,ADD,l,6.0E5/0.4S////////3.88 

MAT,ADD,2,5.76EQ/0.45////////4.66 

D ISP, ADD, 478/474/467/463/344/335/318/309, 1 ,0 

DI8P, ADD, 292/283/266/257/240/231/214/205, 1 ,0 

D I SP , ADD , 1 88/ 1 79/ 162/1 53/79/70/53/44 ,1,0 

DISP, ADD, 27/18/106/1 13/1 17/124/128/135, 1 ,0 

DISP, ADD, 371/378/382/389/393/400/404/41 1,1,0 

DISP, ADD, 415/422/426/433/437/444/448/455, 1,0 

DISP,ADD, 1T17/97T102, 1 ,0 

DEC 

INS 

ANAL»EIGENVALUE 

FILE»S0LID 

SAVE=26,27 

‘TITLE, CIRCULAR  FDN  ON  A  HALFSPACE 

*EIGCN,40,0,40,0,C).0,0.0,1.0E-5 

*LDCASE, 1 

QUI 

QUI 

WRI ,NISA, SOLID. NIS 
WRI , DBS, SOLID. BIN 
END 
N 


3.14 


The  previous  page  contains  the  coaaands  required  to 
generate  this  aodel  using  the  DISPLAY  II  prograa.  With 
the  exception  of  the  following  clarifications/explanations 
all  coaaands  are  explained  in  Chapter  Two. 

CONNAND>NOD.MER. 1T484 

This  coaaand  tells  the  prograa  to  review  all  nodes 
between  1  and  484  (all  the  nodes  of  the  aodel)  and  to 
coablne  any  two  or  aore  nodes  together,  into  one,  that  are 
located  within  the  "default"  distance  of  each  other.  This 
is  accoapllshed  with  the  following  foraat:  "NOD , HER , 3 , 4 " 
where : 

Entry  3  is  the  list  of  nodes  to  be  aerged. 

Entry  4  is  the  "tolerance",  ie.  the  specified 
distance  between  nodes. 


If  a  aodel  was  generated  with  aore  than  one  patch  this 
coaaand  aust  be  used  to  "tie"  the  patches  together.  I f 


this  coaaand  is  not  used  on  a  aodel  created  by  aore  than 


one  patch,  it  cannot  be  analyzed.  Refer  to  pages  4.3-12  to 
4.3-14  of  the  DISPLAY  II  User's  Manual  for  aore 


inf oraation . 


COMMAND>N 

If  no  tolerance  is  specified  in  the  above  coaaand  the 
prograa  will  ask  if  the  user  would  like  to  aodify  the 
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default  tolerance  (Y/N)?  The  author  chose  to  use  the 


default  tolerance  for  this  merge  and  therefore  responded 
with  "N" . 

The  element  type  used  in  this  model  and  analysis  was  NKTP 
3,  an  axisymmetric  solid  element.  This  was  the  first 
analysis  the  author  performed  using  this  type  of  element 
as  well  as  the  first  model  generated  using  multiple 
patches.  While  performing  this  model  generation  and 
analysis  the  following  restrictions,  not  found  in  any  of 
the  NISA  II  User's  Manuals,  were  discovered; 

The  NISA  II  User's  Manual  states  that  "The  element 
must  lie  in  the  right  half  of  the  global  XY  plane  where 
the  global  X-axis  is  the  radial  direction,  and  the  global 
Y-axis  is  the  axis  of  revolution  or  axis  of  symmetry." 

The  author  determined  that  not  only  does  this  element  have 
to  remain  in  the  right  half  of  the  global  XY  plane  but  it 
must  also  remain  in  the  first  quadrant  of  the  global 
coordinate  system,  le  all  nodal  coordinates  must  be 
positive . 

d 

All  the  element  meshes  must  be  generated  starting 
with  the  patch  closest  to  the  origin  of  the  coordinate 
system,  move  left  along  the  X-axis,  then  move  up  to  the 
next  layer  of  patches  and  move  left  along  the  X-axis  from 
the  Y-axis  again.  The  following  figure  contains  5 
patches,  blocks  1  to  5.  To  correctly  generate  a  finite 
element  mesh  from  these  patches,  start  with  patch  1  and 
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finish  with  patch  5,  coapletlng  the  generation  in 
numerical  order. 


V  I 


ANAL=FREQUtNCY 

DAMP=VISCOLJB 

FILE=SQLID 

GENF»OFF 

♦TITLE 

SINUSOIDAL  LOADING  ON  A  CIRCULAR  FDN  ON  A  HALFSPACE 
♦MODE 

1,39,1,1001 
♦DAMP 
1001,0 
♦ADDF 
4. 0,8.0, 1 
♦SPEC 

AMPL I  Tl JDE  SPECTRUM 
1001,0,3,1.0,1. 0 
PHASE  SPECTRUM 
1002,0,3,1.0,0.0 
♦DCFO 

473 , FY , 1 00 1 , 1 002 ,478,1. 30000 
♦RSET 

1001 .344, UY,3S0, 1 

1002. 344,  UX, 360,1 
♦SPOUT 

1001 ,0,0,1001 
1002,0,0,1002 
♦ENDDATA 
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FREQUENCY  RESPONSE  INPUT  OP  A  CIRCULAR  FOUNDATION  ON  A 

HALFSPACE 

All  the  coBMands  listed  In  the  frequency  response  input, 
shown  on  the  previous  page,  have  been  explained  in  Chapter 
Two.  Accordingly  a  line  by  line  description  is  not 
provided.  However  the  following  general  description  is 
provided . 

This  model  consists  of  an  18  foot  diaaeter  circular 
foundation  on  a  halfspace.  The  halfspace  is  aodeled  by 
setting  a  displaceaent  constraint  in  the  X  and  Y 
direction  234  feet  below  the  base  of  the  foundation  and 
setting  a  displacement  constraint  in  the  X  direction  along 
the  center  of  the  foundation  and  234  feet  to  the  right  of 
the  foundation's  center.  There  is  one  force  being  applied 
to  the  center  of  the  foundation,  node  478,  in  the  Y 
direction.  An  analysis  of  the  displacements  in  the  X  and 
Y  directions  for  nodes  244  to  360  is  required.  These 
nodes  are  located  on  the  surface  (top)  of  the  halfspace 
starting  at  the  center  of  the  foundation  and  ending  54 
feet  away. 


CHAPTER  FOUR 


CONCLUSIONS 


The  following  conclusions  can  be  wade  based  on  the  results 
of  this  research  problem: 

1.  There  are  many  different  Finite  Element  Analysis 
programs  available  on  the  market  today.  One  of  the 
programs  available  for  use  on  a  personal  computer  is 
called  NISA  II  and  is  produced  by  Engineering  Mechanics 
Research  Corporation  located  in  Troy  Michigan.  NISA  II  is 
a  very  powerful  program  that  is  relatively  user  friendly, 
once  the  user  becomes  familiar  with  it.  provides  rapid 
analysis,  and  is  fairly  inexpensive  to  use  on  a  personal 
computer.  Additionally  it  provides  answers  in  good 
agreement  with  check  cases. 

2.  Even  though  a  program  is  considered  user  friendly,  as 
in  the  case  of  the  NISA  II  program,  it  still  requires  an 
extensive  amount  of  "user  time"  for  someone  who  is 
unfamiliar  with  the  program  to  become  proficient  in  its 
utilization .  A  large  amount  of  this  user  time  is 
typically  spent  trying  to  solve/understand  roadblocks 
encountered  that  are  not  clearly  explained  in  the  program 
provided  user  manuals.  Once  an  individual  becomes 
proficient  with  a  program,  "reasonably"  correct  answers 
should  be  provided. 

3.  The  NISA  II  User's  Manuals  provide  a  broad 
explanation  on  how  to  use  the  program.  This  is  fine  for 
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soaeone  who  is  familiar  with  the  basics  of  a  finite 
element  modeling  program.  This  report,  on  the  other  hand, 
is  designed  to  provide  a  much  more  basic  user's  manual  on 
how  to  use  the  NISA  II  computer  program.  As  a  result  this 
report  gives  complete  step  by  step  examples  on  how  to 
develop  and  analyze  various  types  of  static  and  dynamic 
problems,  as  well  as  Invaluable  hints  on  where  to  avoid 
problems  not  discussed  in  the  NISA  II  User's  Manuals. 

4.  The  comparison  between  NISA  and  other  published 
values,  of  an  axisymmetric  solid  on  a  halfspace  under  a 
sinusoidal  loading,  show  the  following: 

A.  All  values,  even  those  of  the  published 
solutions,  tend  to  have  rather  large  differences. 

B.  NISA  II  values  tend  to  compare  reasonable  well, 
not  varying  any  greater  than  any  of  the  other 
solutions . 

C.  The  NISA  II  values  were  for  a  "model"  with  rigid 
boundaries.  These  boundaries  were  located  so  far 
away  from  the  loaded  solid,  that  it  was  presumed  they 
would  have  no  effect.  However  as  no  other  finite 
element  solution  was  obtained,  for  a  comparison,  this 
can  not  be  completely  determined. 

5.  When  comparing  the  NISA  II  solutions  to  those 
calculated  using  Lysaer  it  was  noted  that  Lysmer  provided 
displacements  in  excess  of  six  times  the  amount  provided 


by  NISA  II.  The  following  are  possible  reasons: 

A.  The  NISA  II  Model  had  rigid  boundaries,  which  aay 
have  allowed  energy  to  bounce  off,  where  as  the  Lys«er 
finite  element  model  had  energy  adsorbing  boundaries. 

B.  The  Lysmer  curve  of  frequency  verses  displacement 
had  only  one  peak.  The  NISA  II  curve  of  frequency  verses 
displacement  gave  multiple  peaks.  The  reason  for  the 
multiple  peaks  could  be  due  to  the  rigid  boundaries  and 
therefore,  this  needs  to  be  looked  into  more. 

C.  As  a  result  of  the  above  it  is  highly  recommended 
that  an  exact  comparison  between  NISA  and  another  Finite 
Element  program  be  accomplished.  Based  on  discussions 
between  the  author  and  the  author's  advisor  it  is 
understood  that  this  comparison  will  be  done  in  the  near 
future . 


I 
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REQUIRED  EQUIPMENT 


In  order  to  use  the  NISA  II  prograa.  version  88.7,  on  a  PC 
the  following  ainiaua  equipaent  is  required: 

1.  IBM  PC-XT  or  IBM  PC-AT  coapatible  Personal 
Coaputer.  The  author  used  a  HP  386. 

2.  640K  bytes  of  central  aeaory. 

3.  IBM  coapatible  Color  Graphics  Adaptor  (CGA)  and 
aonitor . 

4.  A  20  MB  hard  disk  with  a  360  KB  floppy  drive. 

5.  A  Math  Coprocessor  (Intel  8087  of  80287  chip). 

6.  A  132  character  line  printer. 

7.  MS-DOS  Version  2.1  of  above  operating  systea. 

8.  A  Parallel  port  for  the  protection  device 
supplied  by  EMRC . 

Additionally  the  following  equipaent  is  Optional: 

1.  IBM  coapatible  Enhanced  Graphics  Adaptor  (EGA) 
and  Enhanced  Display  aonitor. 

2.  IBM  compatible  Professional  Graphics  Controller 
(PGA)  and  Professional  Display  aonitor. 

3.  Ink  Jet  printers  (ie  Tektronix  4696)  and  screen 
duap  software  for  EGA  and  PGA. 

The  author  encountered  a  problea  during  the  initial  set  up 


A.  1 


and  operation  of  the  proeraa.  This  problea  was  due  to 
the  MS-DOS  version  (DOS  3.3)  and  all  the  associated 
software  requiring  too  auch  space.  As  a  result  the 
hardware  did  not  have  600K  bytes  of  RAN  available  for  the 
NISA  II  and  DISPLAY  II  prograas.  To  over  coae  this 
problea  the  author  deleted  everything  except  the  aain  DOS 
3.3  froa  the  operating  systeas  start  up  path.  This 
resulted  in  having  about  596K  of  RAM  available  for  the 
NISA  II  and  DISPLAY  II  prograas  to  operate,  which  was 
enough . 
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APPENDIX  B 
RUNNING  TIME 


The  author  used  a  HP  VECTRA  QS/16  to  perform  all  analysis 
on.  The  actual  time  required  to  perform  an  analysis 
depends  more  on  the  type  of  problem  being  analyzed  than 
anything  else.  However  in  an  attempt  to  provide  some  type 
of  guide  to  running  time  the  following  graphs  are 
provided.  As  shown  on  these  graphs,  whether  it  is  the  log 
total  time  or  Just  cpu  time,  the  number  of  degrees  of 
freedom  is  the  controlling  factor  on  length  of  time 
required  to  perform  an  analysis. 

In  addition  to  the  above  "guide”,  the  following  is 
provided.  Two  eigenvalue  analyses  were  performed  on  a 
model  with:  860  DOF.  430  nodes,  127  elements,  and  a  wave 
front  of  84.  One  analysis  used  "accelerated  subspace" 
eigenvalue  extraction  and  the  other  analysis  used 
"Inverse"  eigenvalue  extraction.  The  accelerated 
subspace  method  took  approximately  3993  seconds  (about  67 
minutes)  to  provide  the  first  40  eigenvalues.  The  inverse 
method  took  approximately  7708  seconds  (about  1  1/2  hours) 
to  extract  30  eigenvalues  starting  with  a  frequency  of  2 
cps . 
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NISA  II  EXECUTIVE  COMMANDS 


Input  Dtu  Setup 
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Table  52Ji  Alphabetical  list  of  all  executive  commands  for  all  analyses 


Executive  Gotnmand  Name^^^ 

Description  / 

applicable  Analyiis  Typas^^ 

ACUTOFF 

Range  of  model  coupling 

RA 

ANALYSIS 

Specify  analysis  type 

ALLl,  ALU 

AUTC  CONSTRAINT 

Spurious  noimsl  totarion  conitraint 

ST.EV.BU 

BLANK  COMMON 

Blank  common  storage  unit 

ALLl.  ALU 

DAMPING 

Type  of  damping 

ALU 

DELTTTME 

Tune  increment 

TR 

DIRECTION 

factional  superposition  procedure 

SHO 

EIGEN  EXTRACTION 

Eigenvalue  extraction  method 

EV.BU 

ELEMENT  ECHO 

Echo  element  input 

ALLl 

EXECUTION 

Select  execution/checking  run 

ALLl.  ALU 

ENDTIME 

1  Ending  time 

I  TR 

rlename  I 

Specify  NISA  file  prefix  | 

ALLl.  ALU 

FLOWER 

Lower  Hfflit  of  exciting  fteqnency 

RA.FR 

FRQRDP 

firtnr  fnr  awwyfiwy 

RA.FR 

FSMALL 

Reqioency  cut-off  to  identify  ligid 
body  modes 

ALU 

FUPPER 

Upper  Umit  of  exciting  fteqnency 

RA,PA 

GENFREQUENCY 

fiequacy  points 

RA.FA 

GFACTOR 

Acoeleiaticn  dtre  to  gravity 

AlU 

INITUL  TEMPERATURE 

Define  initiai  tenqsenture 

TH 

INPOLATION 

Type  of  inteipoiatian 

RA.FR.SHO 

INPHASE 

Type  of  interpolation  for  phase  spectra 

FR 

INTEGRATION 

Type  of  integratian  procedure 

RA 

MASS  FORMULATION 

Mass  fotmulsiion  for  eigenvalue 
analysis 

EV 

(1)  Miniimim  atbieviaiioas  are  in  bold  face 

(2)  ST:  Linear  static  ‘  EV:  Eigenvalue  BU:  Bockliaf 

NL:  Nonlinear  stadc  SH:  Steady  state  heat  transfer  TH:  Tnnsieat  heat  transfer 

TR:  Transient  dynamics  FR:  Frequency  response  RA:  Random  vibniioo 

SHO:  Shock  spectnnn 

ALLl:  All  static  and  beat  transfer  analyses  (ST.  EV.  BU.  NL.  SH.  TH) 

ALL2:  All  modal  dynamic  analyses  (TR.  FIL  RA.  SHO) 
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Table  5JL3  Alphabetical  list  of  all  executive  conunands  for  all  analyses  (continued) 


Executive  Command  Name^^^ 

Deaciipcion  4 

\pplicahle  Analysis  Types^ 

MAXCPUTIMB 

Limit  of  CPU  aecooda 

NL 

MRESPONSE 

Control  printout  of  modal  icqxmaes 

FR 

NLTYPB 

Spedfy  noolinear  static  analysis  type 

NL 

NODE ECHO 

Edto  node  Input 

ALLl 

ORTHOTROPIC 

Orthociopic  direction  definidon 

SH.TH 

RESEQUENCE 

Element  resequencing 

ALLl 

restart 

Restart  options 

SH.TH 

SAVEFIUE 

Save  spediic  NISA  files 

ALLl.  ALL2 

SIGMA 

Stefan-Boltzmann  coefficient 

SH.TH 

SORT  STRESS 

Son  highest  stresses 

ST 

STARTINOTIME 

Starting  time 

TR 

STEP 

Total  number  of  steps  for  heat  transfer 

analysis 

TH 

STRESS  SnnFFENINQ 

Stress  optiow 

ST 

WARNDIO 

SetwaolBgflig 

AU.1.AXX2 

(1^  rthwvirtoM  we  hi  bolld  ftce 

(2)  ST^  liawMde  EV:  EifenvalBe  BU:  BackUnf 

NL;  NwiBi— ’Wb  SH:  SiaadyMMbeanMftr  1H:  DaMlMbiK 

SHO.  SbockvecttHi 

ALLl:  AnMic«dlMCttMte«alyaes(ST.EV.BU.NL.SH.1H) 

ALL2:  AH  nodal  dynaaik  analyses  (IB,  FR,RA,  SHO) 
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APPENDIX  E 

THEORETICAL  DISPLACEMENT  CALCULATIONS 
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*  ♦ 

«  DYNAMIC  ANALYSIS  * 

*  using  '  GAZETAS  '  METHOD  * 

*  Earthquake  engineering  voi.L'  nol,lVS3  * 

*  * 

*  * 

*  IBM  version  l.O  * 

*  * 

*  ♦ 


^ »  2  CP^ 

*««■»«■««««««««*««*«■«««*■» 

*  Input-data  ♦ 

♦■it-**********#**********.********************* 


MASS  DENSITY  OF  SOIL  (KIPS-SEC-'-2/FT-'-4>  =  .  0038B 

DEPTH  TQ  BEDROCK  CH]  (FT)  =  234 

FOUNDATION  EMBEDMENT  CD 3  (FT)  »  0 

SHEAR  MODULUS  OF  SOIL  (KSF)  =  206 


Data  for  little  Ckl  and  little  Cc3  i 


HORIZONTAL  DYNAMIC  STIFFNESS  CKH3  -  1 

VERTICAL  DYNAMIC  STIFFNESS  CKV3  -  .97 

ROCKING  DYNAMIC  STIFFNESS  CKR]  •  .94 

TORSIONAL  DYNAMIC  STIFFNESS  CKT3  »  .96 

HORIZONTAL-ROCKING  STIFFNESS  CKHR]  »  1 

VERTICAL  DAMPING  COEFFICIENT  CCV3  -  .84 

HORIZONTAL  DAMPING  COEFFICIENT  CCH3  «  .59 

ROCKING  DAMPING  COEFFICIENT  CCR]  =  .11 

HORIZONTAL-ROCKING  DAMPING  COEFFICIENT  LCHR3  »  0 

TORSIONAL  DAMPING  COEFFICIENT  CCT]  »  9.000001E-02 

POISSON  RATIO  =  .45 

ROCKING  MASS  OF  INERTIA  -  CG  /AXIS  X  (KIP-SEC'"2-FT)  -  119.631 
TORSION  MASS  OF  INERTIA  -  CG  /AXIS  Z  (KIP-SEC^2-FT)  »  208.388 
MASS  OF  FOUNDATION  ♦  MACHINERY  (K1P-SEC''2/FT)  »  5.145 


E.8 


TORSIONAL  b/.CITlNG  t-IOHLNT  K  l-Mf)  =  U 
HORIZONTAL  EXCITING  FORCE  (LIPS)  =  0 

ROCKING  EXCITING  MOMENT  (FT-KIP)  =  0 

VERTICAL  EXCITING  FORCE  (KIPS)  »  50 

EQUIVALENT  RADIUS  OF  FOUNDATION  (FT)  -  9 

OPERATING  FREQUENCY  OF  SYSTEM  (RPM)  »  120 

VERTICAL  PHASE  ANGLE  IN  DEGREES  -  0 

HORIZONTAL  PHASE  ANGLE  IN  DEGREES  =  0 

ROCKING  PHASE  ANGLE  IN  DEGREES  =  0 

TORSIONAL  PHASE  ANGLE  IN  DEGREES  =  0 

HYSTERETIC  MATERIAL  DAMPING  (7.)=  0 

DISTANCE  FROM  CENTER  OF  GRAVITY  TO  BASE  OF  FOUNDATION  (FT) 

*  OUTPUT-DATA  ♦ 

HORIZONTAL  STATIC  STIFFNESS  (KIP/FT)-  9753.0S2 
VERTICAL  STATIC  STIFFNESS  (KIP/FT) =  26172.78 

HORIZONTAL-ROCKING  STATIC  STIFFNESS  <KIP/FT)«  0 
TORSIONAL  STATIC  STIFFNESS  <KIP/FT)»  800926.1 
ROCKING  STATIC  STIFFNESS  (KIP/FT)  *=  7327BG.8 

VERTICAL  DISPLACEMENT  (FT)»  1.862894E-03 

VERTICAL  PHASE  ANGLE  (RAD)«  -.4137561 
TORSIONAL  DISPLACEMENT  <FT)«=  O 
TORSIONAL  PHASE  ANGLE  (RAD) a  -4. B03612E-02 
HORIZONTAL  DISPLACEMENT  (FT)»  O 
HORIZONTAL  PHASE  ANGLE  (RAD)«  1.20483 
ROTATIONAL  DISPLACEMENT  (FT)-  0 

ROTATIONAL  PHASE  ANGLE  (RAD)*  1.20483 
*****  End  a-f  Dynamic  Analyai*  *♦*•* 
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R(3)  - 

9 

R(5)  - 

9 

R<7)  =« 

9 

R(9)  - 

9 

AI  (1) 

-  .2895918 

AI  (3) 

-  .4123002 

AI  (5) 

-  5.399169E-02 

AI  (7) 

-  0 

AI  (9) 

-  4.417502E-02 

KI  (1) 

-  9753.052 

KI  (2) 

-  2824.404 

KI  (3) 

-  25387.6 

KI  (4) 

»  10791.04 

KI  (5) 

»  688816.8 

KI  <6) 

=  39564.24 

KI  (7) 

=  0 

KI  (8) 

=  0 

KI  (9) 

=  768891 

KI (10) 

=  3538 1 . 02 

HR=  8940.587  HI=  282 

64933.67 
NR=  5.e06465E+09 


CR«=-29259.  16  CI=-8473.211  RR= 


10791.04 
PV»-. 4137561 
05-  0 

RA-  760484.4 
9376. 106 
HI-  0 

H6-  0 
H7«  0 
II-  0 


V3=  26839.96 
01-  735983.7  02' 

PO— 4.803612E-02 
RB-  8.555496E-02 
HB-  .3059867 
H2-  8.555496E-02 


NI=  2.22521SE+09 
V4-  .4137561  V5= 

35381.02  03= 


Vl  = 

1 . 862B94E-03 
736833.6  04= 


CA-  30461.35  CB* 


H3-  0 


H4i 


H8-  0  H9-  O 

12-  0  13-  0 

16-  .3659669 

Rl«  0  R2-  .3059867  R3-  0 

R7-  0  RB-  0  R9-  0 

SI-  0  82-  0  83-  0 

PR-  1.20463 

End  of  internal  storage  data 


PH-  1.20463 
14-  1.570796  I5« 

R4-  3.423471  R5* 

84-  1.570796 


757702.9  RI= 
24575. 13  V2= 

4.803612E-02 
3.423471  HA- 
3. 423471  H5-  0 

6.21B248E+09 
0  R6-  0 
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DYNAMIC  ANALYSIS 
using  *  GAZETAS  '  METHOD 
Earthquake  engineering  vol.2  nol.1983 


IBM  version  1 .0 


f »  3  Cp4  MAS^ 


Input  -  data 


* 


MASS  DENSITY  OP  SOIL  CKIPS-SEC“2/FT“4 ) -  .00388 

DEPTH  TO  BEDROCK  [H]  (?T)  -  234 

FOUNDATION  EMBEDMENT  [D]  (FT)  -  0 

SHEAR  MODULUS  OP  SOIL  (KSF)  •  206 

Data  for  little  [k]  and  little  (e] 

HORIZONTAL  DYNAMIC  STIFFNESS  [KH]  -  .99 

VERTICAL  DYNAMIC  STIFFNESS  CKV] 

ROCKINO  DYNAMIC  STIFFNESS  [XN]  -  .9  ‘ 

TORSIONAL  DYNAMIC  STIFFNESS  [KT]  •  .92 

HORIZONTAL-ROCKING  STIFFNESS  [KBR]  -  i 
VERTICAL  DAMPING  COEFFICIENT  (CV]  •  .85 

HORIZONTAL  DAMPING  COEFFICIENT  [CH]  -  .6 

ROCKING  DAMPING  COEFFICIENT  [CR]  •  .1 

HORIZONTAL-ROCKING  DAMPING  COEFFICIENT  [CHR]  •  0 

TORSIONAL  DAMPING  COEFFICIENT  [CT]  -  9.000001E-02 


POISSON  RATIO  »  .45 


MASS  OF  FOUNDATION  «  MACHINERY  (KIP-SEC'2/FT)  -  5.145 

TORSIONAL  EXCITING  MOMENT  (FT'KIP)  -  0 

HORIZONTAL  EXCITING  FORCE  (KIPS)  -  0 

ROCKING  EXCITING  MOMENT  (FT-KIP)  <  0 

VERTICAL  EXCITING  FORCE  (KIPS)  «  50 

EQUIVALENT  RADIUS  OF  FOUNDATION  (FT)  -  9 

OPERATING  FREQUENCY  OF  SYSTEM  (RPM)  -  180 

VERTICAL  PHASE  ANGLE  IN  DEGREES  -  0 

HORIZONTAL  PHASE  ANGLE  IN  DEGREES  -  0 

ROCKING  PHASE  ANGLE  IN  DEGREES  -  0 

TORSIONAL  PHASE  ANGLE  IN  DEGREES  »  0 

KYSTERETIC  MATERIAL  DAMPING  («)-  0 

DISTANCE  FROM  CENTER  OP  GRAVITY  TO  BASE  OF  FOUNDATION  (FT) 


OUTPUT  -  DATA 


HORIZONTAL  STATIC  STIFFNESS  (KXP/FT)- 
VERTICAL  STATIC  STIFFNESS  (KIP/FT)-  tSlTt.ti"* 
HORIZONTAL-ROCKING  STATIC  STIFFNESS  (KIP/FT)-  0 
TORSIONAL  STATIC  STIFFNESS  (KIP/FT)-  SOOSSS.l 
ROCKING  STATIC  STIFFNESS  (KIP/FT)-  79STS3.S 
VERTICAL  DISPLACEMENT  (FT)-  1.837587E-03 

VERTICAL  PHASE  ANGLE  (RAD)-  -.6459611 
TORSIONAL  DISPLACEMENT  ^FT)-  0 

TORSIONAL  PHASE  ANGLE  (RAD)-  -7 . 989973E-02 


HORIZONTAL  DISPLACEMENT  (FT)-  0 


rotational  phase  angle  (RAD) 


.9740001 


End  of  Dynaaic  Analysis  ••••• 


VS-  230. 

4187 

R(l)  - 

9 

R(3)  - 

9 

R(5)  - 

0 

R(7)  - 

9 

R(9)  - 

9 

AI(l)  - 

.4417802 

AI(3)  - 

.6288128 

AI(5)  - 

7.362503E-02 

AI(7)  - 

0 

AI(9)  - 

6.626253E-02 

KI(1)  - 

9655.521 

KI(2)  - 

4308.412 

KI(3)  - 

23555.5 

KI(4)  - 

16379.26 

KI(5)  - 

659505.4 

KI(6)  - 

53951.23 

KI(7)  - 

0 

KI(8)  - 

0 

KI(9)  - 

736853.8 

KI(IO) 

-  53071.52 

HR-  7627.477  HI-  C30 

92716.94 
II-  4.43925aE<i-00 
16979.26  V9< 

PV— .6489611  01< 

09-  0  PO< 


C306.412  CR--26966.87  Cl  — 12925.24  RI-  703699.6  RI< 


27309.6 

662612.8 


II* 

V4« 

02* 


3.00t703l-09 
.34ft311  V8- 

83071.82  03- 


VI-  21727.48  V2- 

1.6378371-03 

664933.6  04-  7.9399731-02 


RA-  709980.9 
.8031660 

RB- 

.130979 

CA-  31719.43  .  C3-  3.86129 

.11-  0 

12- 

. 130979 

13-0 

■  I4-'' 

3.86130^ 

17-  0 

H8- 

0 

19-  0 

Pl- 

.9749001 

II-  0  12- 

16-  .5958963 

0 

13-  0 

14- 

1.570796 

Rl-  0 

R7-  0 

R2- 

R8- 

.5031669 

0 

R3-  0 
R9-  0 

R4- 

3.56129 

31-  0 

PR-  .9749001 

32- 

0 

33-  0 

34- 

1.570798 

18-  0 


R6-  0 


End 
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DYNAMIC  ANALYSIS 
utlne  '  CAZETAS  '  METHOD 
Earthquake  enginterlnc  vol.2  nol.lOSS 


IBM  version  1.0 


m 

0 

0 

0 


♦'Hi 


Input  -  data 


MASS  DENSITY  OF  SOIL  (KIPS-SEC‘2/FT'4 ) -  .00388 

DEPTH  TO  BEDROCK  [H]  (FT)  -  234 

FOUNDATION  EMBEDMENT  [D]  (FT)  -  0 

SHEAR  MODULUS  OF  SOIL  (KSF)  -  SOS 
Data  for  llttlo  [k]  aad  littlo  [o]  t 


HORIZONTAL  DYNAMIC  STIFFNESS  [NH]  -  .SB 
VIRTICAL  OYMANZC  STIFFilSS  IKYJ  ‘ 

ROCKXRO  OYIANXC  STIFFNESS  [KR]  -  .• 

TORSIONAL  DYNAMIC  STIFFNESS  IRT]  -  .BB 
HORXZOHTAL'ROCKIHQ  STIFFNESS  [KHR)  -  1 

VERTICAL  DAMPINO  COEFFICIENT  [CV]  -  .SB 
HORIZONTAL  DAMPING  COEFFICIENT  [CB]  -  .6 

ROCKING  DAMPING  COEFFICIENT  [CR]  -  .2 

HORIZONTAL-ROCKING  DAMPING  COEFFICIENT  [CHR]  •  0 

TORSIONAL  DAMPING  COEFFICIENT  [CT]  -  B.OOOOOlR-02 


POISSON  RATIO  -  .45 


NASS  OP  FOUNDATION  *  MACHINERY  ( K 1 P-SCC ‘2 /PT )  -  0 

TORSIONAL  EXCITING  MOMENT  (FT-KIP)  -  0 

HORIZONTAL  EXCITING  FORCE  (KIPS)  •  0 

ROCKING  EXCITING  MOMENT  (PT'KIP)  •  0 

VERTICAL  EXCITING  FORCE  (KIPS)  -  50 

EQUIVALENT  RADIUS  OP  FOUNDATION  (PT)  -  9 

OPERATING  FREQUENCY  OF  SYSTEM  (RPN)  •  180 

VERTICAL  PHASE  ANGLE  IN  DEGREES  -  0 

HORIZONTAL  PHASE  ANGLE  IN  DEGREES  -  0 

ROCKING  PHASE  ANGLE  IN  DEGREES  -  0 

TORSIONAL  PHASE  ANGLE  IN  DEGREES  -  0 

HYSTERETIC  MATERIAL  DAMPING  («)•  0 

DISTANCE  FROM  CENTER  OP  GRAVITY  TO  BASE  OP  FOUNDATION  (PT)- 


OUTFUT  >  DATA 


BORIZOMTAL  STATIC  STXFfMBSS  (KlF/fT)« 

VERTICAL  STATIC  STXVVVESS  (KIF/FT)*  ' 

HORIZONTAL-ROCKING  STATIC  STIFFNESS  (KIF/FT)-  0 
TORSIONAL  STATIC  STIFFNESS  (KIF/FT)-  SOOBtS.l 
ROCKING  STATIC  STIFFNESS  (KIF/FT)-  TStTSS.S 
VERTICAL  DISPLACENENT  (FT)-  1.742741E-03 

VERTICAL  PHASE  ANGLE  (RAD)-  -.6075967 
TORSIONAL  DISPLACENENT  (FT)-  0 

TORSIONAL  PHASE  ANGLE  (RAD)-  -7 . 190033E-02 


HORIZ^  hL  DISPLACEMENT  (FT) 


0 


ROTATIONAL  PHASE  ANCLE  (RAD) 


1  .06ft473 


•••••  End  of  Oyntaic  Analytii  ••••• 

VS-  230. 4187 
R(l)  -  ® 

R(3)  -  9 

R(5)  -  9 

R(7)  -  9 

R(0)  -  9 

Al(l)  -  .4417S02 

AI(3)  -  .6258128 

AI(5)  -  7 .362503E-02 

AI(7)  -  0 

AI(9)  -  6.626253E-02 

KI(1)  -  9655.521 

KI(2)  -  4308.412 

KI(3)  -  23555. S 

KI(4)  -  16379.26 

KI(5)  -  639505.4 

KI(6)  -  53951.23 

KI(7)  -  0 

KI(8)  •  0 

KI(0)  •  738853.8 

ZI(IO)  -  53071.52 

HR-  9855.521  HI-  4308.412  CR— 28986.57  Cl  — 12925.24  RR-  746405.1  RI- 

92728.94 

■R-  9.189424t«09  HI*  8.88I849|o09  -VI*  28898.9  V2- 

19879.88  V8*  88880.44  V4«  .8878887  V8-  1 .74874ljr-8V 

PV— .8078987  Ol-  788888.8  08*  88071.88  08-  788788.8  ^  04-  7.188088B-02 

08-  0  PO— 7.180088B«08 

RA-  752142.8  R8-  .1888981  CA-  31719.48  C8-  3.88188  lA-  10878.15  H8- 

81-  0  18-  .  1888981  ""|8-  0  If.  o 

H7-  0  H8-  0  19-  0  PI-  1.88#478^ 

II-  0  12-  0  13-  0  .  14-  1.870798  18-  8.998343E-09 

16-  .5013234 

Rl-  0  R2-  .4198996  R3-  0  R4-  3.58129  R5-  0  R6-  0 

R7-  0  R8-  0  R9-  0 

SI-  0  82-  0  S3-  0  84-  1.870798 

PR-  1.069473 

End  of  I'tornal  storafe  data 


E.23 


